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This  volume  comprises  three  papers  written  In  compliance  with 
the  U.S.  Coast  Guard-CRPLi  contract  D0T-CG-842923’-A.  Part  I:  "Geophysical 
and  Geological  Data  Base  Evaluation  for  Loran-C  Ground  Wave  Propagation 
Medium"  comprises  Task  1  under  this  contract.  Part  II:  "Loran-C  Ground 
Wave  Secondary  Phase  Corrections  Over  Nonhomogeneous  and  Irregular  Ground 
Using  Transient  Signal  Propagation  Techniques"  rapresents  Task  3  under 
this  contract.  Finally,  Part  III:  "An  Experimental  Program  to  Measure 
ECD  and  Test  the  Transient  Solution  (With  Addendum)"  comprises  Task  2 
under  this  contract.  Each  paper  is  designed  separately.  However,  the 
common  topic  of  Loran-C  Pulse  Propagation  engulfs  all  three  papers  and 
are  presented  here  as  the  final  documentation  on  the  above  referenced 
contract. 
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ABSTRACT 


The  geophysical  data  available  for  a  select  propagation  path  for  use 
in  the  study  .of  the  propagation  of  a  Loran-C  pulse  is  evaluated  and  recom¬ 
mendations  for  the  initial  estimates  of  ground  geoelectric  section  resis¬ 
tivities  is  presented.  Soil,  geological  data  and  terrain  are  considered 
in  the  impedance  data.  Coupling  between  ground  surface  impedance  and  ter¬ 
rain  is  postulated  through  the  soil  depth.  The  continuity  of  the  data  is 
stressed. 


GEOPHYSICAL  AND  GEOLOGICAL  DATA  BASE  EVALUATION  FOR  LORAN-C 
GROUND  WAVE  PROPAGATION  MEDIUM 

J,  R,  JOHLER 


1.  .  INTRODUCTION 

Operation  of  Loran-C  overland  introduces  a  need  for  using  the  full 
elegance  of  the  ground  wave  theory  that  has  been  developed  during  the  past  70 
years.  Loran-C  is  a  pulse  transmission  but  the  development  of  special  tech¬ 
niques  for  pulses  permits  the  use  of  this  ground  wave  theory  by  Loran-C  engine¬ 
ers  and  scientists.  The  full  exploitation  of  the  ground  wave  theory  of  prop¬ 
agation  requires  a  data  base  independent  of  the  Loran-C  coordinate  data  base-- 
the  geophysical  data  base  for  the  ground,  References  12,  13,  and  14.  It  is 
fortuitous  that  the  ground  wave  theory,  after  70  years  of  development,  Refer¬ 
ence  14,  can  use  geophysical  data  to  predict  loran  coordinates.  This  leads  us 
to  the  interesting  situation  in  which  we  have  found  a  link  between  the  sciences 
of  geophysics  and  geography  on  the  one  hand  with  that  of  radio  wave  propagation 
about  the  terrestrial  sphere  on  the  other  hand.  It  is  recognized  that  this  link 
really  already  existed,  but  not  in  any  great  detail.  Thus,  since  publication  of 
NBS  Circular  573,  Reference  15,  it  has  been  possible  to  represent  the  geophysical 
condition  of  a  ground  wave  propagation  path  for  Loran-C  as  smooth,  homogeneous 
with  electrical  constants, a,  the  conductivity,  e2,  the  dielectric  constant  and  y2 
the  permeability  of  the  ground.  Thus,  starting  with  Sommerfeld  (1909),  Reference 
37,  it  took  almost  30  years  to  get  a  satisfactory  solution  to  this  problem.  But 
the  advent  of  large  scale  and  now  smaller  but  more  efficient  computers  has  led  to 
a  representation  of  the  boundary  condition  for  the  ground  wave  at  the  surface  of 
the  earth  in  a  more  precise  manner.  Thus,  we  are  now  interested  in  the  nature  of 
the  soil,  the  geological  structure  under  the  soil  and  the  topography  of  the  earth's 
surface.  All  of  this  can  be  used  to  refine  Loran-C  prediction  and  we  believe  such 
is  quite  necessary  if  nanosecond  predictability  is  to  be  obtained. 
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2. 


GEOPHYSICAL  STRUCTURES  FOR  GROUND  HAVE  PROPAGATION 


The  land  masses  of  the  earth  are  nonhomogeneous  and  irregular  in 
varying  degrees,  depending  upon  a  particular  geographic  location.  In  fact, 
one  wonders  in  general  if  the  complicated  earth's  surface  can  indeed  be  re¬ 
presented  by  a  suitable  boundary  condition  for  ground  wave  propagation.  One 
of  the  reasons  for  going  to  a  low  frequency  such  as  100  kHz  used  by  Loran-C, 
is  to  reduce  to  a  negligible  value  the  scatter  fields  from  such  objects  as 
large  mountian  ranges  and  to  reduce  the  attenuation  of  the  signal  over  land 
masses,  expecially  those  that  contain  mountainous  regions.  The  ground  wave 
propagation  problem  then  becomes  primarily  one  of  diffraction  in  the  presence 
of  nonhomogeneous  and  irregular  ground.  While  such  a  problem  is  still  formi¬ 
dable,  it  is  tractable  and  we  believe  can  be  reduced  to  an  operational  prag¬ 
matism  in  the  case  of  Loran-C  chains  that  are  now  operated  over  the  land  masses 
of  the  earth. 

This  paper  reviews  the  specific  geographical  data  that  could  be 
obtained  for  an  area  of  the  Western  United  States  in  which  a  considerable  diversity 
of  such  data  is  possible.  A  geodetic  line  between  the  Searchlight,  Nevada  trans¬ 
mitter  (latitude:  35°19'18.18"  north;  longitude:  114°48'17.43"  west)  and  Ft.  Cronk- 
hite,  California  (latitude:  37°50' 29.464"  north  and  longitude:  122°52 ‘41 . 170"  west) 
was  researched  for  available  geographical  data.  This  work  describes  this  research 
and  pertinent  philosophy  applicable  to  duo-disciplinary  enterprise  of  this  type. 

Consider  first  the  nature  of  the  ground  in  any  particular  locality  as 
it  changes  with  depth  below  the  surface  of  the  topmost  layer  of  soil.  Figure  1 
depicts  a  geoelectric  section  of  ground.  The  Loran-C  signal  propagates  in  the 
half  space  above  such  ground.  Thus,  the  waves  entering  the  ground  are  exponentially 
attenuated  with  depth  or  distance  through  the  material  that  comprises  the  ground,  . 
depending  upon  the  conductivity,  a,  or  the  resistivity,  p.  Although  the  ground  is 
layered  geologically,  such  layering  does  not  necessarily  represent  an  electrical 
boundary  like  the  much  clearer  air-ground  boundary  at  the  surface  of  the  ground. 

But  the  ground  does  vary  in  electrical  resistivity, p.j ,  dielectric  constant,  c- 
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t  =  1,  2,  3...,  and  premeabil ity,  y..  with  depth  below  the  surface.  Thus,  the 
ground  is  anisotropic  in  the  sense  that  the  average  resistivity,  p^,  is  different 
in  the  -w-  direction  below  the  surface  than  the  u-  direction  parallel  to  the  sur¬ 
face.  The  basic  quantity  that  can  be  used  in  the  grdund  wave  theory  of  propagation 
above  the  surface,  Z>0,  is  the  ground  impedance  reflected  to  the  ground-air  sur¬ 
face,  Z  =  0,  by  all  of  the  subsurface  material  of  importance.  With  the  aid  of  this 
impedance  we  can  establish  an  electrical  boundary  condition  for  the  ground  wave  at 
the  physically  obvious  air-ground  interface.  But  how  shall  we  model  the  region 
below  this  interface?  There  is  a  tendency  to  use  formulas  copied  out  of  published 
books  and  papers  to  solve  problems  of  this  type.  When  considering  a  new  endeavor, 
it  is  always  best  to  return  to  physical  fundamentals,  since  the  problem  for  which 
a  given  formula  was  derived  may  not  be  precisely  applicable  to  the  new  problem. 

In  truth,  while  our  solution  for  the  problem  at  hand  may  become  complicated,  it  is 
not  difficult.  We  shall  introduce  some  very  simple  concepts  to  treat  the  geoelectric 
nature  of  the  ground. 

The  impedance,  x,  is  a  complex  quantity  normalized  to  the  impedance  of 
free  space,  Z0-  3/7  ohms.  Thus, 

x  =  Z/Z0  .  (1) 

This  quantity.  Reference  13,  can  be  used  directly  in  connection  with  the  terrain 
irrgeularities  to  satis*''  a  Leonto.vitch  boundary  condition  at  the  earth's  surface. 
The  impedance,  Z,  for  vertically  polarized  waves  calculated  from: 

Z  =  W  (2) 

where, 

'  5V1  [o  zl21  *v 

^u  Zzl  0  \i  (3) 

where  a  u,  v,  w  Cortesian  coordinate  system  for  the  electric  fields,  E;u^v  and 


and  the  magnetic  fields,  £f  ...  has  been  employed.  Then, 

u ,  v 


xCi  =  Z/Z0  =  (1  -  T)(l  +  T)-1  Ci 


(4) 


where  T  is  calculated  from: 
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where  q  =  1,  2,  3,  refers  to  the  three  layer  model  of  the  ground  starting  with 
the  top  electrical  horizon.  Here,  w  =  2nf,  where  f  is  the  frequency,  oq  =  conduct¬ 
ivity  of  qth  electrical  horizon;  =  electrical  horizon  dielectric  constant,  and 

w„  is  the  electrical  horizon  thickness.  C-  is  the  cosine  of  the  angle  of  incidence 
q  i 

of  the  wave  on  the  u-v  surface,  and  S..  is  the  sine.  Reference  13,  gives  the  re¬ 
lationship  of  the  impedance  to  the  propagation  simulator  of  both  the  integral 
equation  type  and  the  classical  type.  are  upgoing  and  downgoing  waves 

(subscript)  respectively  in  the  first  and  second  layer  (subscript)  respectively. 


Of  course,  the  matrix  equation  (5)  can  be  enlarged  to  accomodate 
any  number  of  horizons  of  ground  material  depicted  in  Figure  1.  This  is  a  trivial 
but  tedious  extension.  However,  there  seems  to  exist  in  most  circumstances  a  top 
soil  layer  and  a  geological  basement  rock.  The  soil  may  vary  in  thickness  from 
exposure  of  the  basement  rock  to  several  hundred  meters.  A  typical  value  of  soil 
thickness  seems  to  be  20  meters.  In  regions  of  heavy .rainfall ,  the  soil  is  com¬ 
paratively  homogeneous.  However,  in  arid  regions  and  desert  or  semi-desert  regions 

during  the  dry  season  the  top  part  of  the  soil  seems  to  dry  out  and  a  double  hori¬ 
zon  may  exist.  This  process  reflects  an  impedance  change.  Of  course,  in  special 
areas,  multiple  layers  of  comparatively  precise  distinction  may  be  noted.  However, 
it  is  anticipated  thatthe  three  horizon  model  will  serve  for  initial  estimation  of 
the  impedance  of  the  ground.  Indeed,  it  may  be  very  difficult  to  find  geophysical 
data  in  greater  detail. 

The  effect  of  the  soil  on  the  ground  impedance  is  illustrated  theoret¬ 
ically  in  Tables  1  through  3.  Here  we  have  postulated  a  20  meter  soil  thickness. 

The  topmost  layer  of  soil  conductivity  is  parametric  in  values  between  5  mhos/m 
and  .02  mhos/m  in  Tables  1  and  2  and  between  .01  mhos/m  and  .0001  mhos/m  in  Table  3. 

The  thickness  of  the  top  horizon  has  been  varied  in  each  case  between  0  and  20 

meters,  whilst  the  second  horizon  has  been  varied  between  10  and  0  meters  such  that 
the  total  soil  thickness  has  been  constant  at  20  meters.  Thus,  Figure  1,  if  the 
top  layer  of  soil  were  5  mhos/m  (the  conductivity  of  seawater)  at  a  thickness  of 
only  1.6  meters,  the  surface  impedance,  x  =  0.00105  exp  (jO. 79699)  would  obtain. 
Tables  1A,  IB,  and  the  impedance  would  thereafter  be  controlled  by  the  top  layer 
of  soil. 

If  such  a  highly  conducting  top  horizon  did  not  exist,  i.e.,  the  top 
horizon  was  of  conductivity,  o  =  0.0001  mhos/m  and  the  bedrock,  of  infinite  ex¬ 
tent  and  located  20  meters  below  the  top  of  the  soil  with  a  conductivity,  o  =  0.005 
mhos/m  the  combination  would  reflect  a  very  high  impedance  to  the  surface, 
x  =  0.06911  exp  (jl. 16066).  Now  the  phase  of  the  impedance  can  vary  about  45°  or 
0.79699  radians  (for  seawater)  to  +  tt/2  or  +  90°.  It  is  more  common  in  the  natural 
ground  for  the  impedance  to  approach  +  90°.  Under  such  circumstances,  we  say  the 
ground  impedance  is  inductive.  Less  than  45°  we  call  capacitive.  The  case  just  cited 
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with  a  value  of  1.16066  radians  Is  highly  Inductive  ground,  corresponding  per¬ 
haps  to  a  desert  or  arid  condition  of  the  soil  with  a  more  highly  conducting 
basement  rock  underneath.  Most  arid  regions  of  the  world,  however,  do  on  occasion 
have  rainfall.  One  would  then  expect  considerable  change  in  the  surface  imped¬ 
ance.  If  the  wet  soil  obtained  a  conductivity,  o  =  .02  under  such  conditions, 
one  would  expect  the  ground  surface  impedance  to  decrease  to  0.01639  exp  (jO. 79117) 
such  that  the  ground  impedance  phase  also  reduced  to  a  very  much  less  inductive 
value. 


The  effect  of  such  an  impedance  change  upon  the  propagated  wave  is 
profoundly  a  function  of  the  proportion  of  the  total  propagation  path  involved  in 
the  change.  If  the  entire  propagation  path  is  involved  the  change  in  secondary 
phase  correction  is  given  in  our  previous  report,  CRPL^  77-9,  Reference  13.  Thus, 
at  1000  km  there  would  be  a  change  from  10  ys  to  2.5  ys  or  a  change  of  7.5  ys. 

This  example  of  course,  is  an  extremely  great  change.  In  general,  ground  impedance 
changes,  when  they  occur  are  much  smaller  and  not  of  such  an  extent  as  to  influence 
a  1000  km  propagation  path.  The  application  of  detailed  impedance  changes  to  the 
calculation  of  secondary  phase  correction  and  time  difference  has  been  given  in 
Reference  13,  and  hence  will  be  considered  beyond  the  scope  of  this  report. 

Tables  2A  and  2B  give  the  amplitude  and  phase  of  the  impedance  with  a 
very  low  conductivity  bedrock,  o  =  0.00001  mhos/m  and  more  highly  conducting  and 
variable  soil  on  top.  In  this  case. we  note  the  impedancb  is  more  capacitive,  i.e., 
the  phase  angle  becomes  less  than  tt/4  or  45°.  Finally  in  Tables  3A  and  3B  we 
have  illustrated  a  poorly  conducting  bedrock  with  moderate  to  poorly  conducting 
soil  on  top  of  this  bedrock.  Here  we  note  that  the  conductivity  a  =  0.0001  mhos/m 
on  top  tends  to  make  the  impedance  inductive  as  its  thickness  increases,  and  then 
abruptly  turns  capacitive.  If  such  a  species  of  electrical  horizon  occured,  con¬ 
siderable  secondary  phase  correction  changes  would  be  evident  as  can  be  deduced 
from  Reference  13,  provided  a  great  proportion  of  the  propagation  path  were  involved. 


Table  1A.  Ground  Impedance  Calculations  at  100  kHz  Using  Three-Layer  Model, 

Amplitude, 
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Table  2A.  Ground  Impedance  Calculations  at  100  kHz  Using  Three-Layer  Model, 

Amplitude 
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Table  2B.  Ground  Impedance  Calculations  at  100  kHz  Using  Three-Layer  Model, 

Phase 
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Table  3B.  Ground  Impedance  Calculations  at  100  kHz  Using  Three-Layer  Model, 

Phase 
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3. 


DATA  FOR  GROUND  IMPEDANCE  DETERMINATION 


Figures  2  and  3  depict  soil  and  geological  information  available  for 
western  United  States.  The  particular  propagation  path  between  Searchlight, 

Nevada  and  Ft.  Cronkhite,  California  has  been  indicated  on  these  maps.  Tables  4 
and  5  summarize  the  types  of  soil  along  this  propagation  path  based  upon  infor¬ 
mation  obtained  from  Reference  1,  (U.S.  Department  of  Interior)  and  Reference  10, 
(Soviet  Geography),  respectively.  References  2,  3,  4,  5,  6,  7,  8,  and  9,  give 
further  supporting  information  on  the  nature  of  the  soil.  Table  6  gives  the  geo¬ 
logical  information  for  the  basement  rock  pertinent  to  this  propagation  path. 

This  data  can  be  obtained  from  Reference  1,  (U.S.  Geological  Survey). 

The  correlation  of  the  spatial  perturbations  in  the  Loran-C  phase  (or 
LF  phase)  with  geology  was  noticed  some  time  ago  by  Pressy  et  al  (1953),  Reference 
22.  For  continental  United  States  an  extensive  radio  survey  of  ground  wave  field 
strengths  in  the  low  frequency  band  (broadcast  frequencies  540-1640  kHz)  was  con¬ 
ducted  and  published  as  NBS  Circular  546,  Reference  21.  These  data  were  expressed 
in  terms  of  an  effective  ground  conductivity,  oeff.  This  was  accomplished  by 
matching  ground  wave  attenuation  versus  distance  curves  with  those  observed  in  the 
vicinity  of  various  broadcast  stations.  These  observations  have  been  interpreted 
in  many  ways.  One  such  interpretation  was  accomplished  by  the  Federal  Communications 
Commission  scientists.  Reference  23,  in  the  form  of  a  map  delineating  areas  of 
effective  ground  conductivity,  cr^.  What  exactly  is  meant  by  ag^?  The  ground 
wave  theory  of  propagation,  References  13  and  24,  does  not  use  conductivity, 
mhos/m,  directly.  Thus,  the  quantity  that  enters  into  Maxwell's  equations  is  the 
ground-air  interface  surface  impedance.  If  the  ground  is  of  infinite  extent 
in  the  negative  vertical  direction  and  a  conductivity  a,  mhos/m  is  specified  as 
constant  for  this  entire  half  space,  we  have  the  type  of  ground  used  in  classical 
ground  wave  theory  for  deducing  in  Kirby  et  al,(1958).  Reference  21.  The 
dielectric  constant,  "e2",  and  the  permeability,  "y2"  are  also  specified  as  con¬ 
stant.  Then,  the  impedance,  relative  to  free  space,  equation  (1),  is  given  by: 


X  *  ki/k2  Jl  -  (ki/k2)2j  **  , 

where,  w  =  2irf, 

ki  =  (u/c)Je2  -  j  Og^/coO) 


(6) 


Hi  *  1.000338, 
e2  ^  15, 

£o  a  (c2y0 )  F/m, 
c  =  2.997925  (108)m/S, 

Po  =  4tt  (10~7)  H/m, 

j  =  V^TT 

then, 

°eff  %  62  • 

Some  values  of  ground  impedance  of  interest  are: 
for  seawater, 

x  =  0.0011  exp  (jO.79) 


(7) 


°eff  '  6  -  J° 
for  high  impedance  land, 

x  =  0.040  exp  (j0.81) 

oeff  =  0.0034  -  j  0.0003. 

In  general,  eeff  is  a  complex  number,  but  other  uncertainties  about 
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its  precise  value  usually  do  not  justify  worrying  about  the  imaginary  part,  as 
far  as  Reference  21,  is  concerned.  However,  when  we  change  frequency,  the  non- 
homogeneous  nature  of  the  ground  in  the  negative  vertical  direction  can  affect 
°eff  because  x  changes  profoundly  with  frequency  and  the  nature  of  the  subsurface 
electrical  layers  as  discussed  above  in  connection  with  Tables  1,  2  and  3. 

Thus,  as  a  consequence  a  frequency  dependence  in  was  indeed  noticed.  Then 
the  term  was  coined  to  represent  the  extremely  low  frequency  limit  of  0^(01). 
Thus,  since  the  ground  is  made  up  of  a  multiplicity  of  conductivities,  o,  or 
resistivities,  p,  ohm-m,  it  is  more  logical  to  work  with  the  more  general  quantity, 
the  impedance,  x.  It  is  of  interest  to  note  that  even  the  classical  theory.  Re¬ 
ference  15,  has  an  impedance  calculator  installed  to  convert  a  or  into  im¬ 
pedance,  x.  Thus,  References  15,  in  the  case  of  vertical  electric  polarization, 
the  residues,  s,  of  the  residue  series  are  calculated  from  the  roots  t  =  of 
Riccati's  equation. 


96 

TT  -  2  6*e  T  +1  =  0, 


(8) 


where ,  -J- 

6e  =  -J  (M)  x 

where  a  =  6.36739  (106)m,  the  radius  of  the  earth.  Here  x  is  calculated  as 
shown  in  equation  (6). 

Another  interpretation  of  NBS  Circular  546,  Reference  21,  was  instituted 
by  Keller  and  Frischknecht  (1966),  Reference  20.  In  this  work  we  find  a  very 
interesting  correlation  of  the  conductivity  (resistivity)  with  geologic  age  of 
sedimentary  rocks.  Also,  at  extremely  low  frequencies  (60  Hz)  a  correlation  be¬ 
tween  geologic  age  and  resistivity  was  published  by  Card  (1935),  Reference  18, 
and  Lowy  (1911),  Reference  16.  The  absolute  values  for  the  resistivity  of  the 
bedrock  for  a  particular  geological  age  in  these  cases  is  different.  At  such 
drastically  different  frequencies  one  could  expect  a  difference  in  the  resistivities 
deduced  for  various  reasons,  the  most  important  of  which  may  be  the  depth  of  pen¬ 
etration  of  the  longer  wavelength  into  the  ground.  One  would  expect  the  medium 
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frequencies  of  the  broadcast  band  to  penetrate  to  a  considerably  shorter  dis¬ 
tance.  However,  it  is  interesting  to  note  that  broadcast  frequencies  exhibit 
attenuation  rate  changes  with  the  geological  underlayment  that  correlates  with 
geologic  age,  notwithstanding  the  overburden  or  soil  which  tends  to  decorrelate 
with  the  age  of  bedrock  in  its  formation.  More  important  is  the  corresponding 
strong  correlation  of  Loran-C  frequency  attenuation  rates  and  phase  lags  with 
the  age  of  the  geologic  underlayment.  Reference  22.  In  contrast,  the  higher 
frequencies  seem  to  be  correlated  with  soil  types.  References  17  and  19,  and  in¬ 
deed,  Smith-Rose  (1934)  gives  conductivities  of  various  soil  types  in  England. 

He  also  noted  the  frequency  dependence  of  the  ground  wave  attenuation  rates  be¬ 
tween  1  kHz  and  10,000  kHz.  We  are,  as  a  result  of  such  considerations  led  to 
consider  the  multiple  electrical  horizon  type  soil  depicted  in  Figure  1  as  at 
least  a  triple  electrical  horizon  geoelectric  structure.  If  necessary,  we  can 
use  more  sophistication  later,  but  the  availability  of  pertinent  geophysical  data 
does  limit  us  at  present.  The  improvement  to  be  gained  in  this  model  over  the 
classical  half  space  concept  is  in  our  opinion  enormous,  and  may  be  sufficient 
to  achieve  our  initial  prediction  goals  for  Loran-C. 

The  description  to  follow  is  designed  to  illustrate  procedure  or  the 
beginning  of  an  evolutionary  technique  for  making  initial  estimates  of  the  surface 
impedance  of  the  ground  from  the  rather  vast  information  on  geology  and  soil.  It 
is  important  to  note  that  most  of  the  information  on  geology  and  soil  is  irrelevant 
to  the  problem  at  hand,  but  this  information  in  any  particular  problem  must  be 
sifted  to  obtain  the  very  few  "bits"  of  information  required.  In  our  estimation 
model  we  shall  look  for  three  values  of  conductivity  and  two  values  of  soil  depth. 
It  is  anticipated  that  all  of  this  information  may  not  be  available.  It  is  also 
anticipated  that  in  some  cases  it  may  not  be  necessary.  The  procedure  will  be  one 
of  development  of  an  expertise  that  will  require  judgement  and  intuition;  that  will 
be  improved  by  further  experience  and  will,  of  course,  be  as  quantitative  as  is 
possible.  Our  initial  estimates  for  the  soils  will  consider  only  the  dominant 
orders  and  suborders  depicted  in  Figure  2.  There  are  detailed  delinations  avail¬ 
able  of  various  kinds  of  special  soils  that  do  exist  in  any  of  the  great  orders. 

The  amount  of  detail  we  require  for  our  initial  impedance  estimates  will  thus  be 
quite  restricted,  so  that  the  approach  will  be  practical. 


-17- 


Figure  2.  Illustrating  soil  type  delineation  in  South  Western  United  States  with  Searchlight 
Nevada  to  Ft,  Cronkhite,  California  propagation  path  indicated  in  upper  lefthand 
PART  OF  CHART/  REFERENCE  10. 


The  geological  and  soil  characteristics  of  the  propagation  path  have 
been  delineated  in  Tables  4,  5,  and  6.  Table  7  gives  the  estimated  soil  resistiv¬ 
ities  along  the  propagation  path  from  Searchlight,  Nevada  and  Ft.  Cronkhite,  Cali¬ 
fornia. 


Soil  depths  were  estimated  from  References  29  through  35.  The  hardpan 
below  the  soil  was  included  with  soil  in  most  of  the  cases  considered.  However, 
the  actual  depth  of  the  bedrock  was  quite  uncertain  for  the  literature  surveyed. 

It  is  expected  that  many  of  the  values  given  in  Table  7  may  require  adjustment. 

Between  a  distance  of  0  and  30  km,  we  find  an  area  of  thick  Quarternary 
(Q)  deposits  characteristic  of  western  United  States.  Reference  20,  gives  Gila 
conglomerate  with  a  lithology  description  as  alluvium,  enclosed  debris  filled 
basins  known  as  bol son  deposits,  lake  beds,  and  windblown  sands.  The  correspond¬ 
ing  soil  description  (D2-8)  of  the  general  order  extends  over  this  area  between 
0  and  148  km  is  a  warm,  dry  soil  of  the  ARGID  classification  with  a  mean  annual 
temperature  greater  than  47°F.  The  region  is  quite  arid  for  periods  much  greater 
than  3  months  of  the  year.  Alkali  (sodium)  deposits  exist;  clay  has  accumulated; 
there  is  low  organic  content  in  the  soil.  Indeed,  the  region  comprises  a  subtropi¬ 
cal  desert  soil  that  may  be  intersperced  with  both  mountain  serozems  and  typical 
desert  soils.  These  are  brownish  gray  near  the  surface  and  lighter  colored  below 
and  are  underlain  with  a  carbonate  or  hardpan  layer.  We  set  thickness  at  6  meters 
for  the  dry  top  soil  and  10  meters  for  the  moist  underlayment.  The  dry  top  layer 
estimated  at  2,000  ohm-m  and  the  underlayment  at'100  ohm-m,  including  the  Quarter¬ 
nary  alluvium  type  deposits  of  great  depth.  References  18,  20,  25,  26,  and  27, 
may  be  used  to  assist  with  resistivity  estimations.  At  30  km  we  note  Precambrian 
orthogneiss  and  paragneiss  (Xm),  2,500,000,000  to  1,600,000,000  years  old  as  under¬ 
layment  to  the  ARGID  soils.  Such  underlayment  in  western  United  States  has  varied 
composition  and  origin.  Such  Precambrian  gneisses  may  exhibit  a  resistivity  range 
4,000  ohm-meters  to  250  ohm-meters  on  the  average.  In  fact.  Reference  18,  gives  a 
range  of  1,000-14,000  ohm-m  for  the  intersperced  Cambrian  and  Precambrian  geologic 
ages.  We  shall  use  2,000  ohm-m  for  our  initial  estimate. 

Between  30  and  47  km,  we  find  the  Precambrian  orthogneiss  and  paragneiss 
basement  rocks  that  range  from  2,500,000,000  to  1,600,000,000  years  in  geologic 
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TABLE  4 


Soil  Types  Between  Searchlight,  Nevada  and  Ft.  Cronkhite,  Calitornia 
Based  Upon  U.S.  Department  of  Interior,  Reference  1. 


X 

DISTANCE 

KILOMETERS 

DESCRIPTION 

19.9 

0 

D2-8:  Warm,  dry*  soli  of  ARGID  typo  with  mean  annual  tem¬ 
perature  greater  than  47"F.  Ilaplargids  plus  Paleorthids, 
Torrtpsamments,  Paleargids,  and  Calciortlvids,  gently  sloping 
to  steep.  This  soil  is  an  ARIDISOl  that  includes  soils  that 
have  pedogenlc  horizons  and  are  low  In  organic  matter  and 
are  never  moist  as  long  as  3  consecutive  months.  AGR1DS  are 
a  subgroup  of  ARIDISOLS  that  have  a  horizon  in  which  clay 
has  accumulated  with  or  without  alkali  (sodium);  these  soils 
can  be  used  mostly  for  range  and  some  irrigated  crops. 

3.32 

148.953 

D2-7 :  Warm,  dry  soil  of  ARGID  type  with  mean  annual  tem¬ 
perature  greater  than  47*F.  Ilaplargids  plus  Paleorthids, 
Torriorthents,  and  Rock  land,  gently  or  moderately  sloping. 

A. 15 

173.743 

D2-8:  (soo  above). 

16.58 

204.731 

D2-7:  (see  above). 

8.29 

328.533 

U7-2:  Warm,  dry  soil  of  XLRULT  type  with  mean  annual  tern- 

peraturp  greater  than  47°F.  Haploxerulfs  plus  Xerumbrepts 
and  Xerorthents,  steep.  This  soil  if  an  ULTISOL  that  in¬ 
cludes  soils  that  are  low  in  bases  and  have  subsurface  hori¬ 
zons  of  clay  accumulation;  usually  moist,  but  during  the 
warm  season  of  the  year,  some  arc  dry  par'  of  the  time. 
XERULTS  are  a  subgroup  of  ULT I  SOLS  that  ai  :  relatively  low 
in  organic  matter  in  the  subsurface  horizons.  They  exist 
in  climates  with  rainy  winters  but  dry  summers;  during  the 
warm  season  of  the  year,  these  soils  are  continually  dry 
for  a  long  period;  these  soils  can  be  of  use  as  range  or 
woodland. 

*  Dry  soils  lack  moisture  for  plant  growth  for  long  periods. 


19.9  390.434  Ai 2-2:  Warm,  dry  soil  of  XERAl.F  type  with  mean  annual  tem¬ 

perature  greater  than  4/#r.  Ilaplozeralfs  plus  Palexeralfs 
and  Xerorthants  (shallow),  moderately  sloping  to  steep. 

This  soil  is  an  ALT' I  SOL  that  includes  soils  that  arc  medium 
to  Mali  in  bases  (base  saturation  at  pH  8.2)  and  have  gray 
to  brown  surface  horizon  and  subsurface  horizons  of  clay 
accumulation;  usually  moist  but  during  the  warm  season  of 
the  year  some  are  dry  part  of  the  time.  XERALfS  are  a  sub¬ 
group  of  AlFISOLS  that  exist  in  climates  with  rainy  winters 
but  dry  summers;  during*  the  warm  season  of  the  year  these 
soils  are  continually  dry  for  a  long  period;  these  soils  can 
be  used  for  range,  small  grain  and  irrigated  crops. 

4.15  539.027  03-3:  Warm,  dry  soil  of  ARGID  type  with  mean  annual  tem¬ 

perature  greater  than  47°f.  Natrargids  plus  Nadurargids, 
Haplaquolls,  and  Torr iorthents,  gently  sloping. 

2.99  570.015  A12-2:  (see  above). 

4.64  592.341  E7-2:  Warm  soils  of  0RTIILIIT  type  with  mean  annual  tempera¬ 

ture  greater  than  47°T.  Xerorthents  plus  Durixeralfs,  llap- 
loxeralfs,  Xerof luvenls,  and  Palexeralfs,  gently  sloping. 
This  soil  is  an  ENT 1  SOL  that  includes  soils  that  have  no 
pedogenlc  horizon.  OR THE NTS  are  a  subgroup  of  ENllSOLS 
comprising  loamy  or  clayey  INTI SOLS  that  have  a  regular 
decrease  in  organic-matter  content  with  depth;  these  soils 
may  be  used  for  range  or  irrigated  crops  in  dry  regions  and 
for  general  farming  in  humid  regions.  The  E7-2  belongs  to 
the  subspecies  XrRORTMFMTS  (formerly  Rrger.ols,  brown,  or 
Alluvial  soils)  and  comprise  0RTHCN1S  that  ore  in  climates 
with  rainy  winters  but  dry  summers  and  during  the  warm  sea¬ 
son  of  the  year,  they  are  continually  dry  for  a  long  period. 


TABLE  4  (continued) 


13.6  626.908  A12-2:  (see  above). 

1.66  728.S39  H20:  Saltwater  Sen  Francisco  Harbor. 


0.83  740.934 


INO  746.697 


H16-1:  Kara  soil  of  XEROLL  type  with  mean  annual  tempera¬ 
ture  greater  than  47*F.  Haploxerolls  plus  Arglxerolls, 
Chremoxererts,  and  Xerothants  (shallow),  steep.  This  soil 
Is  a  MOLL I SOI  that  Includes  soils  that  have  nearly  black, 
friable,  organic-rich,  surface  .horizons  high  in  bases;  formed 
mostly  In  subhumld  and  semiarld  warm  to  cold  climates. 

XEROLLS  are  H0LLIS0LS  that  are  In  climates  with  rainy  winters 
but  dry  simmers;  during  the  warm  season  of  the  year,  these 
soils  are  continually  dry  for  a  long  period;  they  may  be 
used  for  wheat,  range  and  irrigated  crops.  Under  the 
XEROLLS,  this  soil  is  classified  under  the  great  group 
(or  phase)  of  HAPLOXEROLLS  (formerly  Chestnut  and  Brown  soils) 
that  are  XEROLLS  that  have  a  subsurface  horizon  high  in  bases 
but  without  targe  accumulations  of  clay,  calciua  carbonate, 
or  gypsum. 

Note:  Gently  sloping— slope  mainly  <10X 

Moderately  sloping— 10X  <  slope  <25X 
Steep— slope  >?5X 


-21- 


TABLE  5. 


Soil  types  Between  Searchlight,  Nevada  and  Ft.  Cronkhite,  California 
Based  on  Soviet  Geography,  Reference  10. 


X 

19.9 


3.32 


4.15 


DISTANCES, 

KILOMETERS  DESCRIPTION 

0  Subtropical  Desert  Soils:  intersperced  with  mountain 

Serozems  and  desert  soils.  The  Serozems  (slerozems)  com¬ 
prise  any  of  a  group  of  zonal  soils  that  are  brownish  gray 
at  the  surface  and  lighter  colored  below,  based  in  a  car¬ 
bonate  layer  or  a  hardpan  layer,  and  characteristic  of 
temperate  to  cool  arid  regions  with  mixed  shrub  vegltation. 

148.953  Solonchak:  Salt  marsh;  any  of  an  Intrazonal  group  of 

stongly  saline  soils;  usually  light  colored  and  without 
characteristic  structural  form  and  typically  developed  In 
poorly  drained  arid  and  semiarld  areas  vegitated  mostly  by 
halophytes. 


173.743  Brown  mountain  soils. 


16.58  204.731  Semidesert  soils. 


8.29. 

328.533 

19.9 

390.434 

4.15 

539.027 

2.99 

570.015 

4.64 

592.341 

13.6 

626.988 

1.66 

728.539 

0.83 

740.934 

END 

746.697 

Brown  mountian  soils.  Mountain  meadow  soils. 

Brown  forest  soils  (largely  neutral)  and  gray  forest  mountain 
soils.  Brown  mountian  soils. 

Brown  soils  of  Xerographytic  Subtropical  forest  and  shrub, 
(same) 

(same) 

Turf -carbonate  mountain  soils. 

(same) 

(same) 
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Illustrating  geologic  basement  rock  type  delineation  in  Western  continental  U 
States  with  Searchlight.,  Nevada  to  Ft,  Cronkhite,  California  propagation  path 
indicated.  Reference  11. 


TABLE  6. 

Geological  Formations  on  Searchlight,  Nevada  To  Ft.  Cronkhite,  California 

Propagation  Path 


X 

DISTANCE 

KILOMETERS 

DESCRIPTION 

4.07 

0 

Q:  Quanternary;  large  areas  of  thick  deposits  of  Western 
United  States. 

2.32 

30.391 

Xm:  Orthogr.eiss  and  paragneiss;  Precambrian,  2,500,000,000 
to  1,600,000,000  years  old;  extensive  in  Western  United 
States  and  of  varied  composition  and  origins. 

4.07 

47.714 

Q:  Quanternary  (see  above) 

0.58 

78.105 

uPz:  Upper  Paleozoic;  Permean  separated  in  part,  P;  UPzc 
clastic  wedge  facies  in  Nevada  and  Idaho,  Diamond  Peak, 

Wood  River,  Milligen  and  associated  formations. 

0.53 

82.436 

Xm:  (see  above) 

0.58 

86.767 

uPz:  (see  above) 

0.58 

91.098 

lPz:  Lower  Paleozoic;  Western  United  States  and  Piedmont 
Province;  €,  Cambrian  separated  in  parts. 

0.58 

95.429 

€:  Cambrian 

0.87 

99.760 

2:  Pre  Cambrian;  800,000,000  years  to  Cambrian;  2-Sedimen- 
tary  rocks,  include  (among  other  unites):  in  Northwestern 
United  States:  Windermere  Group;  in  Great  Basin:  Noonday, 
Johnnie  and  other  formations;  in  Grand  Canyon,  Chuar  Group 
of  Grand  Canyon  subgroup. 

1.16 

106.256 

Q:  (see  above) 

0.58 

114.918 

Y:  Pre  Cambrain:  1,600,000,000  to  800,000,000  to  Cambrian: 
Y-sedimentary  rocks:  Includes  (among  other  units):  in 
Southwestern  Untied  States:  Pahrump  Group,  UntCan  Group  of 
Grand  Canyon. 

2.32 

119.249 

Ti:  Tertiary  intrusive  rocks:  Plugs,  stocks,  and  lacoliths 
of  all  tertiary  ages;  some  granite  plutons  of  early  middle 
tertiary  age. 

0.58 

136.572 

Q:  (see  above) 

1.16 

140.903 

Y:  (see  above) 

3.48 

149.565 

Q:  (see  above) 

1.16 

175.550 

iTv:  Lower  tertiary  volcanic  rocks;  Oligocene  and  Cocene; 
some  Cretaceous  included  in  Arizona  and  elsewhere. 

0.58 

184.212 

Xm:  Pre  Cambrian:  Ortnogneiss  and  Paragneiss. 

O.S8 

188.543 

2:  Pre  Cambrian:  Sedimentary. 

0.29 

192.874 

Xm:  Pre  Cambrian:  Orthogneiss  and  Paragreiss. 

0.70 

195.039 

Xm:  Tertiary  intrusive  rocks. 

0.58 

200.266 

Xm:  (see  above) 

2.32 

204.597 

Q:  (see  above) 

0.73 

221.920 

€:  (see  above) 

0.37 

227.371 

2:  (see  above) 
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TABLE  6  (continued) 


0.46 

233.867 

T1:  (see  above) 

2.32 

237.302 

Y:  (see  above) 

1.74 

254.625 

Q:  (see  above) 

0.35 

267.618 

Qv:  Quarternary  volcanic  rocks. 

1.45 

270.231 

uPz:  Upper  paleozoic. 

0.70 

281.058 

uKjb:  Upper  Cretaceous,  Taylor  Group  (=  Montana  Group); 
in  Montana  plains,  separated  into  UKia  and  b  at  base  of 
Judith  River  formation,  Western  Interior  region. 

0.56 

286.285 

uPz:  Upper  Paleozoic 

0.58 

290.616 

1PZ:  Lower  Paleozoic;  Western  United  States  and  Piedmont 
Province. 

1.74 

294.947 

Qv  :  Quanternary  volcanic  rocks. 

2.32 

307.940 

uKjb:  Upper  Cretaceous. 

1.16 

235.263 

Q:  Quarternary. 

5.23 

333.925 

Kg2:  Upper  Cretaceous  granite  rocks;  80,000,000  to  90,000, 
years  old, --Sierra  Nevada,  California. 

4.07 

372.977 

Kg i :  Lower  Cretaceous  granite  rocks;  105,000,000  to 
120,000,000  years  old;  Sierra  Nevada,  California. 

3.48 

403.368 

lMze:  Lower  Mesozoic,  eugeogyncl inal  deposits,  Middle  and 
Lower  Jurassic  and  Triassic,  Western  United  States. 

1.74 

429.353 

uPz;  Upper  Paleozoic. 

0.58 

442.346 

lMze:  (see  above) 

0.S8 

446.677 

Kg!.’  Lower  Cretaceous. 

0.29 

451.000 

lMze:  Lower  Mesazoic. 

1.16 

453.173 

Kgi:  Lower  Cretaceous. 

1.16 

461.837 

lMze:  Lower  Mesazoic. 

2.90 

470.497 

•Kgr.  Lower  Cretaceous. 

16.26 

492.142 

Q;  Quanternary. 

0.87 

613.555 

Tx:  Tertiary;  Paleocene;  includes  Midway  Group  of  Gulf 
Coastal  Plain. 

7.55 

620.051 

uMze:  Upper  Mesozoic  engeosyncl inal  deposits. 

1.16 

676.426 

Tpc:  Tertiary;  Pliocene  continental  deposits  of  Tertiary. 

2.90 

685.088 

uK:  Mesazoic;  Upper  Cretaceous. 

0.58 

706.742 

urn:  Mesazoic;  Ultramafic  rock  of  Mesozoic  Age  in  Western 
United  States 

2.32 

711.073 

Q:  Quarternary. 

1.74 

728.396 

H20;  Seawater:  San  Francisco  Bay. 

0.70 

741.389 

lMze:  Lower  Mesoioic  eugeosgncl inal  deposits. 

END 

746.697 
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TABLE  7 


ESTIMATED  SOIL  RESISTIVITY,  SOIL  DEPTH  AND  GEOLOGICAL  UNDERLAYMENT  RESISTIVITY 
ALONG  PROPAGATION  PATH  BETWEEN  SEARCHLIGHT,  NEVADA  AND  FT.  CRONKHITE,  CALIFORNIA 


DISTANCE,  GEOLOGICAL  SOIL 

KILOMETERS  TYPE  TYPE  p  p  p  W  W 


age.  These  rocks  are  quite  extensive  in  western  United  States  and  are  of 
varied  composition  and  origins.  From  Reference  18,  an  estimate  of  2,000  ohm-m 
was  selected. 

At  a  distance  of  78  km  along  the  propagation  path  under  consideration 
an  upper  Paleozoic  (uPz)  geologic  basement  is  encountered.  From  Reference  26, 
values  of  300-400  ohm-m  are  noted  and  from  Reference  18,  values  between  8  and 
44  ohm-m.  These  limestone,  bituminous  shale,  clays  and  sandstones  of  the  Permean 
ages  are  quite  variable,  perhaps  due  to  the  water  and  saline  content.  A  value 
of  280  ohm-m  has  been  employed  in  our  initial  estimate.  A  slightly  higher  value 
of  300  ohm-m  was  selected  for  the  lower  Paleozoic  (lPz)  formations  encountered  at 
91  km.  The  more  ancient  Cambrian  formations  according  to  Reference  18,  exhibit  val¬ 
ues  between  6-605  ohm-m,  typified  by  the  Cambrian,  Eocambrian,  Arzberger  series 
or  Saxothuringien  formation  found  in  Europe. 

At  95  km,  the  most  ancient  of  the  Paleozoic  rocks  are  found,  i.e.,  the 
Cambrian  (£)  for  which  a  value  of  400  ohm-m  was  selected  as  an  initial  estimate. 

This  formation  is  followed  by  Precambrian  sedimentary  (Z)  rocks  separated  from  the 
Cambrian  (£)  by  800,000,000  years.  Such  rocks  include  (amongst  other  units)  in 
northwestern  United  States:  Windermere  group;  in  Great  Basin;  the  Noonday,  Johnnie 
and  other  formations;  in  the  Grand  Canyon,  Chuar  group  of  the  Grand  Canyon  subgroup. 
In  general  the  unarticulated  Precambrian  range  between  100  and  5,000  ohm-m.  From 
Reference  18,  we  find  that  Precambrian  and  combinations  of  Precambrian  and  Cambrian 
exhibit  resistivities  between  1,000  and  14,000  ohm-m.  Also,  the  Cambrian  and 
Ordivician  exhibit  resistivities  between  100  and  1,000  ohm-m.  One  may  then  esti¬ 
mate  the  upper  Precambrian  between  450  and  6,500  ohm-m.  Similarily,  at  114  km 
we  find  Precambrian  rocks  1,600,000,000  to  800,000,000  years  old  as  sedimentary 
rocks  (Y)  in  southwestern  United  States  belonging  to  the  Paheump  group,  Unkan 
group  of  the  Grand  Canyon  subgroup.  For  this  lower  group  we  ascribe  a  value  of 
2,000  ohm-m. 

At  119  km  Tertiary  intrusive  rocks  (Ti)  are  encountered,  comprising  plugs, 
stocks,  and  laccoliths  of  all  Tertiary  ages.  Some  granite  plutons  of  early  middle 
Tertiary  age  are  also  found.  According  to  Reference  18,  Cretaceous,  Tertiary, 
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Quaternary  and  combinations  of  these  periods  of  geologic  time  exhibit  resisti¬ 
vities  between  2  and  30  ohm-m.  Quarternary-Tertiary-Pleistocene  deposits  are 
quoted  at  a  range  of  7  to  19  ohm-m.  A  value  of  20  ohm-m  was  selected. 

At  a  distance  of  148  km  from  the  transmitter  a  change  in  soil  order  is 
noted.  The  order  D2-7  is  similar  to  the  D2-8  noted  to  this  point.  Thus,  it  is 
a  warm  dry  soil  of  the  ARGID  type  over  gently  or  moderately  sloping  rock  land. 
Table  4.  The  difference  between  this  soil  and  previous  soils  is  more  evident 
in  Table  5  in  which  we  note  the  soil  is  classified  as  Salonchak,  comprising  salt 

marsh;  belonging  to  a  group  of  strongly  saline  soils,  usually  light  colored  and 

without  characteristic  structural  form  and  typically  developed  in  poorly  drained 
arid  and  semi-arid  areas  vegitated  mostly  by  halophytes.  Thus,  only  salt  resist¬ 
ant  plants  can  survive  in  such  soil.  Notwithstanding  the  high  saline  content  of 
the  soil  we  assign  a  surface  value  of  1,000  ohm-m  because  of  the  existent  desert 
condition.  As  we  go  to  greater  depths  a  resistivity  of  100  ohm-m  was  selected. 

It  should  be  noted,  however,  that  rapid  conductivity  changes  can  occur  if  and 
when  it  does  rain  on  the  desert.  During  the  dry  season,  the  water  table,  or  the 
level  at  which  the  soil  becomes  saturated  with  water,  occurs  at  depths  of  2  to  5 

meters.  This  level  is  quite  unclear  as  a  geoelectric  boundary.  It  is  rather  a 

gradual  increase  in  water  content. 

At  a  distance  of  175  km  from  the  Searchlight,  Nevada  transmitter.  Lower 
Tertiary  volvanic  rocks  (1  Tv)  are  found  with  specific  identities  of  Oligocene  and 
Eocene  and  some  Cretaceous.  These  rocks  range  in  value  between  4  and  31  ohm-m, 
Reference  18.  A  value  of  12  ohm-m  was  selected  as  the  initial  estimate.  Beyond 
this  point  the  geologic  formation  and  soil  repeat  the  various  species  already  dis¬ 
cussed,  to  a  distance  of  267  km.  Here  Quarternary  volcanic  rocks  (Qv)  are  noted. 
These  structures  range  between  10  and  200  ohm-m,  Reference  36.  For  these  rocks 
a  value  of  90  ohm-m  was  selected. 

At  a  distance  of  281  km  from  the  transmitter  we  note  Upper  Cretaceous 
(uK3b)  rocks  of  the  Taylor  group  characteristic  of  the  Montanna  plains.  For 
these  rocks  a  range  in  Reference  18,  between  8  and  17  ohm-m  was  noted.  A  value 
of  10  ohm-m  was  selected. 
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A  soil  order  change  at  328  km  from  the  transmitter  is  noted  (U7-2). 

As  noted  in  Tables  4  and  5,  these  are  characterized  as  brown  mountain  soils  or 
mountain  meadow  soils.  This  area,  328  km  to  390  km  ranges  in  elevation  between 
3,300  and  2,300  meters  above  sea  level.  Table  8,  and  slopes  downhill  in  a  west¬ 
ward  direction.  The  desert  conditions  continue  to  prevail,  but  the  rainfall  and 
moisture  conditions  in  the  soil  are  greater.  A  value  of  500  ohm-m  was  selected 
for  the  top  1.5  to  2  meters  of  soil  and  a  value  of  200  ohm-m  was  selected  for  the 
lower  layers  of  soil.  An  order  (A12-2)  soil  with  similar  electric  properties  was 
noted  at  390  km  from  the  transmitter.  These  are  semidesert  soils  as  described 
in  Tables  4  and  5. 

At  distances  between  333  and  390  km  we  note  Upper  Cretaceous  granite 
rocks  ( Kg 2 ) ,  80,000,000  to  90,000,000  years  old  and  Lower  Cretaceous  granite 
rocks  (Kgi),  105,000,000  to  120,000,000  years  of  a  type  that  are  characteristic 
of  the  Sierra  Nevada  range  of  mountains  in  California.  Granite  rocks  tend  to 
exhibit  high  resistivity.  Reference  (36),  and  range  in  value  between  500  and 
2,000  ohm-m  in  the  Mesozoic  geologic  age.  A  value  of  1,850  ohm-m  has  been  selected 
for  these  granite  rocks. 

At  403  km.  Lower  Mesozoic  (IMHe)  eugeosyncl inal  deposits  (a  syncline 
with  volcanic  materials  mixed  with  sediments)  including  Middle  and  Lower  Triasso- 
Jurassic.  Reference  18,  gives  a  range  of  values  for  the  Triasso-Jurassic  between 
4  and  350  ohm-m.  A  value  of  140  ohm-m  was  selected. 

At  a  distance  of  613  km  from  the  transmitter.  Tertiary  (Tx),  Paleocene 
rocks  were  noted.  Reference  18  gives  values  4  to  31  ohm-m  for  this  geologic  age. 

A  value  of  25  ohm-m  was  selected  for  the  initial  estimate. 

At  a  distance  of  676  km  from  the  transmitter  Lower  Tertiary  (Tpc) 

Pliocene  deposits  were  delineated.  These  rocks  are  sedimentary  and  volcanic  in 
nature.  Reference  36  and  37,  gives  a  resistivity  range  of  10-200  ohm-m.  A  value 
of  150  ohm-m  was  selected.  Also  at  685  and  706  Mesozoic  rocks,  Upper  and  Middle, 
(uK,  urn)  ultramafic  (extremely  buoic,  very  low  in  silica,  an  insulator,  and  /rich 
in  iron,  a  conductor).  A  value  of  15  ohm-m  was  selected  for  these  rocks.  Also 
the  Quarternary  deposits  as  we  approach  the  California  coast  were  lowered  to  a 
resistivity  of  75  ohm-m 
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A.  TREATMENT  OF  TERRAIN  DATA 

Figue  4  depicts  the  terrain  as  the  elevation  above  sea  level,  meters, 
as  a  function  of  distance  in  kilometers  between  Searchlight,  Nevada  and  Ft. 
Cronkhite,  California.  Table  8  (in  APPENDIX)  gives  a  computer  listing  of  the 
data.  These  data  are  also  required  by  the  propagation  simulator  of  the  type 
described  in  CRPL^  Report  77-9,  Reference  13.  The  curve  in  Figure  4  of  elevation 
as  a  function  of  distance  from  the  transmitter  does  contain  discrete  points,  but 
it  is  necessary  to  treat  these  data  as  a  continuous  curve  for  use  in  such  a  prop¬ 
agation  simulator. 

Similarily,  the  soil  and  geologic  data  is  taken  at  discrete  intervals 
as  discussed  in  this  report.  It  is  important  to  note  that  this  also  is  not  meant 
to  imply  sudden  discontinuities  in  the  ground  wave  propagation  medium.  Such  sudden 
discontinuities  require  special  treatment  available  from  other  CRPLi  techniques. 

But  the  natural  ground  is  rarely  discontinuous.  In  fact,  the  various  soils  and 
geological  rocks  seem  to  blend  into  a  continuum. 

There  is  another  fundamental  observation  regarding  the  inter-relation¬ 
ship  of  the  grotind  boundary  condition  at  the  surface  and  the  terrain.  Thus,  the 
soil  thickness,  W1  +  W2,  does  not  form  a  discontinuous  step  function  as  might  be 
implied  from  the  tabulations  of  data  given  in  this  report.  Indeed,  the  soil  is 
usually  studied  and  sampled  in  the  fertile  valleys  where  the  soil  is  quite  deep 
and  suitable  for  agricultural  purposes.  As  we  climb  the  slopes  of  mountains,  the 
soil  slides  back  into  the  valleys  or  at  least  tends  to  do  same.  Thus,  on  steep 
mountain  slopes  we  may  find  very  little  soil.  This  immediately  suggests  the 
smoothing  technique  a>  coupling  technique  between  soil  and  terrain. 

In  the  previous  report.  Reference  13,  we  noted  at  a  distance  point, 

X.j,  the  elevation  point  and  the  terrain  derivatives: 

Zi  =  A  +  BX.  +  CXj  ,  (10) 

^2  =  B  +  ZCX^,  .  (11) 
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Meters 


Figure  4.  Geodesic  profile  from  digitized  terrain  data  base 

ILLUSTRATING  ELEVATION  ABOVE  SEA  LEVEL  AS  A  FUNCTION  OF  DISTANCE 
FROM  THE  TRANSMITTER  OVER  A  PATH  BETWEEN  SEARCHLIGHT,  NEVADA  AND 

Ft.  Cronkhite,  California. 
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=  2C,  1  =  1,  2,  3 . 


where  using  data  of  type  given  in  Table  8, 


•[(xf+rxi-i)'(xS"xt-i),].[  (xi+rlti-i,(xrxi-i>"’j 

b-  |(ir2i.1)-c(x;-x1.1)j-[xrx).1J‘ 


A  =  Z.  -  X.  (B  +  CX^ 
Zi  =  A  +  Xi  (B  +  CX^ 


"5X7  =  B  +  zcxi 


Now,  let  and  W2  be  our  soil  data  base  and. 


R  *  W!/S. 


s  =  W,  +  W2  . 


S  (Xi  ,Zi )  =  s  -  s 


where  ^  is  a  function  of  the  sheer  strength  of  the  soil  arid  the  end  points  of 
the  step  of  the  data  interval.  Then  for  each  small  step  X.  in  the  propagation 
simulator  of  Reference  13, 
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(22) 


(Xj )  s  R  St 

V2  (X1,Zi)  *  S  -  Wit  (23) 

where  W  is  the  precise  geoelectric  section  thickness  used  in  the  impedance 
calculation  at  each  point  X^.  Then,  the  surface  impedance  relative  to  free 
space,  x  =  x^,  equations  (4,5) 

xi  =  f  (Zj.X^  (24) 

A  detailed  discussion  of  these  matters  rests  in  the  domain  of  CRPL^ 
computer  technology  and  hence  is  beyond  the  scope  of  this  paper.  However,  this 
brief  summary  indicated  the  notion  of  both  soil  and  terrain  continuum  and  the 
notion  that  these  continuum  are  coupled  to  each  other  and  to  the  impedance  which 
in  turn  must  also  be  continuous.  If  the  geology,  soil  or  terrain  data  are  not 
continuous,  as  might  well  be  the  case  in  certain  very  anomalous  regions  of  the 
earth,  special  treatment  of  the  problem  is  required  and  CRPL^  will  in  these  cases 
employ  other  appropriate  techniques.  This  also  is  beyond  the  scope  of  this  paper. 
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5. 


CONCLUSIONS 


The  propagation  path  between  Searchlight,  Nevada  and  Ft.  Cronkhite, 
California  is  Indeed  a  "Worst  Case"  type  path  if  diversity  of  terrain  and  geology 
is  a  criterion.  It  certainly  is  not  necessarily  a  "Worst  Case"  path  as  far  as 
effects  on  propagation  are  concerned,  because  we  have  found  that  such  diversity 
often  tends  to  "average  out"  in  some  strange  fortuitous  manner  leaving  one  disap¬ 
pointed  as  to  the  observation  of  the  dramatic  propagation  phenome>  on.  It  is, 
however,  an  excellent  choice  for  demonstration  purposes  and  we  recommend  its  use 
to  demonstrate  the  transient  propagation  and  "envelope  to  cycle  difference" 
phenomenon. 

We  conclude  that  the  dominant  feature  of  the  propagation  path  is  probably 
the  terrain  and  the  geological  structure,  but  we  are  also  of  the  opinion  that  soil 
should  be  considered. 

We  would  like  to  stress  the  continuous  nature  of  the  soil,  terrain  and 
geological  data  if  propagation  simulators  like  the  referenced  integral  equation 
technique  are  employed.  Thus,  the  soil  order  boundaries  and  the  geological  struc¬ 
ture  boundary  should  not  be  treated  as  sharp  discontinuities  unless  special  ap¬ 
propriate  techniques  for  such  matters  are  Introduced. 

We  recommend  that  geological  resistivity  initial  estimates  should  be 
based  on  extremely  low  frequency  observations,  since  radio  measurements  show  a 
frequency  dependence  caused  by  soil  and  rock  stratification. 

We  note  that  the  ground  wave  does  not  propagate  in  the  ground.  Thus, 
sharp  boundaries  in  the  ground  are  permissible  because  the  waves  entering  the 
ground  are  highly  damped.  The  lower  stratum  of  ground  do  reflect  a  surface  im¬ 
pedance  locally  that  is  significant. 

The  number  of  geoelectric  horizons  of  ground  can  be  increased  to  a  greater 
number  than  three  for  more  accuracy  in  the  initial  estimates,  if  geophysical  data 
can  be  obtained  for  this  purpose. 


-34- 


REFERENCES 


U.S.  Department  of  Interior  (1969),  Distribution  of  principal  kinds  of  soils: 
orders,  suborders  and  great  groups;  National  Cooperative  Soil  Survey  Classifica¬ 
tion  of  1967  (U.S.  Geological  Survey,  Washington,  D.C.  20242). 

Soil  Survey  Staff  (1960),  Soil  classification,  a  comprehensive  system,  7th 
approximation  (U.S.  Dept.  Agriculture  Soil  Conserv.  Service,  Washington,  D.C. 
20242);  see  also:  Supplement  to  soil  classification  system  (7th  Approximation). 

Thorp,  J.  and  G.O.  Smith  (1949),  Higher  categories  of  soil  classification— order, 
suborder,  and  great  soil  groups;  Soil  Sci.,  67,  No.  2,  117-126. 

U.S.  Dept,  of  Agriculture  (1938),  Soilds  and  men,  yearbook  of  agriculture  (U.S. 
Dept,  of  Agriculture,  Washington,  D.C.  20402). 

Oakes,  H.,  and  J.  Thorp  (1951),  Dark-clay  soils  of  warm  regions  variously  called 
Rendzina,  black  cotton  soils,  regur,  and  tirs,  (Soil  Sci.  Soc.  America  Proc., 

15,  347-354.  (definition  of  Grumusols). 

Tavernier,  R.,  and  G.O.  Smith  (1957),  The  concept  of  braumerde  (Brown  forest 
soil)  in  Europe  and  the  United  States;  Advances  in  agronomy,  £;  217-281(New  York, 
N.Y.,  Academic  Press  Inc.)  (definition  of  Sols  Brans  Acides). 

Harper,  W.G.  (1957),  Morphology  and  genesis  of  calcisols,  Soil  Sci.  Soc.  America 
Proc.,  21^,  No.  4,  420-424.  (definition  of  calcisols). 

Cline,  M.G.  (1955),  Relationships  among  soils  of  Hawaii,  U.S'.  Dept,  of  Agri¬ 
culture  Soil  Conserv.  Service  Soil  Survey,  Territory  of  Hawaii,  ser  1939,  No. 25, 
67-88  (definition  of  Hzdrol  Humic  and  Humic  Gerruginous  Latosols). 

McClelland,  O.E.,  C.A.  Mogen,  W.M.  Johnson,  F.W.  Schroer,  and  J.S.  Allen  (1959), 
Chernozems  and  associated  soils  of  eastern  North  Dakota— Some  properties  and 
topographic  relationships.  Soil  Sci.  Soc.  America  Proc.,  23,  No.l,  51-56. 

Soviet  Geography  (1965),  Physical -Geographic  atlas  of  the  world,  1964  (Review 
and  translation,  American  Geographical  Society,  Broadway  at  156th  Street, 

New  York,  N.Y. ) 

King,  P.B.,  H.M.  Beikman,  G.J.  Edmonston  (1974),  Geologic  map  of  the  United 
States,  scale  :  :2, 500, 000  (U.S.  Geological  Survey,  Washington,  D.C.). 

Burch,  L.B.,  R.H.  Doherty  and  J.R.  Johler  (1976),  Loran  calibration  by  predic¬ 
tion,  NAVIGATION,  JOurn.  of  Instit.  of  Navigation,  23,  No.  3,  Fall,  195-200. 

Johler,  J.R.  (1977),  Loran-C  propagation  error  corrections  over  nonhomogeneous, 
irregualr  ground  using  the  integral  equation  technique,  CRPL^  Report  77-9, 
U.S.C.G.  Contract  D0T-CG-74629-A,  (Colorado  Research  and  Prediction  Laboratory, 
Inc.,  1898  South  Flatiron  Court,  Boulder,  CO  80301). 


(14)  Samaddar,  S.N.,  The  theory  of  Loran-C  ground  wave  propagation --a  review 
(U.S.  Coast  Guard,  G-D0E-4.TP54;  Washington,  O.C.  20590;  presented  at 
Helsinki  Conference  of  URSI;  to  be  published). 

(15)  Johler,  J.R.,  W.J.  Kellar,  and  L.C.  Walters  (1956),  Phase  of  low  radio  fre¬ 
quency  ground  wave,  NBS  Circular  573  (Supt.  of  Doc.,  U.S.  Gov.  Print.  Off., 
Washington,  D.C.  20402). 

(16)  Lowy,  H.  (1911),  DielektrizitStskonstante  und  Leitfahigkeit  der  Gesteine; 
Annualen  der  Physik,  36,  125-133. 

(17)  Smith-Rose,  R.L.  (1934),  Electrical  measurements  on  soil  with  alternating 
currents;  0.  Inst.  Elec.  Engrs.,  75,  221-237. 

(18)  Card,  R.H.  (1935),  Correlation  of  earth  resistivity  with  geologic  structure 
and  age.  Trans.  Amer.  Inst.  Elec.  Engrs.,  54,  1153-1161, (See  Figure  l,p  387). 

(19)  Barfield,  R.H.  (1934),  Some  measurements  of  the  electrical  constants  of  the 
ground  at  short  wave  lengths  by  the  wave-tilt  method;  J.  Inst,  of  Elec.  Engrs., 
75,  214-220. 

(20)  Keller,  W.J.  and  F.C.  Frischknecht  (1966),  Electrical  Methods  in  Geophysical 
Prospecting  (Pergamon  Press,  New  York,  N.Y.,  see  Table  10,  p.  40). 

(21)  Kirby,  R.S.,  J.G.  Harman,  J.G.  Capps,  and  R.N.  Jones  (1954),  Effective  radio 
ground  conductivity  measurements  in  the  United  States,  NBS  Circular  546  (Suptd. 
of  Doc.,  U.S.  Gov.  Print.  Off.,  Washington,  D.C.  20402). 

(22)  Pressey,  B.G.,  G.E.  Ashwell ,  and  C.S.  Fowler  (1953),  The  measurement  of  the 
phase  velocity  of  ground  wave  propagation  at  low  frequencies  over  a  land  path, 
Proc.  Inst.  Elec.  Engrs.,  100,  pt.  Ill,  73-84. 

(23)  Fine,  H.  (1954),  An  effective  ground  conductivity  map  for  continental  United 
States,  Proc.  IRE,  42,  No.  9,  1405-1408. 

(24)  Johler,  J.R.  (1971),  Loran  radio  navigation  over  irregular,  inhomogeneous  ground 
with  effective  ground  impedance  maps.  Telecommunications  Research  and  Engine¬ 
ering  Report  22  (Suptd.  of  Doc.,  U.S.  Gov.  Print.  Off.,  Washington,  D.C.  20402). 

(25)  Parkhomenko,  E.I.  (1967),  Electrical  Properties  of  Rocks,  (Plenim  Press,  See 
Table  19,  p.  111). 

(26)  Heiland,  C.A.  (1946),  Geophysical  Exploration,  (Prentice-Hall,  see  Tables  65-67, 
pp.  660-664) 

(27)  Wait,  J.R.  (1959),  Overvoltage  Research  and  Geophysical  Applications  (Pergamon 
Press,  see  Table  2,-  p.  106). 

(28)  Soil  Survey  Staff  (1951),  Soil  Survey  Manual,  U.S.  Dept,  of  Agriculture,  Soil 
Survey  Staff,  Handbook  No.  18  (Suptd.  of  Doc.,  U.S.  Gov.  Print.  Off.,  Washington, 
D.C.  20402). 


-36- 


Cosby,  S.W.  and  E.B.  Watson  (1928),  Soil  Survey  The  Gilroy  Area  California, 

U.S.  Dept.  Agr.  Series  1923,  No.  19  (Suptd.  of  Doc.,  U.S.  Gov.  Print.  Off., 
Washington,  D.C.  20402). 

Carpenter,  E.J.  and  Cosby,  S.W.  (1939),  Soil  Survey  of  Contra  Costa  County, 
California,  U.S.  Dept,  of  Agr.  Series  1933,  No.  26  (Suptd.  of  Doc.,  U.S.  Gov 
Print.  Off.,  Washington,  D.C.  20402). 

Carpenter,  E.J.  and  Cosby,  S.W.  (1934),  Soil  Survey  of  the  Svisum  Area, 
California,  U.S.  Dept,  of  Agr.  Series  1930,  No.  18  (Suptd.  of  Doc.,  U.S.  Print. 
Off.,  Washington,  D.C.  20402). 

Watson,  E.B.  and  Storie,  R.E.  (1928),  Soil  Survey  of  the  Bishop  Area,  California, 
U.S.  Dept,  of  Agr.  Series  1924,  No.  3  (Suptd.  of  Doc.,  U.D.  Dept.  Print.  Off, 
Washington,  D.C.  20402). 

Storie,  R.E.  and  Owen,  B.C.  (1940),  Soil  Survey  of  the  Visalia  Area,  California, 
U.S.  Dept,  of  Agr.  Series  1935,  No.  16  (Suptd.  of  Doc.,  U.S.  Gov.  Print.  Off., 
Washington,  D.C.  20402). 

Storie,  R.E.  and  Trussell,  D.F.  (1937),  Soil  Survey  of  the  Barstow  Area,  Cali¬ 
fornia,  U.S.  Dept,  of  Agr.  Series  1933,  No.  8  (Suptd.  of  Doc.,  U.S.  Gov.  Print. 
Off.,  Washington,  D.C.  20402) 

Carpenter,  E.J.  and  Youngs,  F.O.  (1928),  Soil  Survey  of  the  Las  Vegas  Area, 
Nevada,  U.S.  Dept.  Agr.  Series  1923,  No.  8  (Suptd.  of  Doc.,  U.S.  Gov.  Print. 

Off.,  Washington,  D.C.  20402). 

Keller,  G.V.  (1966),  Electrical  properties  of  rocks  and  minerals.  Section  26, 
Handbook  of  Physical  Constnats,  p.  562. 

»• 

Sommerfeld,  A.  (1909),  "Uber  Die  Ausbreitung  der  Wellen  in  der  Drahtlosen 
Telegraphie,"  Ann.  Phys.,  vol .  28.  pp/  665-736:  March. 


TABLE  8 

Elevation  Above  Sea  Level  As  A  Function  Of  Distance  Along  Geodesic  Between  Searchlight 
v.  .  2.  Nevada  and  Ft,  Cronkhite,  California. 


; 


—  -  - 


52# 

531 


8  CONTINUED 


f>7.Zi 


TABLE  8  CONTINUED 


TABLE  8  CONTINUED 


CONTINUED 


1868  SOUTH  FLATIRON  COURT 
BOULDER.  COLORADO  80301 
P.0  BOX  1056 

BOULDER.  COLORADO  80306 
TELEPHONE:  303/444-1530 


COLORADO  RESEARCH  AND  PREDICTION  LABORATORY,  INC. 

SPECIALISTS  IH  ELECTROMAGNETICS  AND  COMPUTER  SIMULATION 


LORAN-C  GROUND  WAVE  SECONDARY  PHASE  CORRECTIONS 
OVER  NONHCMOGENEOUS  AID  IRREGULAR  GROUND 
USING  TRANSIENT  SIGNAL  PROPAGATION  TECHNIQUES 


CRPLj  REPORT  78-12 
FEBRUARY  1,  1979 
REVISED  MAY  1.  1979 
J.  R.  Johler 


Report  on  work  sponsored  by  Department  of  Transportation, 
United  States  Coast  Guard,  Washington,  D.C.  20590 
Under  Contract  No.  D0T-CG-842923A,  dated  August  15,  1978 


FOREWORD 


This  paper,  CRPL^  Report  78-12,  comprises  Task  3B  or  Final 
Report  of  Deliberates,  under  Contract  D0T-CG-842923A,  loran-C  Pulse 
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ABSTRACT 


A  new  rigorous  technique  for  the  computation  and  prediction  of 
loran  time  difference  coordinates  from  Loran  receiver  geographic  coordin¬ 
ates  and  visa  versa  is  introduced  by  application  to  a  specific  loran  over¬ 
land  situation.  The  Loran-C  standard  antenna  current  pulse  has  been  used 
to  describe  the  propagation  of  the  transient  signal  over  ground  that  is 
nonhcmogeneou8  and  irregular.  It  is  demonstrated  that  the  less  sophisti¬ 
cated  and  traditional  continuous  wave  approach  to  the  calculation  of 
secondary  phase  correction  (or  the  propagation  time  correction  relative 
to  the  speed  of  light)  may  be  in  error  over  nonhomogeneous  and  irregular 
land  masses  due  to  the  dispersion  of  the  pulse.  In  this  paper  we  employ 
this  technology  to  the  problem  of  resolving  the  relationship  between  the 
antenna  current  pulse ,  the  near  field  pulse  and  the  far  field  pulse. 

Whilst  a  topographic  data  base  has  been  incorporated  into  the  technique 
for  the  Continental  United  States,  there  exists  a  need  to  introduce 
geologic ,  soil  and  hydrology  data  bases  (in  addition  to  topographic  data 
bases  for  other  parts  of  the  world)  if  this  technology  is  to  be  used  in 
a  prediction  scheme. 

We  have  found  as  a  result  of  the  analysis  presented  in  this 
paper  that  the  rigorous  definition  of  the  envelope  to  cycle  difference 
requires  the  introduction  of  a  concept  of  a  tagged  point  in  time  on  the 
leading  edge  of  the  pulse  envelope  in  addition  to  the  conventional  tagged 
point  in  time  on  the  cycle.  The  behavior  of  the  difference  between  these 
tagged  points  in  time  is  a  function  of  the  propagation  mechanism.  It  is 
suggested  by  the  results  obtained  that  spatial  propagation  error  compensa¬ 
tion  in  the  receiver  may  be  possible  over  very  irregular  ground.  However , 
it  is  recorrmended  that  any  such  compensation  be  implemented  in  the  context 
of  the  theory  presented  to  guarantee  uniqueness. 

The  envelope  to  cycle  difference  (ECD)  of  the  standard  pulse 
manifests  a  near  field  phenomenon  because  its  value  decreases  with  in¬ 
creasing  distance  from  0.42  ye  at  100  and  135  hn  to  a  value  of  0.39  at 
200  km  and  a  value  of  0.19  at  500  km  in  the  presence  of  average  ground 
electrical  constants.  Severe  terrain  features  impose  local  perturbations 
in  the  ECD. 
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LORAN-C  GROUND  WAVE  SECONDARY  PHASE  CORRECTIONS 
OVER  NOMKKENBOUS  AND  IRREGULAR  GROUND 
USING  TRANSIENT  SIGNAL  PROPAGATION  TECHNIQUES 

BY 

J.  Ralph  Johler 


1.  Introduction  To  The  Loran-C  Signal 

The  Loran-C  radio  navigation  chain  is  a  configuration  of  three  or 
more  pulse  transmitters  located  at  precisely  known  geographic  positions  depicted 
in  Figure  1  as  M,  SW,  SX.  These  transmitters  emit  a  radio  pulse  that  has  a 
frequency  composition  cresting  between  90  kHz  and  110  kHz.  The  time  of  trans¬ 
mission  of  each  pulse  is  precisely  synchronized  in  Universal  Time  with  the  aid 
of  cesium  oscillator  clocks  (Hefley,  1972),  Reference  (1). 


SX 


Figure  1.  Loran-C  Chain  depicting  the  Master  (M),  Secondary-W, 

AND  SECONDARY-X  WITH  CORRESPONDING  GEODESIC  DISTANCES  TO  THE 

Observer  (0),*  d^  d^  and  the  baseline  (constant)  geodesics 
dBW  AND  dBX' 


One  of  the  pulse  transmitters  is  designated  the  Master  (M)  and 
the  others  are  called  Secondary-W,  Secondary-X,  Secondary-Y  and  Secondary- 
Z,  SW,  SX,  SY,  SZ.  In  the  past  the  notation  SM,  SI,  S2,  ...  has  been  used. 

The  minimal  chain  for  navigation  requires  only  two  secondaries  as 
shown  in  Figure  1,  as  SW  and  SX.  A  navigator  using  such  a  system  can  locate 
his  position  at  0,  Figure  1,  by  making  a  time  difference  measurement  y, 
with  a  loran  receiver.  Thus,  he  may  use  the  Master— Secondary-W  and  the 
Master— Secondary-X  time  differences,  y. ,  i  =  1,  2  respectively.  If  the 
loran  signals  transmitted  over  the  various  propagation  paths,  d<-w,  dsx,  and 
d^,  traveled  at  the  speed  of  light,  c  =  2.997925(10®)  m/s,  the  equation, 

y^  *  a  constant,  (1) 


would  describe  a  spatial  arrangement  of  hyperbolas  over  the  surface  of  the 
earth,  for  each  i  *  1,  2.  Any  two  hyperbolas  comprising  one  from  each  set, 
i  *  1  and  2,  would  intersect  almost  uniquely  at  the  observer's  position,  0. 
In  general,  equation  (1)  is  not  precisely  true.  The  hyperbolas  are  ever  so 
slightly  distorted  by  propagation  times  that  vary  from  the  speed  of  light 
propagation  time, 

.)  -  t„  (d  „)  (2) 


ri,k 


61  *?■ 


’d1,k  -  d„,kl  *  lc  (d1,l 


where  i  *  1,  2,  and  dQ  =  dM, 


for  the  observer,  k  =  1,  2,  3, 


in  Figure  1.  The  points 
designates  different  geographic  co- 


dj  —  and  d2  ^SX 


ordinates  such  as  <j>k>  Ak>  for  the  latitude  and  longitude.  8i  are  constants 
called  the  emission  delay  of  the  secondaries.  These  constants,  once  selected 
for  convenience  to  avoid  overlapping  of  pulses,  are  held  by  the  chain  manage¬ 
ment.  The  quantities,  tc>  are  the  distortion  corrections  or  the  secondary 
phase  corrections.  Thus,  the  distortion  can  be  corrected  with  a  spatially 
dependent  correction  that  is  usually  expressed  as  tc  in  microseconds  or  <)>c 
in  radians. 


The  precise  nature  of  the  antenna  current  pulse  is  depicted  in 
Figure  2.  This  pulse  has  been  defined  by  the  U.  S.  Coast  Guard  (1974) 
Reference  (2).  The  time  of  arrival  measurement  at  the  receiver,  0,  is  • 
usually  accomplished  on  the  zero  crossing  of  the  third  cycle.  Each  trans¬ 
mitter  emits  pulse  groups  of  eight  pulses  except  the  Master  which  emits  a 
ninth  pulse  for  identification.  Each  group  is  spaced  in  time  so  that  over¬ 
lap  does  not  occur  and  after  some  time  has  elapsed  another  sequential  group 


n 


■■  ■■  ■ 


of  pulses  is  transmitted  at  the  Loran-C  chain  group  repetition  rate.  This 
repetition  rate  varies  such  that  there  exists  40,000  to  99,990  microseconds 
between  Loran-C  chain  pulse  groups  and  gives  the  interesting  illusion  of 
line  spectra  in  frequency  domain  of  the  pulse.  The  position  of  the  various 
spectra  in  frequency  is  dependent  upon  the  repetition  rate.  If  we  consider 
only  a  single  pulse  as  depicted  in  Figure  2,  our  repetition  rates  become 
zero  and  the  line  spectra  merge  into  a  continuum  in  the  frequency  domain. 
Thus,  the  very  powerful  analysis  tool  of  the  Fourier  Transform  Integral  be¬ 
comes  applicable  because  one  can  then  integrate  over  the  complex  (amplitude 
and  phase  are  used)  spectrum.  This  approach  will  be  employed  in  this  paper. 
Finite  repetition  rates  can  be  introduced  if  required  using  discrete  fre¬ 
quencies  under  the  amplitude  and  phase  envelopes  for  a  particular  pulse 
repetition  rate.  Thus,  the  analysis  to  follow  is  completely  general. 

We  perform  the  analysis  with  a  recently  conceived  theoretical 
propagation  analysis  scheme.  This  analysis  scheme  has  been  programmed  on 
a  large  scale  computer  for  use  in  a  proposed  experimental  program  to 
measure  ECD  and  test  the  transient  solution  described  in  this  paper, 
Reference  (21),  (CRPLi  Report  78-10).  This  scheme  uses  a  modular  con¬ 
struction  shown  in  Figure  3.  Thus,  for  example,  the  propagation  transform 
simulator  can  be  the  more  advanced  integral  equation  technique  modified 
for  transient  analysis  and  originally  reported  in  Reference  (9).  Also, 
new  and  more  advanced  techniques  under  study  can  be  inserted  at  this  point 
in  the  analysis  as  a  result  of  the  modular  arrangement  of  the  analysis 
scheme. 

Whilst  the  development  of  such  advanced  computer  technology  is 
important,  it  cannot  be  used  for  precision  general  Loran-C  coordinate 
prediction  without  the  basic  geophysical  data  bases  shown  in  Figure  3. 

The  soil,  hydrology  and  geologic  data  bases  require  considerable  develop¬ 
ment  as  digital  data  bases  at  this  point  in  time.  Thus,  geophysical  data 
base  development  in  a  form  suitable  for  Loran-C  has  not  at  all  kept  pace 
with  the  geophysical  data  requirements  of  ground  wave  theoretical  develop¬ 
ment.  The  soil,  hydrology  and  geologic  data  base  modules  have  been  im¬ 
provised  as  discussed  in  Reference  (20).  This  specially  improvised  data 
base  will  be  used  in  this  paper  to  exercise  the  analysis  scheme. 
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2. 


Signal  Propagation  Time. 


An  ideal  signal  for  radio  navigation  would  appear  to  be  a  contin¬ 
uous  wave  (CW),  Figure  4,  that  can  be  represented  as  a  sinusoid  as  a  function 
of  time,  t.  The  frequency  of  this  wave  is  defined  as:  f  =  oj/2tt  Hertz.  The  wave 
is  uninterrupted  in  time  between  negative  and  positive  infinity  and  the  solution, 
E(t),  along  the  time  coordinate,  t,  for  Maxwell's  equations,  takes  the  form, 

E(t)«  exp(+ jcot)  =  cos(cot)  +  j  sin  (cot),  (3), 

where  j  =  /-l.  Figure  4  can  be  reconstructed  using  the  real  part  of  this  solu¬ 
tion. 

Re  exp  (+  jut)  =  cos  (cot).  (4) 


|  AMPLITUDE 


Figure  H.  Sinusoid  as  a  function  of  time,  t,  depicting  amplitude 

OF  CYCLE  Re  E(t)  AND  ENVELOPE  |E  (t) |  OF  A  CONTINUOUS  WAVE. 


The  imaginary  part  does  not  exist  per  se  physically,  but  may  be  used  to 
synthesize  the  signal  envelope,  |E(t)|, 

|E(t) |®  /cos2  (tut)  +  sin2  (wt)  =  1.  (5) 

We  then  say  the  signal  amplitude  is  unmodulated,  because  its  envelope  is  con¬ 
stant  with  all  time. 

Such  a  signal  can  be  propagated  via  the  ground  wave  propagation 
medium  over  the  distances.  Figure  1,  d<.y,  d^  and  dM  and  from  the  zero  crossing 
of  the  cycles,  phase  information  can  be  obtained  so  that  the  propagation  time, 
equation  (2), 


<n,  d1k)/c  *  tc  (d1k),  (6) 

can  be  deduced.  Since  this  wave  train  has  infinite  length,  it  has  an  infinite 
number  of  zero  crossings,  so  without  further  knowledge  we  could  not  decide  which 
cycle  to  use  in  our  calculation  of  position.  This  basic  problem  in  radio  nav¬ 
igation  is  often  called  "lane  ambiguity"  and  at  100  kHz  the  lane  widths  are  10 
microseconds,  (1  cycle). 


This  problem  is  resolved  by  interrupting  the  wave  train  at  precise 
times  such  that  the  signal  becomes  the  pulse  specified  in  Figure  2.  Mathematically, 
the  pulse  is  specified  as  a  source  function,  F$(t),  such  that  its  Fourier  trans¬ 
form,  fs(u>),  is  given  by: 


V»> 


exp  ( - joot)  Fs(t)  dt. 


(7) 


This  mathematical  approach  to  pulse  synthesis  (1)  permits  introduction  of 
propagation  media  involving  complicated  boundary  conditions;  (2)  permits 


synthesis  of  most  any  pulse  shape  from  continuous  wave  solutions;  and  (3) 
permits  ready  adaptation  to  large  scale  computer  analysis  of  signals.  Such 
approaches  have  been  exploited  In  the  past  by  the  author.  References  (3, 

4,  5,  6,  7,  8,  and  9)  using  a  considerable  number  of  pulse  shapes.  Including 
pulses  from  nuclear  detonations  (EM-pulses).  Thus,  Inversion  of  the  source 
transform,  fs(u>),  with  the  propagated  electric  field  for  a  time  harmonic  source 
transform  £(<u,d),  allows  one 'to  reconstruct  the  transient  field,  E(t',d), 


E(t 


00 

'•d)  =  h  j 


exp  (jut')  £(u,d)/  |  F$(t)  exp(-jut)dt 


>  do)  , 


(8) 


where. 


t'  »  t  -  md/c  (9) 

and  m  =  1.000338  for  air  at  the  surface  of  the  ground.  Here  we  have  re¬ 
cognized  the  spatial  parameter,  d,  representing  distance  from  the  transmitter. 
Equation  (8)  is  valid  for  all  linear  amplitude  analysis  schemes.  Thus,  it 
clearly  applies  to  the  propagation  of  the  Loran-C  ground  wave  pulse.  Indeed, 
we  can  employ  this  analysis  to  calculate  the  propagation  time  corrections 
(secondary  phase  corrections)  for  most  any  condition  of  the  ground  over  which 
the  signal  propagates. 

The  selection  of  the  source  function,  Fs(t),  is  quite  arbitrary. 
Pulses  like  the  Loran-C  pulse  can  be  constructed  from  sinusoides.  This 
technique  uses  complex  source  functions.  Reference  (3).  Suppose  a  carrier 
frequency,  uc>  modulated  by  a  frequency,  up,  has  the  form. 

Re  F$(t)  =  Re  exp  (-vkt) 

=  exp(-ckt)  cos(ukt),  ((k  t<  ~) 

=  0,  (t  <  0)  .  (10) 

Thus,  a  complex  frequency,  vk,  has  been  employed  such  that. 
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where 


\  s  ck  +  H 


k  =  1,  2,  3 . 


(ID 


Vi  -  Ci  +  j(i>i 

v2  *  c2  +  ju>2 

V,  =  c3  +  jo>3 - 


U)l  =  w. 

oj2  =  oj  +  2  a) 

c  p 

w3  =  w  -  2u>  . 

c  p 


Here,  ck,  is  a  damping  constant.  Application  of  equation  (10)  to  equation 
(8)  leads  to: 


E(t',d) 


1 

2? 


*  j("k 


-i 

5(<*>,d)  exp(ju>t')  do). 


(12) 


-00 


The  field  of  a  continuous  wave,  5(w,d)  is  defined: 


E(u>,d)  =  |£(u>,d)|  exp 


or 


C(u>,d)  =  j  |?(w,d)|  exp 


where  $c  is  the  secondary  phase  correction  in  radians.  The  use  of  t'  in 
the  transformation  of  equation  (9)  eliminates  the  factor  exp  [-jkjd]  in 
subsequent  formulations  (see  for  example,  equation  (39)  et.  seq.).  The 
factor  j  arises  from  the  assumed  transmitter  model,  i.e.,  the  Hertzian 
dipole,  equation  (6  AI),  and  is  not  regarded  as  part  of  the  secondary 
phase  correction. 

Upon  separation  of  real  and  imaginary  parts  and  introducing: 
♦'c  =  ♦c  “  7 

we  find: 

00 

E(t',d)  B  ^  |5(w,d)| 
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Using  this  quite  general  form  of  a  damped  sinusoid,  we  can  now  introduce  the 
modulation  frequency,  uy  as,  for  example,  a  sine  squared  type  of  modulation: 


Re  F$(t) 


Re  Fs(t) 


=  exp  (-Cit)  sin2  (wpt)  sin  (u»ct) 
=  h  exp  (-Ci t)  sin  (u)ct) 

-  H  exp  (-Cit)  sin  (u>c  +  2uip)t 


-  h  exp  (-Ci t)  sin  (wc  -  2wp)t. 


(15) 


Therefore, 


Fs(t)  =  (j/2)  exp(»-V|t)  -  (j/4)  exp  (-v2t) 


-  (j/4)  exp  (-v3t), 


(16) 
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where 


Vi 

*  Cl 

+ 

J“c  =  Cl 

+  j<*>! 

V2 

*  Cl 

+ 

j(“c  +  2 

V  3 

C2 

+  jd)2 

Vs 

*4 

O 

N 

+ 

J(U)C  -  2 

u>p).  * 

C3 

+  j<»>3 

Then, 

E(r,d)  =  (j/2)  E(t',d,vi )  -  (o/4)  E(t',d,v2)  -  (j/4)  E(t',d,v3)  (18) 

Figure  5  depicts  a  ground  wave  pulse  propagated  to  a  distance  of 
1000  km  from  the  source  over  average  ground  with  a  surface  impedance,  x  =  0.33 
exp  (j0.777)  using  the  equations  (14)  and  (15)  and  the  constants, 

C!  =  2500 

fp  =  2500  Hz  =  u)p/2ir 

f  =  100  kHz  =  (u>/2t0/103  . 

Both  the  envelope  ±  |E(t',d)|  and  the  cycles  under  the  envelope,  i?eE(t',d) 
are  depicted  as  a  function  of  local  time,  t'.  The  amplitude  is  in  volts 
per  meter  for  a  one  ampere-meter  source  current  dipole.  This  pulse  is 
similar  to  the  Loran-C  pulse  depicted  in  Figure  2.  Whilst  we  have  shown 
the  entire  pulse,  it  is  necessary  to  magnify  detail  of  the  pulse  to  study 
the  effects  of  propagation  via  the  ground  wave.  Ionospheric  waves  distort 
this  pulse  in  practice  and  the  reader  is  referred  to  Reference  (8)  for  the 
ionospheric  wave  pulse.  This  latter  pulse  is  of  no  consequence  to  Loran-C 
under  normal  conditions,  since  the  pulse  is  measured  for  navigation  timing 
purposes  at  the  last  zero  crossing  of  the  third  cycle.  This  occurs  before 
the  arrival  of  the  ionospheric  wave. 

It  can  be  concluded  from  equation  (16)  that  the  sine  squared 
pulse  modulation  is  merely  the  sum  of  three  sinusoids  with  the  complex 
frequencies  vlt  v2,  and  v3  (the  carrier,  the  upper  and  the  lower  sidebands). 

It  is  Interesting  to  note  that  this  is  perhaps  the  purest  form  of  amplitude 


modulation.  Thus,  any  reduction  In  the  upper  or  lower  s 'deband  relative  to  each 
other  or  to  the  carrier  can  result  in  a  concommi ttant  phase  modulation.  Since 
the  sidebands  are  not  propagated  with  equal  vigor,  some  phase  modulation  can  be 
Imposed  by  the  time  harmonic  transform  of  the  field,  5(w,d).  Computer  simulation 
studies  of  such  modulation  has  been  accomplished  depicting  both  ground  waves  and 
Ionospheric  waves.  Reference  (5,  8,  9). 


TIME,  MICROSECONDS 

Figure  5.  Illustrating  computer  synthesis  of  a  ground  wave  pulse 

BY  SUPERPOSITION  OF  DAMPED  SINUSOIDS  AT  A  DISTANCE  OF  1000  KM  FROM 
THE  TRANSMITTER  OVER  AVERAGE,  SMOOTH  GROUND. 


3.  Relationship  Between  Near  Field  and  Far  Field  Envelope 

and  Cycle  Propagation  Time. 


The  U.S.  Coast  Guard  (1974),  Reference  (2),  has  given  a 
standard  form  of  pulse.  Figure  2,  such  that  the  Loran-C  modulation  can 
be  written: 

Re  F$(t)  ■  bjt2  exp  (-Cit)s1n(wct),  (t>0),  (19) 

*  0  (t<0) 
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where  Ci  *  (2/65) ( 10s )  and  bj  *  1.753(10*)  for  1  ampere-meter  current  pulse 
crest.  For  500  kilowatts  radiated  power  at  the  pulse  crest,  bi  *  9. 576( 101 3 ) . 
This  leads  to  a  complex  source  function, 

Fs(t)  *  jbit2  exp  (-Vjt)  (20) 

*  bit2  exp  (-cit)  |sin(wct)  +  j  cos  («ct)j  (t^O), 

fs(u)  *  2jbi(vi  +  ju>)"*  «  2jbj  ci  +  j((oc  +<*>)]"*.  ,21) 


This  modulation  can  be  completely  represented  by  a  slight  modification  of 
equation  (14)  as  follows: 


E(t',  d)  *  4*>i 


_  _3_ 

[ci2  +  (uc+  u)2J  2  cost  ut'  -  4'  +  3  tan’lJj(u>c  +  u)/ci~j  } 


+  Cc?  +  (“c '  2  cos^a)t" +  +  3  tan_1C-K  ■ 


+j  {  l  cf  +  (w  +  U>)2]  2  sin  tut'  -  4>  '+  3  tan_1[-(a)  +  io)/Ci]} 


+  [cf  +  U  -  u>)2]  2  sin{-wt'  +  ♦'  +  3  tan"  [-(«c  -  wj/ci]}  >  du.  (22) 


In  this  specific  case  a  single  damped  sinusoid  is  all  that  is  required  to  rep¬ 
resent  the  modulation. 

In  general,  given  the  arbitrary  modulation,  F$(t),  that  we  have 
described  above,  equation  (8)  can  be  written, 

00 

E(t',d)  *  57  /  EU.d)  exP  (Jut')  f$  (u>)  dw,  (23) 

-  00 

and  a  periodic  modulation  envelope  has  a  transfonv.. 


fs(«) 


F$  (t)  exp  (-jut)  dt. 
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where  T,  Is  the  period  of  the  envelope.  The  Loran-C  receiver  can  tag  points 
In  time  on  both  the  pulse  cycle  and  the  pulse  envelope.  This  can  be  simulated 
mathematically.  The  equations  for  amplitude  Independent  tagged  points  In  time 


|E(t',d)|  -  Ce^  |E(t'.d)|  =  0,  and 

3pr|E(t',d)|  =  0, 


(25  a) 


(25  b) 


where  Is  the  constant  that  determines  a  particular  tagged  point  in  time 
on  the  envelope.  The  inflection  point  can  be  determined  by  the  2nd  derivative. 
The  cycle  zero  crossings  are  of  course  the  roots  of: 

Re  E(t',d)  =  0.  (26) 

Such  quantities  can  be  calculated  from  the  above  described  theory.  Thus,  the 
roots  of  equation  (25) 

t-  •  Tc  ,  (27) 

and  the  roots  of  equation  (26), 

t'  *  tc  (28) 

give  the  envelope  and  cycle  propagation  time  corrections  respectively.  The 
envelope  to  cycle  discrepancy  or  the  envelope  to  cycle  difference,  r,  (ECD) 
is  then  the  difference  between  envelope  and  cycle  propagation  corrections, 

T  and  t  respectively,* 

r  =  Tc  -  tc  •  (29) 

We  select  in  our  loran  receivers  a  zero  cycle  crossing  that  is  closest  to 
T  ,  the  envelope  propagation  time  correction.  Thus,  CQ  in  equation  (25)  is 
set  so  that  Tc  is  close  to  the  last  zero  crossing  of  the  third  cycle.  This 
quantity,  when  derived  at  various  geographic  positions  for  constant  Ce, 
is  a  measure  of  the  pulse  dispersion  by  the  propagation  medium.  The  propaga¬ 
tion  time  of  the  envelope,  ue,  and  the  cycle,  u,  are  given  by: 


ue  =  Tc  +  hid/c, 
u  =  tc  +  nxd/c  . 


(29  a) 
(29  b) 


The  envelope  derivatives  can  be  calculated  from: 


E(t',d)  =  ^  E(t',d)  j  cos  (arg  ^  E(t',d) 

-  arg  E(t',d)}  . 


Equations  (25  a,  26,  27,  28  and  29)  comprise  a  rigorous  definition 
of  envelope  to  cycle  difference  (ECD). 
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Also,  It  follows  that: 

AZ 

cos  (arg  ^pr  E(t',d) 

-  arg  E(t',d)}  . 

+  { | E( t",d)  |  }ml{  E(t',d)|2}  sin2{arg^  E(t',d) 

-  arg  E(t',d)}  .  (31) 

The  derivatives  indicated  by  (30)  and  (31)  are  calculated  from  the  analytic 
function  for  the  transient  field,  equation  (14),  for  example. 


d2 

dF’ 


E(t',d) 


E(t',d) 


_d 

dt 


.  < 


E(t',d)  =  ^  J  |5(u)id)| 

o 

'[ck  +  ^“k  +  U)^2]sin  "  ^c^1^  +  +  w)/ckj) 


+ 

ck  +  ^wk  " 

%  r 

sin  {-cot'  +  <j>'(<Did)  +  tan  1  - 

■  H  '  a))/ck] 

+  j 

ck  +  W  + 

-f  *5  | 

cos  {ut'  -  (^(uijd)  +  tan1 

-  (a>k  +  w)/ckj 

} 

1 

“*5 

-(ajk  -  cj)/ck 

ll 

ck  +  (“k  ”  w)2 

cos  {-(lit'  +  d>'(coxd )  +  tan1 

1 

ojdo),  (32) 


and. 


d2 

W* 


00 

E(t',d)  =  J?  /  l5(“‘d)l 


-jck2+  (uk  +  oj) 2J  cos  (o»r  -  <p'  (w,d)  +  tan  j^-(a>k  +  w)/ckj 


-  ck2+  (wk  +uj): 
ck  +  ^wk  "w^2 


l  r  "^5  “  ^  r 

-  ck  +  (<»>k  ■  <*02  cos  t"  <*>t'+<J>'  (u>,d)  +  tan  |-(“k  -w)/ck  ) 


sin  {cut'  -  <p'(u,d)  +  tan  j^-(wk  +  uO/cJ  } 
sin  (-ut'  +  <i>'(w,d)  ►  tan'1  |j-(u>k  -  w)/ckJ  } 


(33) 
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u)2db)  . 


The  effect  of  the  finite  bandwidth  of  the  Loran-C  receiver  is 
ordinarily  negligible.  However,  filters  of  various  types  are  often  employed. 
These  filters  distort  the  pulse  and  affect  the  measured  propagation  times, 
u,  ug.  The  effect  of  a  filter  can  be  readily  assessed  with  the  aid  of  this 
analysis  by  introducing  into  equations  (22)  or  (23),  the  transform  of  the  filter, 
fr(w): 


E(r,d) 


C(u>,d)  exp  (jut')  fr(<i)) 


f$  (<jd)  dw. 


(34) 


A  filter  typically  has  the  form: 

fr(u>)  =  (jojai  [-a?+  juai  +  at"]"1}!1  (35) 

where,  w  =  2"  (10s),  on  =  2.5  (10s),a2  *  4.1  (107). 

The  quantity  n  =  1,2,3...  can  be  increased  to  reduce  the  bandwidth  of  the 
filter  and  assess  the  effect  of  bandwidth  upon  the  pulse  dispersion* 


In  equation  (28)  we  have  rigorously  deduced  the  secondary  phase 
correction,  tc.  This  quantity  plus  the  propagation  time  at  the  speed  of  light 
in  air,  md/c,  is  u,  the  true  propagation  time  of  the  cycle.  It  is  this 
quantity  that  the  Loran-C  receiver  measures  with  precision  (usually  nanosecond 
precision).  It  is  important  to  note  that  t  is  a  function  of  distance  from  the 
transmitter,  d,  and  time,  t'  , 

tc  =  tc  (d,r).  (36) 

It  has  been  standard  engineering  practice  to  use  the  last  zero  cros¬ 
sing  of  the  third  cycle  under  the  pulse  envelope  a.s  depicted  in  Figure  4.  Thus, 
this  point  in  time  is  approximately  30  microseconds  from  the  beginning  of  the 
pulse.  It  can  never  be  precisely  this  amount  because  the  pulse  is  subject  to 
dispersion  by  the  propagation  medium.  The  dependence  of  tc  upon  t'  is  curious. 
Thus,  in  the  past  it  has  been  common  practice  to  derive  the  secondary  phase  cor¬ 
rection  from  the  continuous  wave  carrier  frequency,  f  =  100  kHz,  Figure  3.  This 
is  not  an  exact  secondary  phase  correction.  However,  various  numerical  studies 

*  An  extensive  study  of  loran  pulse  dispersion  due  to  a  receiver  filter  that  in 
turn  was  convolved  with  the  ground  wave  propagation  medium  as  shown  in  equation 
(34)  is  given  in  Reference  1 5 
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References  (3,  5)  have  validated  the  approximation  of  equation  (31): 

tc  -  tc(d).  (37) 

at  least  for  the  cases  of  smooth  homogeneous  ground  and  nanosecond  pre¬ 
cision.  In  the  case  of  irregular  and  nonhomogeneous  ground  this  approxi¬ 
mation  may  require  further  scrutiny. 

Indeed,  whilst  in  previous  works.  References  (3,  5,  9)  we  found 
very  weak  dependence  of  the  envelope  to  cycle  difference,  i.e.,  the 
response  of  the  dispersion  measure  r  to  changes  in  the  secondary  phase 
correction, 

r  =  r  (tc);  (38) 

a  strong  dependence  was  noted  locally  in  the  presence  of  perturbations 
due  to  mountains  (ridges)  and  ground  impedance  nonhomogenities.  Under 
such  conditions  one  wonders  about  the  validity  of  the  approximation  in 
equation  (37).  These  matters  will  be  expounded  using  the  full  transient 
solution  developed  for  this  project.  Thus,  our  secondary  phase  correc¬ 
tions  are  the  rigorous  equation  (36)  and  we  shall  use  the  third  cycle  of 
the  pulse  like  Loran-C  receivers  ordinarily  do. 

The  field  representation  close  to  the  transmitter  for  a  damped 
sinusoid,  Reference  (3),  can  be  written  as  a  transform  for  the  Sommerfeld- 
Norton  surface  wave.  References  (11,  12): 

£(s,d).  f$(s)  =  C  Zifa  exp(s2a)erfc(s/x) 

+  1  .  1 

to(S+v)  t0:is(s+v) 

where  s  =  jw.  Let  I0t  =±  1  ampere-meter^  the  dipole  current  moment 
corresponding  to  a  1  ampere  current  in  a  1  meter  antenna,  then, 

C  =  2I0«02/(4tt  ed3)  =  2(10-7)/d, 

e  =  coGi  =  CoHi2* 

and  neglecting  the  dielectric  constant  C2,  of  the  ground,  the  numerical 
distance  factor,  ot,  is: 

*  The  eign  reversal  on  the  antenna  current  on  successive  pulses  is  known 
as  phase  coding  in  Loran-C  navigation  parlance. 


a  =  n  3d/(2apoC3), 

i 


and  the  propagation  time  in  air  at  the  surface  of  the  ground,  t0,  is 

to  =  Hid/c. 

Now,  a  complex  inversion  integral. 


E(t',d)  =  (s,d). 


is  implied  where. 


f (s)  =  ZF(t 


00 

"•I 


exp(-st')  F(t')dt'. 


Thus,  after  considerable  ado, 

E(t',d)  =  vC  exp(-vt')  ^  -1  -v  /fra  exp(v2a) 

•  erfc(-v/a)-erfc  — -  v/a  1  +  — - - —  \ 

2/a  /_  tov  to2v2  J 

+  vC  ( — - — )  +  C  —  +  v  exp  |  — — |  . 

|to2v2/  _2a  __  \  4a  / 

At  very  short  distances,  d,  the  field  can  be  written: 


5(s,d).  f  (s)  =  C  {—5—  +  - - -  +  — - 

s  )  S  +  V  t0(s+v)  t02s(s+v) 


Whereupon  the  transform  inversion  gives: 

E(t",d)  =  v  C  exp  (-vt')  (  -1  +  -i- 

I  t0v 


t02v2 


t02v2 


+  Cui(t'); 


where 


u0(t)  =  1  t>0, 
■  0  t<0, 


u0(t)  =  f  Ui(t)  dx. 


and 


Thus,  Ul(t),  is  the  Dirac  function,  which  can  be  written: 


Uj(t')  =  Lim 
ot+o 


h +  v]  exp(is — }  * 


Independently,  equation  (42)  reduces  to  equation  (44)  for  ot  =  o  and  equa 
tion  (44)  becomes  independent  of  the  electrical  constants  of  the  ground. 
This  implies  a-**>  or  d-K).  Thus,  equation  (39)  can  be  calculated  from  the 
rigorous  expression.  Reference  (10): 

Pi  3-ju>niX2d/2c 


1  -  k,2/k22) 

. 


Pi  =  Hkx3/2k22))  1  -  kj2/k22)  d  =  u2a  , 


(45) 


As  d  grows  small  the  first  two  terms  of  equation  (39)  approach 


s  _  .  _  v 
S+V  ”  s+v  * 


and  it  reduces  ultimately  to  equation  (43).  In  equation  (42), 
erfc(x)  l 

Alternatively,  using  the  U.S.  Coast  Guard  Loran-C  pulse  specified  in 
equation  (19),  equation  (39)  can  be  rewritten: 


5(s 


,d).  fs(s)  =  bjC  ^  exp  (s2a)  erfc(s^) 


L 

,  1  .  1 
to (S+v) 3  to^S+v)* 


or  equation  (43)  becomes 


£(s,d).  fs(s)  bxC  <f  ^s+v)3  +  t0(s+v) 3  +  ‘T77sTs+vPj  * 


The  inverse  transform  of  equation  (47)  then  is: 


E(t',d)  =  bx  vC  exp  (-vO  |  (-}* 

+  (v  ■  uhr)  +  (W}  >  +  bivC  Wj  ‘ 


(46) 


(47) 


(48) 
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This  equation  (48)  comprises  the  near  field  representation  of 
the  U.S.  Coast  Guard  Loran-C  current  pulse  specified  by  equation  (19). 
There  exists  an  electrostatic  or  non-slnusoidal  component  to  this  pulse. 
Here,  v  a  cx  +  ja)c  where, 

Cl  *  (2/65 ) ( 10s ) , 

u>c  *  2ir(l05). 


The  modulus,  ±  |E(t',d)| , would  give  the  pulse  envelope  as  a 
function  of  time.  The  "real  part  of"  the  analytic  function,  jteE(t'.d) 
would  give  the  cycles  under  the  pulse  envelope  as  a  function  of  time.  Here 
it  should  be  emphasized  that  time  is  local  time,  t',  equation  (9).  Thus, 
causality  specifies  that:  t'  >  0.  This  fundamental  attribute  of  electro¬ 
magnetic  theory  was  Investigated  by  Sommerfeld  and  Brlllouln  (1914), 
References  (13  and  14). 

The  finite  conductivity  of  the  ground  in  the  vicinity  of  the 
transmitter  can  be  introduced  into  equation  (48)  through  the  factor,®, 
and  we  find  the  inverse  transform  of  equation  (46)  is. 


E(r  ,d)  =  blV  C  exp(-vt')  ‘  2W-]  +  [£  +  4av^“t^Vj 


-2a  (1+2  Va+  4a2  v2)  - 


1 


+  (-2avi/a)t'  +  2a(l  +  4a2v2  +  ZVS.)  exp(v2a)  exp 


[(27S  *  v,/5)2| 


-v/x)t'2  +  (1  +  8a2v2  +  2/x)t'  -  (8a2v  +  +  6av/i  +  16a*v2/jt)] 

v  j 


Ml 

•  /ffa  exp  (v2a) 


t' 


erfc(-v»/o T)  -erfc(^j-  -  v/x) 


+  b^C 


t0V 


(49) 


A  more  exact  expression  for  equation  (49)  Is  derived  in  Appendix  I. 

Finally,  the  rigorous  solution  of  the  U.S.  Coast  Guard  pulse  for 
all  distances  from  the  transmitter  in  the  presence  of  nonhomogeneous  and 
Irregular  ground  is  given  by  equation  (22).  A  solution  is  obtained  using 
special  numerical  techniques  on  the  computer.  The  rigorous  solution  can. 
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of  course,  be  checked  numerically  at  short  distances  against  computations  on 
equation  (49)  to  give  confidence  In  the  numerical  mastery  of  the  problem. 

The  rigorous  solution  to  the  problem  currently  uses  In  the  transform  5(u>,d) 
calculation,  a  Volterra  Integral  equation  of  a  type  discussed  In  our  pre¬ 
vious  CRPLj  Report  77-9,  Reference  (17).  This  rigorous  approach  to  the 
calculation  of  the  transform  has  been  greatly  Improved  recently  and  Is  more 
efficient  numerically  than  the  classical  Sommerfeld-Norton  formula,  equa¬ 
tion  (46). 


Equation  (25)  establishes  a  tagged  polnt-ln-time,  t',  rigorously 
on  the  pulse  envelope,  at  either  an  arbitrary  point  for  which  a  value  of 
Ce  Is  selected  or  at  the  special  Inflection  point  prescribed  by  the  second 
derivative  with  respect  to  time.  Whilst  such  tagged  points  In  time  have 
been  Implemented  In  the  past,  modern  loran  receivers  do  not  attempt  an 
exact  envelope  tag.  In  fact,  it  Is  probably  not  necessary  under  most  prac¬ 
tical  circumstances.  The  AN/BRN-5  Loran  receiver  serves  as  an  example  of 

A 

modern  receiver  envelope  detection.  Reference  (16).  In  this  receiver,  r  T 
such  that. 


T  =  f 


r  1 

Re  E(t;,d)/Re  E(tf,d) 


or 


-16 


-i  -i 


Re  E(t2,d)/Re  E(tr,d)  +  0.7 


+  20. 


(50) 


A 

where,  T  Is  expressed  In  microseconds  and,  using  the  last  zero  crossing  of 
the  4th  cycle, 

t2  •  tc  +  37.5 


ti  -  ti  -  5  =  t£  +  32.5 

A 

T  Is  measured  using  the  last  zero  crossing  of  the  4th  cycle  as 
the  value  of  t"2  +  2.5.  As  a  consequence  of  the  secondary  phase  correction 
and  pulse  cycle  dispersion  t"2  +  2.5  will  not  be  precisely  40  microseconds. 
Then  the  envelope  to  cycle  difference,  r.  Is  thus, 

T  «  T  (51) 
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If  Re  E(t'i,d)/Re  E(t'ltd)  Is  ^  1.33,  then  T  -  12.12  ys  and 

a 

changes  In  T  will  be  proportional  to  r.  This  Is  not  exact  because,  the 
root,  t'  •  T0,  of. 

Re  E(t",d)  *  0  , 

or  the  crest  of  the  half  cycles  t"  a  ti  and  t*  are  not  exact  at  t_  +  37.5 
and  tc  +  32.5  respectively  because  the  pulse  may  be  stretched  In  time  due 
to  dispersion.  Thus,  strictly  speaking,  t  Is  not  constant  along  the 
leading  edge  of  the  pulse  due  to  the  propagation  mechanism.  However,  this 
procedure  is  accurate  enough  to  resolve  cycle  or  lane  ambiguity. 

The  various  techniques  used  by  receiver  manufacturers  to  asertaln 
envelope  time,  Tc  or  T,  would  enter  Into  any  experiment  to  validate  the 
pulse  dispersion  theory  presented  In  this  paper.  Therefore,  the  theory 
should  be  validated  using  cycle  zero  crossings  and  crests  along  the  leading 
edge  of  the  pulse.  From  such  measurements  the  rigorous,  T  ,  or  the  other 

^  V 

envelope  tags  like  T  can  be  synthesized. 

It  is  possible  to  observe  the  Loran-C  pulse  at  an  area  monitor 
and  adjust  the  antenna  current  pulse  for  the  ideal  pulse  at  this  point 
In  the  service  area.  Suppose  the  distance  to  the  monitor  is  di.  Then, 
the  antenna  current  pulse  can  be  deduced: 


Fs(t)  =  H  \ 
—00 


f(w,d,) 


[^(w.dJJ  exp  (jut)  du 


(52) 


where  the  current  pulse  is. 


and  the  envelope  Is, 


Re  F$(t), 


|Fs(t)|  • 


Thus,  T 

f  (u,di)  -  j  FCt",di )  exp  (-jut')  dt' 
0 


* 

F(t',< 

4 


(53) 


-  22  - 


where  F  (t'.d)  Is  the  observed  pulse  at  a  distance,  dlt  from  the  trans¬ 
mitter.  If  F  (t',d)  looks  like  an  Ideal  pulse,  the  transform  must  have  a 
spectrum  like  equation  (21), 

f  (u>,  dx)  =  b2 
and,  according  to  equation  (20), 


Ci  +  j(«>c  +  w) 


Ett'.dJ  «  jb2t'2exp  (rvat')  . 


(54) 


where  b2  Is  the  amplitude  constant. 

Furthermore,  the  transient  field  at  some  other  distance,  say 
d  =  d2,  can  also  be  predicted: 


E(t',d2) 


1 

2tt 


J 


f(w,d1)  S(o>,d2) 


•t-i 


C(w,di) 


exp  (jut')  dw 


(55) 


Thus,  from  an  observed  pulse  at  some  point,  di  distance  from  the  transmitter, 
one  can  predict  the  antenna  current  waveform.  Re  F$(t),  or  the  field  at  some 
other  distance,  E(t',d),  d  s  dlf  d2,  ...  .  Here  again  the  analytic  function, 
E(t',d)  comprises  the  pulse  cycles. 

Re  E(t',d)  , 

and  the  pulse  envelope. 


±|E(t',d)|  . 

It  is  also  quite  possible  that  the  area  monitor  pulse  does  not 
conform  to  the  ideal  pulse.  In  this  case  also  the  actual  pulse  can  be  digi¬ 
tized  and  the  transform  calculated  numerically: 

T 

f (<u,di )  *  j  exp  (-juit')  Re  F(t',d,)  dt',  (56) 

0 

where  T  is  the  period  of  the' pulse  envelope.  The  complex  source  function, 
F$(t)  can  be  deduced  from  distant  observations  using  equation  (52). 

Therefore,  equations  (52)  through  (56)  are  important  for  validation 
of  the  theory  of  propagation  of  the  Loran-^C  pulse. 
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Let  us  now  consider  the  relationship  between  the  very  special 
antenna  dipole  current  moment  pulse  given  in  Figure  2  and  specified  by 
equation  (19)  at  some  distance  (far  field)  from  the  transmitter.  Figure  6 
illustrates  such  a  field  in  volts  per  meter  as  a  function  of  local  time, 
t't  defined  by  equation  (9).  This  is  the  far  field  that  would  be  generated 
at  a  distance  of  100  km  over  smooth,  homogeneous  spherical  ground  with  a 
conductivity,  a  =  .005  mhos/m.  The  Intensity  is  given  In  absolute  units 
for  an  antenna  current  moment  that  reaches  a  crest.  Figure  2,  of  54630 
amperes  in  a  one  meter  antenna  at  65  microseconds  or  218.5  amperes  in  a 
250  meter  antenna  uniformly  distributed  over  the  length  of  the  vertical 
structure.  There  exists  small  corrections  at  very  short  distances  for 
antennas  with  nonuniform  current  distribution,  especially  those  that 
approach  A/4  wavelengths  and  it  is  a  comparatively  simple  matter  to 
Introduce  most  any  such  current  distribution  as  demonstrated  in  Reference 
(15)  for  a  monopole.  The  top  loading  of  Loran-C  antennas  also  make  the 
problem  of  source  representation  somewhat  more  complex.  Such  a  problem 
Is  beyond  the  scope  of  this  paper  and  this  research.  Using  the  formula: 

Pr  =  1.6(10-13)  u>2  (I0t)2/«0.  *  1 .67548(10"4)  watts,  radiated  power 

where  B0  -  377  ohms,  the  impedance  of  free  space,  and  I 0t  *  1. 

Then,  for  Pr  =  500,000  watts, 

V'Pf;(I7~  "V~1?67S^TID~V)  =  5*463(104)  Ampere-meters  =  94.75  db. 

Table  3  gives  numerical  values  for  the  cycle  zero  of  the  pulse  at 
100  km.  It  is  necessary  to  add  2.5  microseconds  which  has  been  subtracted 
out  by  the  radiation  mechanism  from  the  antenna  and  to  subtract  the  cycle 
number  times  10  microseconds  to  get  the  secondary  phase  correction  for  each 
cycle.  Thus,  in  the  case  of  the  first  cycle, 

9.03  +  2.5  -  10  =  11.53  -  10  =  1.53. 

For  cycle  number  1.5, 

13.82  +  2.5  -  15  =  16.32  -  15  =  1.32. 
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Tables  1  and  2  give  the  corresponding  near  fields  for  0.1  km  and 
1  km  respectively.  A  computer  oscillogram  of  the  near  field  at  0.1  km  is 
also  given  in  Figure  8,  depicting  both  the  cycles.  Re  E(t",d)  and  the 
envelope  |E(t',d)|.  The  third  cycle  error  relative  to  a  continuous  wave 
signal  varies  between  -0.08  vs  and  0.11  vs.*  It  is  also  evident  from 
Tables  1  through  4  that  the  distortion  is  an  early  time  phenomenon.  The 
distortion  is  not  obvious  from  an  oscillogram  of  the  pulse  such  as  depicted 
in  the  computer  graphs  of  Figures  6  and  7,  d  =  100  and  500  km  respectively. 

It  is  necessary  to  magnify  the  detail  of  the  pulse  to  observe  the  distortion. 
That  the  distortion  exists,  is  evident  from  the  assumed  antenna  current 
pulse,  the  zeroes  of  which  are  precisely  at  10,  15,  20,  25,  ...  vs. 

The  advantages  of  using  a  pulse  in  preference  to  a  CW  (continuous 
wave)  signal,  as  previously  pointed  out  in  the  introduction,  far  outweighs 
the  distortion  discussed  here.  This  distortion,  however,  should  make  us 
very  cautious  in  our  treatment  of  non homogeneous  and  irregular  ground  as 
proposed  in  Reference  (19)  as  an  ultimate  prediction  scheme.  It  is  also 
obvious  that  such  errors  as  discussed  above  tend  to  cancel  in  a  time 
difference  coordinate  observation.  But  since  a  spatial  dependence  in  the 
error  exists,  such  errors  do  not  cancel  exactly. 

The  effect  of  nonhomogeneous  and  irregular  ground  was  illustrated 
for  purposes  of  this  paper  by  introducing  a  specially  prepared  geophysical 
and  geological  data  base  discussed  in  CRPL-j  Report  78-9,  Reference  (20). 

Thus,  a  preliminary  operation  of  the  computer  configuration  shown  in  Figure 
3  was  accomplished  for  a  propagation  path  between  Searchlight,  Nevada 
and  Tecopa,  California,  a  distance  of  approximately  137  km.  The  newly 
acquired  topographic  data  base  for  Continental  United  States  was  employed 
for  this  specific  propagation  path. 


*  It  is  of  special  interest  to  observe  that  the  third  cycle  error 
relative  to  a  continuous  wave  (CW)  at  0.1  km.  Table  1,  is  -0.08 
ys  (microseconds).  At  1  km.  Table  2,  the  error  is  0.11  vs.  At 
100  km.  Table  3,  the  error  is  0.09  ys.  Finally,  Table  4,  at  500 
km,  the  error  is  0.03  ys.  Thus,  the  phase  modulation  of  the  pulse 
over  smooth,  homogeneous  ground  appears  to  be  a  near  field  phenomenon. 


-  25  - 


I .  -  -  ..  ..v .. ..  — • — — BHBI 


Time,  Microseconds 

FIGURE  6.  Theoretically  propagated  Loran-C  pulse  at  a  distance  of  100  km  from  the  antenna  source  curren- 
pulse  given  in  Figure  2,  assuming  the  ground  to  be  smooth  and  homogeneous  with  a  conductivity. 
°  =  .005  mhos/m  or  ground  kith  an  impedance,  Ixl  *  .055. 


Re  E(t',d) 


TABLE  1 


ONE  TENTH  KILOMETER  NEAR  FIELD  CYCLE  ZERO  AND 
SECONDARY  PHASE  CORRECTION,  MICROSECONDS, 


FOR  VARIOUS  ZEROES  OF 

THE  CYCLE,  Re  E(t',d) 

=  0 

Cycle 

Number 

Local  Time,  t',  ys 
of  Cycle 

Zero 

Cycle 

Secondary  Phase 
Correction,  tc,  ys 

Error 

Relative  to 
CW  Secondary 
Phase  Correction 

1 

11.86 

4.36 

-0.33 

1.5 

16.94 

4.44 

-0.25 

tc  (CW)  ■ 

2 

22.04 

4.54 

-0.15 

4.69 

2.5 

27.07 

4.57 

-0.12 

3 

32.11 

4.61 

-0.08 

3.5 

37.12 

4.62 

-0.07 

4 

42.15 

4.65 

-0.04 

4.5 

47.16 

4.66 

-0.03 

5 

52.17 

4.67 

-0.02 

5.5 

57.18 

4.68 

-0.01 

6 

62.18 

4.68 

-0.01 

Time,  Microseconds 

FIGURE  7 .  Theoretically  propagated  Loran-C  pulse  at  a  distance  of  500  km  from  the  antenna  source  curi 

PULSE  GIVEN  IN  FIGURE  2,  ASSUMING  THE  GROUND  TO  BE  SMOOTH,  HOMOGENEOUS  WITH  A  CONDUCTIVITY. 
°  *  .005  MHOS/M  OR  A  GROUND  IMPEDANCE,  Ixl  =  .033. 


% 


Re  E(t',d) 

/  ±lE(t:.d)| 


TABLE  2 


ONE  KILOMETER  NEAR  FIELD  CYCLE  ZERO  AND 
SECONDARY  PHASE  CORRECTION,  MICROSECONDS, 

FOR  VARIOUS  ZEROES  OF  THE  CYCLE,  Re  E(t',d)  =  0 


Cycle 

Number 

Local  Time 
t  ,  ys , 
of  Cycle 

Zero 

Cycle 

Secondary 

Phase 

Correction,  t  ,  ys 

V 

Error 

Relative  to 

CW  Secondary 
Phase  Correction 
tc-tc(CW) ,  ys 

1 

8.94 

1.44 

0.96 

1.5 

13.75 

1.25 

0.27 

2 

18.68 

1.18" 

0.20 

1.5 

23.62 

1.12 

0.14  t  (CW)  =  i 

3 

28.59 

1.09 

0.11 

3.5 

33.56 

1.06 

0.08 

4 

38.54 

1.04 

0.06 

4.5 

43.52 

1.02 

0.04 

5 

48.51 

1.01 

0.03 

5.5 

53.50 

1.00 

0.02 

6 

58.49 

0.99 

0.01 

.98 
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TABLE  3 


ONE  HUNDRED  KILOMETER  FAR  FIELD  CYCLE  ZEROES 
AND  SECONDARY  PHASE  CORRECTION,  MICROSECONDS, 
FOR  VARIOUS  ZEROES  OF  THE  CYCLE,  Re  E(t',d)  =  0 


Cycle 

Number 

Local  Time,  t',  ys 
of  Cycle 

Zero 

Cycle 

Secondary  Phase 
Correction,  tc,  ys 

Error 

Relative  to 

CW  Secondary 
Phase  Correction 
tc-tc(CW),  ys 

1 

9.03 

1.53 

0.48 

1.5 

13.82 

1.32 

0.27  tc(CW)  = 

2 

18.73 

1.23 

jg  1.05  ys 

2.5 

23.68 

1.18 

.13 

3 

28.64 

1.14 

.09 

3.5 

33.62 

1.12 

.07 

4 

38.60 

1.10 

.05 

4.5 

43.59 

1.09 

.04 

5 

48.58 

1.08 

.03 

5.5 

53.57 

1.07 

.02 

6 

58.56 

1.06 

.01 
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TABLE  4 


FIVE  HUNDRED  KILOMETER  FAR  FIELD  CYCLE  ZEROES 
AND  SECONDARY  PHASE  CORRECTION,  MICROSECONDS, 
FOR  VARIOUS  ZEROES  OF  THE  CYCLE,  Re  E(t',d)  =  0 


Cycle 

Number 

Local  Time,  t',  ys, 
of  Cycle 

Zero 

Cycle 

Secondary 

Phase  Correction, 

V 

Error 

Relative  to 

CW  Secondary 
Phase  Correction 
VycW),  ys 

1 

10.66 

3.16 

0.19 

1.5 

15.57 

3.07 

0.10  tc(cw)  = 

2 

20.54 

3.04 

°*07  2.97  ys 

2.5 

25.51 

3.01 

0.04 

3 

30.50 

3.00 

0.03 

3.5 

35.49 

2.99 

0.02 

4 

40.48 

2.98 

0.01 

4.5 

45.48 

2.98 

0.01 

5 

50.47 

2.97 

0.00 

5.5 

55.48 

2.97 

0.01 

6 

60.46 

2.96 

-0.01 
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DISTANCE,  KILOMETERS 


Figure  9  depicts  the  Irregular  nature  of  the  ground  over  the 
Searchlight,  Nevada  to  Tecopa,  California  propagation  path  as  a  profile 
of  ground  elevation  in  meters  relative  to  sea  level  as  a  function  of 
distance  in  kilometers  from  the  transmitter.  The  dotted  line  across  the 
figure  references  the  elevation  of  the  transmitter.  The  corresponding 
secondary  phase  correction  for  a  continuous  wave  signal,  t  (CW),  using 
ground  impedances  calculated  from  the  geophysical  data  base,  is  also 
depicted  on  this  graph  for  reference.  Since  this  report  is  oriented 
towards  the  transient  solution,  its  understanding  and  characteristics, 
this  should  not  be  interpreted  as  a  prediction.  Prediction  requires  much 
greater  detailed  scrutiny  of  the  geophysical  data  base.  While  the  second¬ 
ary  phase  correction  variations  follow  in  large  measure  the  main  features 
of  the  terrain,  there  also  exists  minor  variation  due  to  the  ground  im¬ 
pedance  changes  along  the  propagation  path. 

Using  the  antenna  current  pulse  specified  by  equation  (19), 
together  with  equation  (22),  the  cycle  zeroes  at  a  distance  of  135.64 
km  from  the  transmitter  were  examined  as  displayed  in  Table  5.  This 
table  gives  the  local  time,  t'>  in  microseconds,  the  cycle  secondary  phase 
correction,  tc,  ps  and  the  error  relative  to  assuming  a  continuous  wave 
(CW)  secondary  phase  correction,  tc(CW)  =  1.53  ps,  as  a  function  of  cycle 
number  or  various  zeroes  of  the  cycle.  We  note  again  the  rather  interesting 
phenomenon  of  the  secondary  phase  correction  approaching  the  continuous 
wave  secondary  phase  correction  at  the  cycle  number  6.  At  the  cycle 
number  3,  where  most  loran  receivers  accomplish  a  precision  tag,  the 
error  is  0.11  ps.  Note  that  this  error  is  due  to  the  propagation  mechanism 
and  may,  therefore,  be  expected  to  display  a  spatial  distribution.  Thus, 
at  this  point  in  the  argument  one  can  conceive  the  existence  of  a  110 
nanosecond  (0.11  ps)  error  using  continuous  wave  secondary  phase  correc¬ 
tions!  This  can  never  be  satisfactory  if  precision  predictions  are  con¬ 
templated. 

Let  us  carry  our  numerical  survey  of  the  problem  further  and 
consider  Table  6.  Here  we  find  smooth,  homogeneous  ground  with  an 
average  ground  Impedance,  |x|  =  .033,  also  exhibits  an  error  of  90 
nanoseconds.  Furthermore,  in  the  previously  referenced  Tables  1,  2,  3, 

4,  corresponding  to  distances  of  0.1,  1,  100,  and  500  km  over  average. 
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smooth,  homogeneous  ground,  we  get  errors  of  -80,  110,90  and  30  nanoseconds 
respectively.  It  is  apparent  that  whilst  the  irregular  and  nonhomogeneous 
ground  may  have  a  disturbing  effect  upon  the  cycle,  the  error  can  also 
exist  over  smooth,  homogeneous  ground  in  a  quite  deceptive  manner. 

It  is  also  evident  that  this  phenomenon  is  a  function  of  the 
local  time,  t and  the  problem  could  be  resolved  by  using  the  cycle  number 
6  zero  instead  of  the  cycle  number  3  zero.  But  this  cycle  may  be  contamina¬ 
ted  by  sky  wave  which  would  cause  for  greater  error  than  is  now  under  dis¬ 
cussion.  Hence,  it  is  necessary  to  sample  early  at  the  third  cycle  and 
accommodate  this  error  in  the  prediction  scheme. 

Tables  7  and  8  give  this  error  at  other  distances  along  the  pro¬ 
pagation  path  depicted  in  Figure  9.  It  is  surprisingly  constant  between 
the  distances  selected  and  this  may  be  fortutious,  since  the  error  then 
tends  to  cancel  in  a  time  difference  measurement,  at  least  to  the  resolu¬ 
tion  of  many  loran  receivers.  This,  however,  should  not  be  relied  upon 
in  a  prediction  scheme  since  fortuitous  cancellation  of  error  is  not  a 
general  principle. 

Tables  9  and  10  further  enhance  the  numerical  survey  of  the 
transient  solution  for  the  zeroes  of  the  third  cycle.  A  tagged  point  in 
time  on  the  envelope  using  equation  (25  a)  and  an  arbitrary  value, 

Cg  =  2.904(10" 5 ) .  This  sets  the  envelope  tag  at  approximately  30  ps  at 
d  =  0.1  km.  It  is  evident  from  Table  9  that  the  envelope  to  cycle  difference, 
r,  in  microseconds  can  be  readily  calculated  for  the  various  yalges  of  the 
ground  parameters  (smooth,  homogeneous  ground,  SHG,  or  irregular,  nonhomogeneous 
ground  (ING)  and  at  various  distances.  It  is  noted  that  the  enyelope  tag, 

T  ,  mimics  the  cycle  tag,  t  ,  and  as  a  result  the  envelope  to  cycle  dif- 

L  v 

frence,  r,  tends  to  reflect  only  weakly  the  effects  of  spatial  changes  in 
the  propagation  medium.  Of  course,  a  more  severe  topography  would  give 
stronger  envelope  to  cycle  differences.  Whilst  this  is  suggestive  of 
using  a  form  of  receiver  automatic  compensation  for  spatially  distributed 
propagation  errors,  we  offer  the  opinion  that  any  attempt  at  such  com¬ 
pensation  should  be  accomplished  in  the  context  of  the  theory  of  transient 
propagation,  If  any  degree  of  uniqueness  is  to  be  guaranteed  by  such 
compensation. 


TABLE  5 


FAR  FIELD  CYCLE  ZEROES  AND  SECONDARY  PHASE 
CORRECTION,  MICROSECONDS,  IN  THE  PRESENCE  OF 
NONHOMOGENEOUS  AND  IRREGULAR  GROUND  BETWEEN 
SEARCHLIGHT,  NEVADA  AND  TECOPA,  CALIFORNIA  (135.64  km), 
FOR  VARIOUS  ZEROES  OF  THE  CYCLE,  Re  E(t',d)  *  0 


Cycle 

Number 

Local  Time,  t',  us, 
of  Cycle 

Zero 

Cycle 

Secondary 

Phase  Correction, 
tc,  us 

Error 

Relative  to 

CW  Secondary 
Phase  Correction 
tc-tc(cw),  us 

1 

9.56 

2.06 

0.53 

1.5 

14.34 

1.84 

0.31 

2 

19.24 

1.74 

0.21  tc(CW)« 

2.5 

24.17 

1.67 

0.14  1.53 

3 

29.14 

1.64 

0.11 

3.5 

34.11 

1.61 

0.08 

4 

39.09 

1.59 

0.06 

4.5 

44.07 

1.57 

0.04 

5 

49.06 

1.56 

0.03 

5.5 

54.05 

1.55 

0.02 

6 

59.04 

1.54 

0.01 

V 

E 


TABLE  6 


FAR  FIELD  CYCLE  ZEROES  AND  SECONDARY  PHASE 
CORRECTION,  MICROSECONDS,  FOR  SMOOTH,  HOMOGENEOUS 
GROUND  WITH  IMPEDANCE,  |x|  -  0.033  AND  A  DISTANCE, 
d  ■  135.64  km 


Cycle 

Number 

Local  Time,  t',  ys, 
of  Cycle 

Zero 

Cycle 

Secondary 

Phase  Correction, 
tc.  ys 

Error 

Relative  to 

CW  Secondary 
Phase  Correction 
tc-tc(CW),  ys 

1 

9.19 

1.69 

0.45 

2.5 

14.01 

1.51 

0.27  tc(CW)  = 

2 

18.92 

1.42 

0.18  1.24 

1.5 

23.86 

1.36 

0.12 

3 

28.83 

1.33 

0.09 

3.5 

33.81 

1.31 

0.07 

4 

38.79 

1.29 

0.05 

4.5 

43.78 

1.28 

0.04 

5 

48.77 

1.27 

0.03 

5.5 

53.76 

1.26 

0.02 

6 

58.75 

1.25 

0.01 

! 
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TABLE  7 


FAR  FIELD  CYCLE  ZEROES  AND  SECONDARY  PHASE 
CORRECTION,  MICROSECONDS,  IN  THE  PRESENCE  OF 
NONHOMOGENEOUS  AND  IRREGULAR  GROUND  AT  A  DISTANCE 
OF  99.736  km  FROM  THE  SEARCHLIGHT,  NEVADA  TRANSMITTER 


Cycle 

Number 

Local  Time,  t',  ys 
of  Cycle 

Zero 

Cycle 

Secondary  Phase 
Correction,  tc>  ys 

Error 

Relative  to 

CW  Secondary 
Phase  Correction 
tc-tc(CW),  ys 

1 

9.16 

1.66 

0.55 

1.5 

13.92 

1.42 

0.31  tc(CW)  = 

2 

18.82 

1.32 

0.21  1.11 

2.5 

23.75 

1.25 

0.14 

3 

28.71 

1.21 

0.10 

3.5 

33.68 

1.18 

0.07 

4 

38.66 

1.16 

0.05 

4.5 

43.65 

1.15 

0.04 

5 

48.63 

1.13 

0.02 

5.5 

53.63 

1.13 

0.02 

6 

58.62 

1.12 

0.01 
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TABLE  8 


FAR  FIELD  CYCLE  ZEROES  AND  SECONDARY  PHASE 
CORRECTION,  MICROSECONDS,  FOR  SMOOTH,  HOMOGENEOUS 
GROUND  WITH  IMPEDANCE  |x|  =  0.033  AND  DISTANCE, 
d  =  99.736  km 


Cycle 

Number 

Local  Time,  t',  ys 
of  Cycle 

Zero 

Cycle 

Secondary 

Phase  Correction 

V  vs 

Error 

Relative  to 

CW  Secondary 
Phase  Correction 
tc-tc(CW),  ys 

1 

9.02 

1.52 

0.47  tc(CW)  = 

1.5 

13.82 

1.32 

0.27  1.05 

2 

18.73 

1.23 

0.18 

2.5 

23.67 

1.17 

0.12 

3 

28.64 

1.14 

0.09 

3.5 

33.62 

1.12 

0.07 

4 

38.60 

1.10 

0.05 

4.5 

43.59 

1.09 

0.04 

5 

48.58 

1.08 

0.03 

5.5 

53.57 

1.07 

0.02 

6 

58.56 

1.05 

0.01 
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TABLE  9 


CYCLE  AND  ENVELOPE  LOCAL  TIME  t' ,  FOR  VARIOUS  DISTANCES 
OVER  SMOOTH,  HOMOGENEOUS  GROUND  (SHG)  AND  IRREGULAR, 
NONHOMOGENEOUS  GROUND  (ING)  SHOWING  DERIVATION  OF 
ENVELOPE  TO  CYCLE  DIFFERENCE  (ECD)  AND  SECONDARY  PHASE 
CORRECTION  OF  ENVELOPE  AND  CYCLE,  T£  AND  tc  RESPECTIVELY 

r  = 

3rd  Cycle  Envelope  e 

Distance  Type  Zero  Local  Tagged  Point  2.904(10“5) 


d,  km 

Ground 

Time,  t',  ys 

In  Time, 
t  T  ys 

ECD 
r,  ys 

V 

V  ^ 

0 

Antenna 

Current 

30.00 

30.68 

0.68 

0.00 

0.68 

0.1 

SHG 

34.61 

29.98 

-4.63 

4.61 

-0.02 

1 

SHG 

31.09 

29.33 

-1.76 

1.09 

0.67 

2.9921 

SHG 

30.37 

30.48 

0.11 

0.37 

0.48 

2.9921 

ING 

30.53 

30.51 

-0.02 

0.53 

0.51 

9.9736 

SHG 

30.48 

30.70 

0.21 

0.48 

0.70 

9.9736 

ING 

30.38 

30.78 

0.40 

0.38 

0.78 

40.8917 

SHG 

30.73 

31.19 

0.46 

0.73 

1.19 

40.8917 

ING 

30.90 

31.28 

0.40 

0.90 

1.30 

99.7359 

SHG 

31.14 

31.56 

0.42 

1.14 

1.56 

99.7359 

ING 

31.21 

31.69 

0.48 

1.21 

1.69 

118.6857 

SHG 

31.24 

31.66 

0.42 

1.24 

1.66 

118.6857 

ING 

31.67 

32.04 

0.39 

1.65 

2.04 

135.6408 

SHG 

31.33 

31.75 

0.43 

1.33 

1.76 

135.6408 

ING 

31.64 

32.13 

0.49 

1.64 

2.13 

200. 

SHG 

31.65 

32.04 

0.39 

1.65 

2.04 

500. 

SHG 

33.00 

33.19 

0.19 

3.00 

3.19 

SHG  A  Smooth,  homogeneous  ground,  |x|  =  .033 

ING  A  Irregular,  non homogeneous  ground  as  specified  in  CRPLi  Report 
78-9,  Reference  (20) 

ECD  a  Envelope  to  Cycle  Difference,  r  =  T£-t  . 
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TABLE  10 


COMPARISON  OF  VALUES  OF  ECD,  r,  SECONDARY  PHASE  CORRECTIONS 
FOR  CYCLE,  tc,  ENVELOPE,  T.,  AND  CONTINUOUS  WAVE,  tc(CW),  IN 
THE  PRESENCE  OF  SMOOTH,  HOMOGENEOUS  GROUND  (SHG) ,  AND 
IRREGULAR  NONHOMOGENEOUS  GROUND  (ING) 

SMOOTH,  HOMOGENEOUS  GROUND  IRREGULAR,  NONHOMOGENEOUS  100  kHz  CONTINUOUS 


DISTANCE, 
d,  kM 

(SHG) 

ECD, 

r,  us 

|x|  = 
Cycle 
Tc.  os 

.033 

Envelope 

V  **s 

GROUND  (ING) 

ECD  Cycle  Envelope 

T,  ms  tc,  ms  Tc,  ms 

WAVE  (CW), 

SHG 

tc(cw) 

ING 

0 

0.68 

0.00 

0.68 

0.00 

0.68 

0.00 

0.1 

-4.63 

4.61 

-0.02 

4.61 

-0.02 

4.69 

1 

-1.76 

1.09 

-0.67 

0.98 

2.9921 

0.11 

0.37 

0.48 

-0.02 

0.53 

0.51 

0.43 

0.43 

9.9736 

0.21 

0.48 

0.70 

0.40 

0.38 

0.78 

0.38 

0.28 

40.8917 

0.46 

0.73 

1.19 

0.40 

0.90 

1.30 

0.63 

0.80 

99.7359 

0.42 

1.14 

1.56 

0.48 

1.21 

1.69 

1.05 

1.11 

118.6857 

0.42 

1.24 

1.66 

0.39 

1.67 

1.66 

1.16 

1.52 

135.6408 

0.42 

1.33 

1.75 

0.49 

1.64 

2.13 

1.24 

1.53 

200. 

0.39 

1.65 

2.04 

1.58 

500. 

0.19 

3.00 

3.19 

2.97 

1000. 
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It  is  perhaps  important  to  note  in  Tables  8  and  9  that  the 
constant  Ce=  2 . 904( 10“ 5 )  can  be  arbitrarily  changed  to  tag  other  points 
on  the  pulse  envelope.  But  as  the  propagation  medium  is  changed,  the 
constant  should  indeed  be  kept  constant  so  that  the  envelope  dispersion 
of  a  tagged  point  in  time  can  be  studied. 

It  is  at  first  inspection  of  Table  10  rather  perplexing  to 
note  that  the  secondary  phase  correction  for  the  envelope,  T  ,  decreases 
from  a  value  of  -0.02  us  at  0.1  km  to  a  value  of  -0.67  ys  at  1  km.  Thus, 
equation  (29  a) 

ua  -  u0  =  3.003097  -  0.65  =  2.350971  ys. 
e2  ei 

The  velocity  of  the  envelope  Vg, 

V  =  =  3.828205(10®)  =  1.276151c,  m/s. 

e 

Thus,  if  clocks  were  set  using  the  envelope  tagged  point  in  time  at 
di  =  0.1  km  and  d2  =  1  km,  less  3.003097  ys,  the  clock  set  at  d2  would 
read  earlier  in  time  by  0.65  ys  than  the  clock  set  at  di,  when  the  two 
clocks  are  brought  together  and  compared!  The  apparent  speed  with  which 
the  envelope  traveled  was  1.276951c,  where  c  is  the  speed  of  light.  The 
paradox  lies  in  the  propagation  mechanism  in  the  near  field.  Thus,  in 
this  region  the  attenuation  of  the  higher  frequencies  by  the  ground  is 
nil.  Thus,  there  is  a  linear  increase  or  enhancement  in  radiation 
efficiency  with  frequency  at  very  short  distances  such  that  the  leading 
edge  of  the  pulse  is  enhanced  (rises  faster  as  the  distance  is  increased) 
up  to  the  point  at  which  the  electrical  constants  of  the  ground  attenuate 
the  higher  frequencies  exponentially.  Reference  (7).  Thus,  causality 
is  not  violated  because  t'  >  0  as  specified  at  the  onset  of  this  analysis, 
equations  (19)  and  (9). 

This  discussion  has  assumed  that  the  loran  receiver  has  used 
an  antenna  sensitive  to  the  vertical  electric  field,  £r,  discussed  in 
Appendix  I,  equation  (18  AI).  Loop  antennas  are  also  used  in  Loran-C 
to  replace  the  vertical  whip  antenna.  In  the  far  field,  the  secondary 
phase  correction  for  the  horizontal  magnetic,  H^,  field,  equation  (19  AI) 
is  almost  identical  with  the  vertical  electric  field,  £r.  However, 
the  near  field  is  somewhat  altered.  Thus,  for  example,  the  near  field 
for  equation  (39)  can  be  rewritten  using  equation  (19  AI)  for  this 
field, 
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tf(s,d)*f$(s)  =  go_1C{  /So  exp  (s2a)  erfc(s/x) 

+  to (s+v)  i  •  (57) 

It  is  thus  a  simple  matter  to  rewrite  equation  (49)  leaving 
out  the  last  term  of  equation  (49)  before  performing  the  integration 
specified  by  the  inversion  integral.  The  general  computation  scheme 
described  heretofore  has  also  been  developed  to  study  the  effects  of 
loop  antennas  as  well  as  antennas  receptive  to  vertical  electric 
polarization. 

The  discussion  above  also  assumed  the  electrical  characteris¬ 
tics  of  the  ground  were  constant  with  time.  These  electrical  constants 
may  however  change  with  time  due  to  climatic  changes,  weather  and 
especially  rain  fall.  It  is  evident  that  estimates  of  such  changes 
can  be  used  in  the  theoretical  simulation  discussed  heretofore  to 
account  for  any  such  changes. 


Conclusions  and  Recommendations 


4. 


A  full  transient  propagation  computer  simulation  scheme  has 
been  assembled  and  exercised  for  a  limited  number  of  computations  to 
demonstrate  some  principles  of  transient  behavior.  The  authors  would 
recommend  computation  of  additional  examples  leading  to  a  test  of  the 
validity  of  the  theory  presented. 

It  is  recommended  that  an  enhanced  soil,  hydrology  and 
geologic  data  base  be  developed  in  digital  form  suitable  for  use  in 
computers  for  the  Western  United  States,  especially  in  the  region  of 
Searchlight,  Nevada  and  Death  Valley  California. 

It  is  concluded  that  the  computer  simulation  computations 
presented  in  this  paper  give  results  of  significance  to  the  precision 
prediction  and  interpretation  of  the  secondary  phase  correction  and 
the  envelope  to  cycle  difference.  The  dependence  of  the  envelope  to 
cycle  difference  upon  the  spatial  dependence  of  secondary  phase  correc¬ 
tion  is  weak  in  most  cases  and  automatic  receiver  compensation  for 
spatially  distributed  propagation  errors  should  be  employed  in  the 
context  of  theory  to  guarantee  uniqueness,  if  such  compensation  is 
attempted  indeed. 

It  is  evident  that  the  spatially  distributed  propagation 
parameters  can  vary  with  time  due  to  such  matters  as  rain  fall  and 
weather  conditions  affecting  the  ground  electrical  characteristics. 
These  effects  can  be  studied  numerically  with  the  aid  of  the  analysis 
developed  for  this  paper. 
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6. 


APPENDIX  I:  The  Induction  and  Electrostatic  Fields 


Equation  (39)  for  the  Sommerfeld-Norton  surface  wave  represents 
a  rigorous  solution  to  the  propagation  of  a  continuous  wave  from  vertically 
polarized  electrical  point  source  on  the  surface  of  a  plane  ground, 

e  (.,d)  *  (“  *  jir  *  T3SE7P-]  •  (i  «> 

where  the  attenuation  function 

W  ■  -  j  /apt  exp  (-io2pi)  erfc  (j/pl)j  ,  (2  AI) 

and  the  numerical  distance,  px  =  to2a  is  given  by  equation  (45).  Equation 
(48)  was  derived  from  this  rigorous  formula  by  convolution  with  the  trans¬ 
form  of  the  U.S.  Coast  Guard  standard  pulse, 

f$(«)  =  pi  +  j  (u>c  +  u»)  J  (3  AI) 

The  attenuation  function,  W,  is  related  to  the  scalar  potential  n 
which  is  a  solution  to  the  homogeneous  wave  equation: 

(V2  +  K?)  n  =  0  .  (4  AI) 

The  value  of  W  can  be  calculated  from: 

W  =  n/(2llpri)  ,  (5  AI) 

where,  for  vertical  electric  polarization, 

s  npri  =  IotUoC(4TrbjuD)"1  exp  [j(u>t  -  kxd)J,  (6  AI) 

and,  in  a  spherical  coordinate  system,  r,  6,  4>» 
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D  ■  r2  +  b2  -  2br  cos  0, 
b  =  a  +  hx, 

r  =  a  +  h2, 
a  «  6.36739(10*)  m, 
c  =  2.997925  (10*)  m/s, 
e  =  d/a, 

l0l  =  dipole  current  moment.  Ampere-meters,  which  can  be 
taken  as  ±1,  to  simulate  Loran-C  phase  coding, 

and  hi, 2  are  the  altitudes  of  the  transmitter  and  receiver  above  a  spherical 
(earth)  datum  surface.  As  the  geodesic  distance  d  -*•  0, 


W  -  1  , 


(7  AI) 


n  -f  eo  .  (8  AI) 

Thus,  if  the  antenna  current  is  specified,  and  the  distance,  d, 
is  sufficiently  small,  so  that  (2  AI)  applies,  the  significance  of  (7  AI) 
lies  in  the  second  term  of  equation  (2  AI).  This  is  the  term  containing 
the  ground  electrical  constants  in  the  parameter  px.  Thus,  equation  (7  AI) 
implies  that  the  near  field  is  independent  of  the  ground  electrical  constants. 

In  the  spherical,  polar  coordinate  system,  r,  0,  <f>,  ^|  =  0,  there 
exists  three  TM-waves  (tranverse  magnetic)  that  can  be  calculated 
from  the  potential,  n  =  n(e)  W: 


«r--  <r  s1n  0)"‘  af  [s1"e3ln]  • 

^9  =  r_*  gear  [r  n]  ’ 

V*;  it  ("I  • 


(9  AI) 


C  =  (4Trd)~  Io^PoC  exp  (jwt)  =  29.97925  exp  (jut)  Volts. 


Then, 


II  =  Cs( - jk  1  Db ) ~ 1  Wd  exp(-jkiD),  Volt-meters, 


?r  "  -Cs 


cose 
r  sin0 


—  n  +  — 

30  r 


(10  AI) 
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.  exp  (-ikiP) 


M_1  |  exp  (-jkjD)l  fbr  sin  elf.  1  1 

36*  b  [  D  J  [~B - ■J[1+TOTJ  • 

E  •  C  h  .1  •  [exP.  (-ikiD)~i  .  [  br_sin_e]  _1 

♦  j  WojwBj  D  j  j  D  j  j  1  +  jlcIU 

5r  =  jkl  C  [gxP  (-IktD)  |^2_cps_e  +  3_b_rjSin2  ej 

^  T  2  cos  e  kbr  sin2e  .  -3br  sin2  ell 
+  j  [  k|Da  +  D2 -  +  - J[  Vo1ts/m- 

-*  (W  * 


Cr  .  Jkl  c  |g-xP  Hk»D?  /  z£c 0S_e  +  3_b 

L  J  IL 


For  W  f  1, 


Er  -  Cs(rffA  [*>  si  "<e)  *  "<e)  si  W] 

”  r[(sl  wj  (  sl  n*6>)  *  “  sF  n<C)  +  [sl  n<e’)  [j|« 


+  n(e)  3§4-  W  >  •  Volts/m. , 


■  k  2 1  f  *» 

a.  =cs  ^  “  sl  :n<e>  *  "<<!)  sl  "j  toperes/m  . 


Now.  the  residues  are  calculated  and  summed  In  a  series 
W  A$exp  ^e]-  I AsexpjjiniBs0  exp  fRe&sej 


where. 


A$  A  normalization  function 
6.  *  Reg.  +  j  Img  , 


j|w  -Bsexp[Bse], 

!f  w  =  es  exP  Lese]  • 
es  *  -j  (kja)  »■  a  *■  T$  , 

a  '''.SS  ,  ta'v  .7,  Cs  s  0) 


Thus,  ts  are  the  roots  of 


where 


-  2y2  t  tJ?0 


y  *  -1(k1a)~3'  x"1 


(17  AI) 


For  short  distances,  on  the  ground. 


and 


D  -*•  d,  the  geodesic 
sin  e  -*■  e 

COS  0  1 

0  =  d/a 
hi  =  h2  =  0 


Cr  =  JtoC  1  +  (jwtof1  +  (joito)'2 


W, 


=  J“*o  -C  |1  +  (jwto) 


[l  +  (jwto)_1J  W. 


(18  AI) 
(19  AI) 


Apparently  the  field,  £r  In  the  spherical  coordinate  system  emerges 
with  the  attenuation  function  W  as  multiplicative  factor  to  the  quantity 
(20  AI),  F1, 

Fi  =  [l  +  (juto)-1  +  (juto)'2]  ,  (20  AI) 

while  It  Is  an  additive  term  In  equation  (1  AI)  due  to  Norton  (1937), 

Reference  (11).  This  is  apparently  of  no  numerical  consequence  In  Norton's 
paper,  but  it  somewhat  alters  our  equation  (42)  and  (49),  which  we  derived 
from  Norton's  surface  wave. 
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Upon  rewriting  equation  (46),  for  the  U.S.  Coast  Guard  pulse, 
according  to  equation  (18  AI), 

C(s,d)*  f$(s)  =  biC  +  tfl(s+v)*  +  t02s(s+v)3j 

•  -  Si/ira  exp(s2a)  erfc  (s/5)j  .  (21 

This  results  in  a  probably  more  accurate  but  more  complicated 
expression  for  the  field  than  is  given  in  equation  (49)  In  an  Interned late 
form  It  is  as  follows: 

E(t',d)  =  bivCexp(-vt')  ^  exp  (v2a) 


(21  AI) 


•{[('?  +  2t^  +  7t^ 


t' 

•  j  x  exp  ^-(-^  -  v^a)2j  dx 
0 

t' 

+  jj1  ‘  t^v  “  tvV]  t'~  (  v  +  t7vTjl  *  J  x2exp  ["(27a  “V,^)2J  dT 
L  -■  0 

t'  1 

+  ["  \  +  2t^  +  2t7^]  *  |  t3  exP  •  w«)j  d4 

t'  ° 

+  blVC  [t^V]  ‘  j  exp  [-t2/4“]  ^ 


Again,  after  considerable  ado,  we  find,  compared  to  Equa¬ 
tion  (49),  the  final  form  to  be  the  following: 


E(t',d)  «  bxvCexp(-vt')  |j-  £(1  -  ^ 

+  (v  +t^+^^{1"t^"t^Jr})t'  +  (T^ 

-  4owS  (i  ♦  t^r]  -  {&»  ♦  SaVj  {l  -  ^  -  *j^}| 

2ov,S  (l  -  T&})*'*  (w  *4awS  1+w] 

-  /&>  ♦  8a  v]  h  -  fL  *  t4^)|  •  exp(v2a)  exp  |-  <j£  -  WS)1] 

♦[I'"55  {>-*-  tA})  {$  w} 

*  [l  *  8a2v2}  {l  -  &  -  )  t'  ♦  (»B  [r&j 

-  jl  +  8a2v^|  -  j6aWS  +  8  o*v*^S^ 

*{  1  '  tjv  ‘  tJ*v})]  ^  exp(v2a)  [erfc(-vyS)-erfc(2^-  -  vvS)]j 
+  bivc|  ^4^-  Jl  -  exp(-t'2/4a)jj  .  (23  AI) 


As  we  approach  the  transmitter  and  study  the  near  field  both 
equations  (49)  and  (23  AI)  become  more  accurate  and  we  have  derived  these 
equations  as  special  and  alternative  near  field  representations. 
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APPEFOIX  II:  Oscillograms  of  Loran-C  Pulse 


Figures  2  All  through  11  All  give  oscillograms  automatically 
graphed  by  the  computer  for  the  cases  discussed  In  the  text.  The  cycle. 
Re  E(t%d)  and  the  envelope,  ±|E(t',d)|  of  the  analytic  function  calcula¬ 
ted  by  the  computer,  E(t',d),  are  both  shown.  Magnified  sections  of  the 
pulse  are  also  shown  in  certain  cases.  The  very  special  distances  for 
which  transforms  were  available  turned  out  to  be  non-integer  distances 
when  the  propagation  path  to  Tecopa,  California,  from  Searchlight, 

Nevada  was  divided  Into  uniform  increments  of  distance. 
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FIGURE  1 A  II:  Theoretically  Propagated  Loran-C  Pulse  at  a  Distance  of  2.9921  m  from  the 
Current  Pulse  Given  in  Figure  2,  assuming  the  ground  to  be  smooth,  homogeni 
Impedance,  Ixl  -  .053.  A  magnified  section  of  the  pulse  is  also  shown, 


Re  E(t",d) 


Source  Current 


Theoretically  Propagated  Loran-C  Pulse  at  a  Distance  of  9,974  km  from  the  A 
Source  Current  Given  in  Figure  2,  in  the  Presence  of  tic  Irregular  and  Non- 


FIGURE  5  A  II:  Theoretically  Propagated  Loran-C  Pulse  at  a  distance  of  *10,8917  i 
Current  Pulse  given  in  Figure  2,  assuming  the  ground  to  be  smoot 
Impedance,  1x1  =  .053,  A  Magnified  Section  of  tic  Pulse  is  also 


|E(t',d)|  Re  E(t',d) 


AMPLITUDE,  VOLTS  PER  METER 


Theoretically  Propagated  Ldran-C  Pulse  at  a 
Current  Pulse  Given  in  Figure  2,  Assuming  ti 
Impedance,  Ixl  *  ,055,  A  magnified  section 


I(pM)3|  (pM)3  aa 


II:  Theoretically  Propagated  Loran-C  Pulse  at  a  distance  of  118, 6857  m  from  the  Antenna  Source 
Current  Rjlse  Given  in  Figure  2,  Assuming  the  Ground  to  be  smooth,,  homogeneous  with  an 
Impedance,  Ixl  =  .053,  A  magnified  section  of  the  pulse  is  also  shown, 


(P‘-T)3  9M  |(P‘.T)3| 
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APPENDIX  III:  Topographic  Data  Base 


A  topographic  data  base  has  been  introduced  Into  the  computa¬ 
tion  scheme  for  the  continental  United  States.  This  data  base  Is  Important 
because  It  enters  into  the  estimates  of  ground  Impedance  and  is  also  used 
In  the  ground  wave  propagation  simulator  to  establish  the  phase  perturba¬ 
tions  due  to  the  irregular  nature  of  the  ground.  The  effects  of  terrain  on 
wave  propagation  are  obvious  and  well  established  in  previously  published 
works.  References  (8,  9,  17,  22).  However,  the  effect  upon  ground  impedance 
Is  more  subtle,  depending  upon  the  nature  of  the  soil  overburden  and  its 
relationship  to  the  steepness  of  the  terrain.  This  has  been  briefly  dis¬ 
cussed  in  Reference  (20)  and  is  indeed  a  subject  that  must  be  explored  in 
considerable  depth  for  any  particular  region  before  an  attempt  at  unique 
prediction  is  undertaken. 

To  illustrate  the  importance  of  geophysical  data  bases,  several 
elevation  profiles  are  displayed  herein.  A  sample  of  the  terrain  data  base 
is  displayed  in  Figure  1  A  III.  The  solid  lines  in  the  figure  displays 
ground  elevation  as  a  function  of  distance  along  the  geodesic  between 
Searchlight,  Nevada  and  Ft.  Cronkhite,  California.  In  this  plot,  ground 
elevation  relative  to  sea  level  is  given  in  meters  while  distance  along 
the  geodesic  Is  given  in  kilometers  and  covers  a  total  range  of  200  kilo¬ 
meters.  This  plot  extends  just  beyond  Death  Valley,  California. 

The  following  coordinates  define  the  points  under  consideration: 

Location  Latitude  Longitude 

Searchlight  35°19'  18.18"  N  114°48'  17.43"W  (transmitter) 

Tecopa  35°49‘  4.975"N  116°11 *  5.193"W 

Death  Valley  36°03‘  32.448"N  116°50'  19.467"W 

Ft.  Cronkhite  37°50‘  29.464"N  122°52'  41.170"W 

The  elevation  as  a  function  of  distance  between  Searchlight  and 
Tecopa  Is  shown  as  a  short  dashed  line.  Similarly,  the  elevation  between 
Searchlight  and  the  Death  Valley  reference  point  is  shown  by  the  dash-dot 
line.  Although  the  Tecopa  and  Death  Valley  points  are  approximately  on 
the  geodetic  between  Searchlight  and  Ft.  Cronkhite,  they  are  not  exactly 


so  because  the  terrain  differs  considerably  between  these  curves.  The 
transmitter  datum  and  the  sea  level  datum  are  also  shown.  It  is  of 
Interest  to  note  that  the  region  of  Death  Valley  is  below  sea  level. 

Figure  2  A  III  depicts  elevation  relative  to  sea  level  as  a 
function  of  distance  along  a  geodesic  from  the  Master  station  at  Fallon, 
Nevada  39°33'  06.62"N,  118°49M  56.37"W  to  Searchlight,  Nevada  and  Tecopa 
and  Death  Valley,  California.  Since  these  paths  are  not  along  a  similar 
geodesic  considerable  difference  Is  apparent  between  them.  In  the  Loran-C 
concept  of  using  time  difference  measurements  to  determine  location,  the 
path  to  the  Master  station  is  equally  as  Important  as  the  path  to  the 
secondary  station. 


FIGURE  1  A  III:  Elevation  profiles  along  geodesics  from  Search,  i  gw,  Nevada  to  Tecopa. 
Death  Valley  and  Ft.  Cronkhite,  California. 
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FOREWORD 


This  paper  comprises  Task  2:  Report  on  the  design  of  an  experi¬ 
ment  to  test  the  transient  solution  for  the  U.  S.  Coast  Guard,  Contract  No. 
D0T-CG-842923A.  The  intent  of  this  report  is  to  describe  a  measurement  pro¬ 
gram  and  to  provide  a  basis  for  use  by  the  U.  S.  Coast  Guard  in  establish¬ 
ment  of  definitive  program  goals.  This  experimental  design  utilizes  work 
previously  performed  on  Coast  Guard  Contract  No.  D0T-CG-74629-A  to  review 
Integral  Equation  Computer  Program  for  Loran-C  propagation  and  other  work 
performed  on  this  present  contract,  namely  the  evaluation  of  the  Loran-C 
pulse  dispersion  for  which  a  computer  program  simulation  of  the  ground 
wave  pulse  has  been  configured. 
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ABSTRACT 


The  essential  elements  of  this  combined  theoretical  experimental 
program  include  the  following  items  in  the  time  sequence  that  would  be 
followed: 

(1)  Model  all  baseline  and  transmitter  to  monitor  paths  for 
terrain  and  geophysical  content  of  these  paths; 

(2)  Predict  with  the  propagation  transient  solution  the  ECDs 
and  phase  corrections  that  should  be  observed  locally 
and  remotely  in  the  vicinity  of  each  transmitter  and 
specifically  on  the  baseline  extensions  of  each  trans¬ 
mitter; 

(3)  Measure  both  the  pulse  envelope  and  the  pulse  cycle  propa¬ 
gation  time  on  the  baseline  extensions  for  both  the  local 
and  all  remote  signals; 

(4)  Compare  these  spatially  measured  quantities  with  the 
predicted  values  and  make  adjustments  in  the  modeled 
paths  necessary  to  establish  spatial  fiducials  for  each 
path.  Identify  pulse  envelope  and  cycle  propagation  times 
with  the  envelope  to  cycle  difference  or  ECD  and  the 
secondary  phase  correction; 

(5)  Utilize  this  fiducial  information  as  an  origin  for 
predicting  changes  due  to  temporal  variations; 

(6)  Compare  the  predicted  temporal  variations  with  long 
period  envelope  and  cycle  propagation  time  measurements 
made  at  selected  transmitter  sites;  and 

(7)  Evaluate  the  relationship  between  the  temporal  envelope 
propagation  time  variations  and  the  temporal  cycle 
propagation  time. 

An  addendum  to  this  report  outlines  a  modest  program  using  existing 
data  that  could  be  carried  out  to  demonstrate  the  feasibility  for  embarking  on 
the  above  outlined  program.  Although  existing  data  is  not  capable  of  verifying 
the  ECD  portion  of  the  program,  the  secondary  phase  correction  verification 
effort  is  presumed  to  be  the  most  fundamental  aspect  of  this  entire  program. 
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AN  EXPERIf'ENTAL  PROGRAM  TO  MEASURE  ECD 


AND 

TEST  THE  TRANSIENT  SOLUTION 


R.  H.  Doherty 
A.  R.  Cook 

I .  SUT1ARY  OF  EXPERIMENTAL  PROGRAM 

Baseline  extension  measurements  (BLE)  are  an  important  component 
of  the  calibration  procedure  for  Loran-C  chains.  Thus,  such  measurements 
are  a  chain  management  tool  used  both  in  the  initial  calibration  of  the 
chain  and  in  the  later  maintenance  of  the  chain  calibration.  These  measure¬ 
ments  are  used  as  a  criteria  to  establish  transmitter  emission  delays 
necessary  for  reliable  service  area  coverage.  In  many  cases  heretofore  and, 
particularly,  those  cases  in  which  the  service  area  is  emphasized  as  the  prime 
Loran-C  environment,  calibration  procedures  are  relatively  simple  and  straight¬ 
forward.  In  these  cases,  propagation  path  impedance  estimates  can  be  uni¬ 
form,  atmospheric  refractive  index  effects  can  be  minimal  and  smooth  earth 
propagation  theories  are  often  applied  with  confidence.  When  the  service 
area  involves  land  masses  or  inhomogeneous  propagation  paths  the  problem  is 
much  more  complicated.  Thus,  the  loran  signal  is  propagated  over  ground  that 
exhibits  an  electrical  impedance  that  varies  with  distance  from  the  transmitter 
and  indeed  is  subject  to  geology,  hydrology,  terrain,  and  other  geophysical 
parameter  variations  of  the  earth  surface.  These  applications  require  extreme 
care  in  the  determination  of  the  operational  characteristics  of  a  loran 
chain.  This  care  can  be  ensured  in  the  context  of  an  appropriate  theory  of 
propagation  of  the  groundwave  to  determine  the  adequacy  of  baseline  extension 
measurements  of  phase  and  amplitude  over  the  leading  edge  of  the  pulse.  Thus, 
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both  envelope  and  cycle  propagation  time  (the  secondary  phase  correction  and 
the  ECD)  can  be  observed  simultaneously  and  compared  with  theoretical  predic¬ 
tions. 

With  establishment  of  Loran-C  as  a  primary  navigation  system  for 
general  purpose  usage,  stringent  chain  calibration  requirements  will  become 
Increasingly  Important.  To  support  the  U.S.  Coast  Guard  (USCG)  mission  to 
provide  for  Loran-C  as  a  reliable  general  purpose  radio  navigation  system, 
the  following  effort  should  be  undertaken  to  provide  theoretical  context  and 
validation  for  baseline  extension  and  service  area  measurements  in  the  pre¬ 
sence  of  nonhomogeneous  and  irregular  ground.  Additionally,  the  problem 
should  demonstrate  the  feasibility  of  extending  the  methods  for  the  treat¬ 
ment  of  all  service  area  data.  The  purpose  of  the  program  is  to  develop 
methods  to  ameliorate  calibration  accuracies  and  to  minimize  the  necessity 
for  data  collection  in  the  chain  service  area.  The  program  tasks  should 
establish  precise  relationships  between  baseline  extension  measurements  and 
service  area  calibrations  and  produce  criteria  and  guidelines  that  could  be 
used  for  the  collection  of  pertinent  calibration  data  for  both  secondary 
phase  corrections  and  ECD  variations. 

Accordingly,  this  experimental  plan  to  measure  secondary  phase 
corrections  and  ECDs  and  to  resolve  Loran-C  measurement  inconsistencies 
between  baseline  extension  and  service  area  calibration  methods  is  hereby 
submitted  to  the  U.S.  Coast  Guard.  The  pertinent  features  of  the  program 
are  the  following  items: 

Item  1.  Model  all  West  Coast  Loran-C  chain  baseline  paths  and 
paths  from  the  transmitters  to  the  area  monitors.  The 
modeling  should  include  determination  of  elevations 
and  geophysical  parameters  for  all  points  along  each 
path  (e.g.,  see  Reference  7  and  Addendum  to  this  report). 

Item  2.  Utilize  propagation  predictions  and  solutions  to  the 
transient  analysis  to  predict  envelope  and  cycle  pro¬ 
pagation  time  (secondary  phase  corrections  and  ECDs) 
for  all  of  these  paths.  Even  the  transmitted  ECD  that 
would  be  measured  close  to  the  transmitting  stations 
should  be  predicted  from  the  fact  that  ECD  is  held 
constant  at  the  monitor  and  something  other  than  0 
ECD  is  necessary  at  the  transmitter  to  yieid  this 
condition. 
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Item  3.  Document  and  report  to  the  U.S.  Coast  Guard  all  secondary 
phase  correction  and  ECD  predictions  for  all  close  in  and 
remote  transmitting  stations  for  each  transmitter  and 
each  area  monitor  on  the  West  Coast  Loran-C  chain  before 
measurements  of  ECD  are  made  on  the  baseline  extensions. 

Item  4.  Equip  a  mobile  unit  with  suitably  modified  equipment 

appropriate  to  perform  reliable  master-secondary  and  secon¬ 
dary-secondary  ECD  measurements  as  well  as  secondary  phase 
correction  measurements  on  the  baseline  extensions  of  all 
of  the  transmitting  stations. 

Item  5.  Conduct  a  measurement  program  to  validate  the  theoretical 
development  of  Items  1  and  2  to  verify  the  predictions 
produced  by  Item  2. 

Item  6.  Reevaluate  the  modeling  of  Item  1  and  the  predictions  of 
Item  2  as  a  result  of  any  differences  between  the  pre¬ 
dictions  and  the  measurements.  Modify  the  path  model  or 
the  prediction  algorithms  until  measurements  and  pre¬ 
dictions  agree  within  the  limits  dictated  by  possible 
temporal  changes  that  might  occur  during  the  measurement 
period.  Use  the  modified  path  models  for  fiducial  input 
requirements  for  use  in  Item  7. 

Item  7.  Utilize  the  propagation  programs  described  in  Item  2  along 
with  potential  temporal  parameter  inputs  to  predict  the 
possible  effect  of  temporal  variations  in  both  envelope 
and  cycle  propagation  time  (ECD  and  secondary  phase 
correction). 

Item  8.  Compare  long-term  measurements  of  both  ECD  and  secondary 
phase  changes  made  at  selected  transmitter  sites  with 
variations  predicted  in  Item  7.  Wherever  possible  actual 
weather  changes  should  be  used  to  compare  predicted  effects 
to  measured  variations  in  both  ECD  and  secondary  phase. 

Item  9.  Evaluate  relationships  between  temporal  changes  in  ECD 

and  temporal  change  in  secondary  phase.  Demonstrate  these 
relationships  from  both  the  theoretical  and  the  experimental 


viewpoint.  Thoroughly  document  how  relationships  of  this 
nature  can  be  used  to  improve  the  operational  accuracy  of 
the  Loran-C  system.  Also  indicate  the  necessary  changes, 
if  any,  that  would  be  required  in  the  present  chain  opera¬ 
ting  procedure.  This  study  should  also  demonstrate  the 
area  monitor  function  relative  to  temporal  variations  in 
both  ECD  and  secondary  phase.  Also  the  effect  of  spatially 
distributed  area  monitors  should  be  addressed  from  the 
temporal  variations  standpoint. 


II.  INTRODUCTION  MC  BACKGROUND 


During  the  past  year,  Colorado  Research  and  Prediction  Laboratory, 
Inc.  (CRPL-j)  studied  (under  Contract  DOT-CG-74629-A)  Loran-C  Signal  Propaga¬ 
tion  Prediction  Techniques.  This  effort  resulted  in  the  establishment  of 
improved  Loran-C  pulse  propagation  error  prediction  techniques  for  non- 
homogeneous  irregular  ground  parameters  using  integral  equation  methodologies 
(hereafter  referred  to  as  a  propagation  simulator  that  uses  a  full  wave 
transient  solution  for  the  propagation  of  a  pulse  of  most  any  shape).  In 
addition  to  this  most  recent  accomplishment,  CRPL-j  personnel  have,  in  the 
past,  provided  various  studies  and  reports  concerning  the  Loran-C  propagation 
prediction  problems. 

This  experimental  plan  is  concerned  with  utilizing  technology  unique 
in  prediction  capability  to  address  the  problem  of  the  relationship  between 
calibration  and  ECD  measurements  on  the  baseline  extension  versus  the  service 
area.  In  such  a  program,  the  effects  of  nonhomogeneous  and  irregular  ground 
and  the  modeling  of  the  radio  wave  propagation  mechanism  on  a  geophysical 
data  base  are  of  paramount  importance.  Recognizing  that  the  Loran-C  signal 
is  a  transient  or  pulse  and  not  a  CW  or  continuous  wave  signal,  this  study 
will  evaluate  the  concept  of  "Calibration  by  Prediction"  using  the  full 
pulse  wave  transient  solution. 

The  experimental  plan  is  oriented  toward  measurements  from  the 
West  Coast  Loran-C  Chain  (See  Figure  1).  This  orientation  is  used  because 
the  existing  U.S.  Coast  Guard  transient  analysis  program  under  Contract  DOT- 
CG-842923-A  involves  West  Coast  chain  paths  and  recent  analysis  has  suggested 
some  anomalies  in  this  chain  (Reference  1).  Also  orienting  the  measurement 
program  toward  a  specific  chain  and  specific  locations  clarifies  the  ex¬ 
perimental  descriptions.  Actually  the  experiment  could  De  conducted  utilizing 
any  Loran-C  chain  and  possibly  the  temporal  variations  would  be  better  illu¬ 
strated  using  transmitters  in  the  East  Coast  or  Great  Lakes  chains. 

Relationships  between  ECD  and  the  secondary  phase  corrections  have 
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recently  been  theoretically  proposed  as  a  simple  way  to  correct  for  temporal 
changes  in  secondary  phase  correction  (Reference  2,  pages  154,  155).  Re¬ 
ference  3  (pages  7,  8,  28  and  30)  describes  measurements  where  these  temporal 
variations  in  ECD  were  used  to  explain  differences  in  the  phase  sampling  on 
two  different  receiver  types.  The  theoretical  evaluation  of  these  relation¬ 
ships  (References  2  and  4)  showed  the  relation  of  ECD  to  secondary  phase 
corrections  was  more  a  function  of  earth  curvature  than  of  surface  impedance 
values.  Since  temporal  variations  associated  with  lapse  rate  (the  rate  of 
change  of  index  of  refraction  with  height)  are  changes  of  the  earth  curva¬ 
ture,  it  may  be  possible  to  find  a  functional  relationship  that  would  allow 
temporal  secondary  phase  corrections  to  be  adjusted  from  measurements  of  the 
temporal  ECD  changes.  It  is  intended  to  evaluate  this  possibility  as  the 
final  phase  of  the  present  experiment. 

Complete  descriptions  of  the  specific  tasks  (See  Figure  2)  are 
provided  in  Section  III  of  this  report.  A  bar  graph  depicting  expected  task 
time  scheduling  is  provided  in  Figure  3  of  Section  IV.  The  total  program  is 
anticipated  to  cover  a  time  period  of  18  months  with  an  estimated  labor 
expenditure  of  approximately  four  and  a  half  man  years  or  54  man  months. 

Based  upon  experience  of  experimental  and  theoretical  loran  signal 
propagation  analysis  and  loran  grid  predictions  with  attendant  measurements, 
the  features  listed  in  the  Summary  have  been  selected  as  criterional  to 
successful  completion  of  the  program.  A  developed  intimacy  with  propagation 
phenomena  is  specifically  necessary  to  accomplish  the  required  work  within 
stringent  schedule  requirements.  Intimate  familiarity  with  the  subject  would 
preclude  a  learning  phase  and  therefore  insure  that  all  of  the  participants 
will  make  meaningful  contributions  to  the  program  within  the  shortest  possible 
time  frame.  The  objective  of  these  efforts  should  be  to  provide  the  U.S. 

Coast  Guard  with  a  timely  and  necessary  theoretical  context  and  appropriate 
experimental  techniques  to  reduce  overall  Loran-C  navigation  errors  and 
provide  a  system  that  expedites  safe  and  efficient  transportation  both  in 
the  coastal  confluence  zone  and  elsewhere  throughout  the  service  area. 

This  report  contains: 

(1)  A  definition  of  the  problem  and  historical  basis  for  the 
tasks  to  be  performed; 


(2)  The  technical  approach  which  should  be  used  to  relate  ECD  and 
secondary  phase  measurements  with  spatial  and  temporal  propa¬ 
gation  variations  and  to  establish  statistically  viable 
criteria  for  elimination  of  measurement  data  while  maintaining 
reliable  accuracy; 

(3)  An  outline  of  a  method  for  correlating  ECD  and  secondary  phase 
measured  on  the  baseline  extensions  with  operational  character- 
1st1c4  of  the  chain  Including  service  area  measurements,  (Note: 
the  baseline  extension  Is  defined  as  a  continuation  of  a 
geodesic  connecting  two  transmitters  and  extending  beyond  the 
transmitter  at  either  end  of  the  geodesic.  See  dashed  lines 

in  Figure  1); 

(4)  A  definitive  work  statement  relating  specific  problems  with 
task  requirements,  task  performance,  and  personnel  require¬ 
ments;  and 

(5)  A  program  management  plan  detailing  the  program  structure  work 
flow  and  personnel  that  should  be  committed  to  the  program. 

This  program  should  be  fully  responsive  to  the  needs  of  the  U.S. 
Coast  Guard  during  the  conduct  of  the  experiment  by  making  available  the 
range  of  Loran-C  experience  of  the  entire  professional  staff  engaged  In  the 
theoretical  and  experimental  program. 


FIGURE  2  FLOW  DIAGRAM  OF  AMTICIPATED  EXPERIItNTAL  PROGRAM 
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III.  EXPERIffNTAL  PROGRAM 


All  references  to  measurement  of  ECD  or  prediction  of  ECD  in  this 
program  are  interpreted  to  mean  prediction  or  measurement  of  24  points  on 
the  leading  edge  of  the  Loran-C  pulse  made  using  a  programmable  computer 
oriented  receiver,  e.g.,  Austrom-5000.  These  24  points  are  the  zero  crossings 
of  the  first  6  cycles  both  positive  and  negative  and  the  amplitude  of  the 
pulse  half  way  between  each  pair  of  zero  crossings.  This  will  give  measure¬ 
ments  at  four  known  points  on  each  cycle  of  the  pulse  from  0  to  60  micro¬ 
seconds  (ps)  on  the  pulse.  These  measurements  would  be  similar  to  measurements 
made  by  Fehlner,  et.  al.  (Reference  5).  R.  L.  Frank  (Reference  6)  discussed 
successive-cycle  phase-distortion  using  these  measurements  and  some  preliminary 
predictions  made  by  Johler,  et.  al.  (Reference  2,  Appendix  II,  Tables  A2  to 
A6).  i.is  discussion  and  analysis  is  very  closely  related  to  ECD  and  secondary 
phase  correction  problems. 

By  measuring  and  predicting  24  points  along  the  leading  edge  of  the 
Loran-C  pulse  a  complete  dispersion  comparison  can  be  made.  Receiver  band 
width  characteristics  can  be  included  in  these  predictions  as  described  in 
Reference  14  (pages  12  and  13).  Also  from  either  the  predictions  or  the 
measurements,  ECDs  as  determined  by  various  receiver  manufacturers  can  be 
calculated.  However,  if  ECDs  were  measured  by  any  one  technique  utilized  by 
a  specific  receiver  manufacturer,  it  is  not  readily  apparent  that  an  ECD  could 
be  predicted  for  a  different  manufacturers  technique. 

Therefore,  in  this  experimental  program  definition,  ECD  measurements 
and  ECD  predictions  will  imply  measurements  or  prediction  of  the  24  specified 
points  on  the  leading  edge  of  the  Loran-C  pulse.  The  secondary  phase  pre¬ 
dictions  should  relate  to  the  CW  solution  of  the  propagation  equations  and 
the  secondary  phase  correction  measurement  is  the  difference  between  the 
signal  propagation  at  the  speed  of  light  in  a  standard  atmosphere  and  the 
measurement  made  at  the  standard  zero  crossing  of  third  cycle  of  the  Loran-C 
pulse. 


Figure  2  is  a  block  diagram  of  the  proposed  experimental  program. 

The  various  blocks  in  this  block  diagram  are  addressed  in  the  following 
tasks  of  this  section. 

III.l  Task  1:  rbDELiNG  the  Paths 

The  principal  prerequisite  to  successful  conduct  of  the  total  pro¬ 
gram  is  an  adequate  modeling  of  appropriate  propagation  paths  in  the  chain 
service  area.  For  these  experiments,  propagation  paths  identified  by  the 
base  lines  are  expected  to  provide  the  maximum  amount  of  meaningful  data  for 
expected  effort.  These  paths  are  depicted  in  Figure  1.  Reference  to  the 
figure  reveals  a  total  of  fourteen  paths  necessary  to  adequately  describe 
the  system  and  validate  the  theoretical  transient  solution.  See  Addendum  to 
this  report  for  more  information  for  some  of  these  paths. 

To  facilitate  the  efficient  use  of  propagation  simulators  and  in 
compliance  with  the  concepts  presented  in  CRPL^'s  Report  78-9  (Reference  7) 
entitled  "Geophysical  and  Geological  Data  Base  Evaluation  for  Loran-C  Ground 
Wave  Propagation  Medium,"  propagation  path  descriptions  will  include  terrain, 
geological  substructure,  and  soil  composition  parameters  including  hydrology. 
These  geophysical  data  parameters  when  convoluted  with  accepted  electro¬ 
magnetic  ground  wave  propagation  theory  are  herein  denoted  as  the  geopara- 
metric  propagation  technique. 

A.  Terrain  Data 

In  the  case  of  terrain,  the  data  necessary  to  provide  adequate  pro¬ 
pagation  path  description  requires  a  resolution  to  account  for  relatively 
gross  irregularities  of  the  earth's  surface.  Previous  studies  have  demon¬ 
strated  that  the  natural  ground  over  which  the  Loran-C  signal  propagates 
tends  to  become  inductive  (References  8  and  9).  Irregularities  in  the 
earth's  surface  and  g*-oss  inhomogeneities  beneath  the  surface  often  reflect 
Inductive  ground  surface  impedances  along  the  propagation  path.  Experience 
demonstrates  that  extremely  small  detail  is  neither  required  technically  or 
desirable  from  an  economic  viewpoint.  In  general,  gross  irregularities  appear 
to  be  adequately  described  with  a  data  resolution  on  the  order  of  one  kilometer 
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spacings  (see  References  7  and  15).  In  those  special  cases  of  extremely 
rugged  terrain,  a  resolution  of  up  to  200  meter  spacing  may  be  required. 

Sources  for  terrain  data  are  available  from  CRPL^  computerized 
Continental  United  States  terrain  data  base  (see  Addendum  to  this  report). 
To  presciently  determine  the  necessary  data  resolution  for  a  particular 
area  requires  experienced  judgment  of  potential  effects  of  terrain  on  the 
signal  propagation.  Based  upon  this  experience,  the  CRPL^  terrain  data 
base  appropriate  for  loran  signal  analysis  has  been  assembled. 

B.  Geological  Subsurface  Structure  Data 


A  second  important  characteristic  of  a  Loran-C  propagation  path 
is  contained  in  the  geological  subsurface  structure.  It  has  long  been  ob¬ 
served  that  broadcase  frequencies  exhibit  attenuation  rate  changes  with  the 
geological  underlayment  that  correlates  with  geologic  age.  Similarily, 
there  is  a  strong  correlation  of  Loran-C  frequency  attenuation  rates  and 
phase  lags  with  the  age  of  the  geologic  subsurface  structure. 

The  purpose  of  this  subtask  is  to  identify  and  assemble  a  data 
base  of  geological  subsurface  parameters  required  to  describe  the  under¬ 
layment  along  all  propagation  paths  relevant  to  this  program.  In  general, 
the  data  will  be  obtained  from  geological  maps  of  the  service  area.  To 
accomplish  this  work  will  require  significant  effort  to  sift  through  a  vast 
amount  of  geological  information  to  obtain  a  small  amount  of  relevant  data. 
The  procedure  will  be  one  of  development  of  an  expertise  that  will  require 
judgment  and  intuition.  The  learning  process  is  ameliorated  when  performed 
in  the  context  of  previous  and  allied  propagation  data  base  assembly  efforts. 
We  at  CRPL-j  have  begun  assembly  of  these  type  of  bases. 

The  geological  underlayment  modeling  for  the  propagation  path  will 
be  a  digitization  of  data  points  along  the  path  at  intervals  comparable  to 
other  data.  The  data  density  and  modeling  technique  at  any  point  is 
dictated  by  expected  effects  of  the  area  on  the  signal  propagation.  For 
example,  areas  that  exhibit  singular  discontinuities  in  geological  under¬ 
layment  are  expected  to  require  special  modeling  techniques  for  resolution. 


Here  again,  previous  experience  in  similar  anomonlies  will  provide  criteria 
for  judgments  for  initial  conductivity  assignment  and  later  modeling  adjust¬ 
ments.  Although  this  subtask  represents  a  considerable  investment  of  effort, 
it  is  expected  to  provide  enormous  benefits  toward  high  prediction  goals  for 
Loran-C. 

C.  Soil,  Soil  Depths,  and  Hydrology 

The  final  major  consideration  for  characterization  of  the  Loran-C 
signal  propagation  path  is  the  overburden  or  soil  parameters.  In  general, 
soil  tends  to  decorrelate  with  the  age  of  bedrock  in  its  formation.  Soil 
thicknesses  are  quite  variable  over  the  surface  of  the  earth  and  can  vary 
from  exposure  of  the  basement  rock  to  depths  of  several  hundred  meters. 

An  electromagnetic  ground  wave  propagating  over  the  earth's  sur¬ 
face  is  influenced  by  the  geoelectrical  properties  of  the  soil.  Waves 
entering  the  ground  are  exponentially  attenuated  by  the  resistivity  of  the 
surrounding  material.  These  electrical  parameters  can  vary  with  the  depth 
or  thickness  of  the  soil  through  which  the  wave  penetrates.  Accordingly, 
for  reliable  description  of  soil  characteristics  indigeneous  to  geographical 
areas  along  propagation  paths  it  is  necessary  to  configure  area  soil  physio¬ 
graphies  in  layers  characterized  by  appropriate  geoelectrical  parameters.  A 
predominate  consideration  in  the  assignment  of  geoelectrical  parameters  to 
defined  soil  strata  is  the  area  hydrology.  Water  tables  rise  and  fall,  and 
the  corresponding  change  in  moisture  content  can  cause  the  parameters  to 
vary.  These  considerations  impact  the  number  of  strata  or  geology  horizons 
necessary  to  adequately  describe  a  geographical  region. 

Soil  data  can  be  retrieved  from  soil  maps  and  other  soil  literature 
available  in  geology  libraries.  Other  sources  may  be  available  in  certain 
government  agencies,  particularly  those  of  the  Department  of  Agriculture. 

Here  again  there  is  an  enormous  amount  of  literature  to  be  sifted  to  obtain 
a  small  amount  of  relevant  and  required  data.  As  before,  experience  in 
assembling  these  types  of  data  bases  is  desirable  to  economically  and  effi¬ 
ciently  perform  this  subtask. 

The  data  will  be  digitized  at  discrete  intervals  along  the  path. 


The  spacing  within  the  Intervals  should  correspond  to  spaclngs  of  the  other 
data  basss.  In  our  experience,  soil  descriptions  can  be  "sketchy"  In  some 
geographical  locations  and,  for  example,  the  assessment  of  soil  depth  diffi¬ 
cult  to  obtain.  In  this  case,  reliable  path  descriptions  require  specialized 
modeling  techniques.  Here  again  experience  can  play  an  Important  part  In 
the  successful  completion  of  the  total  program. 

1 1 1. 2  Task  2:  Spatial  Propagation  Prediction 

A  secondary  phase  correction  prediction  should  be  made  for  each  of 
the  14  modeled  paths  using  a  propagation  simulator  for  nonhomogeneous  irreg¬ 
ular  ground  paths.  This  geoparametric  propagation  technique  should  utilize 
the  following  imputs:  establish  terrain  data,  initial  estimates  of  con¬ 
ductivity  for  soils,  initial  estimates  of  soil  depths  and  hydrology,  and 
initial  estimates  of  conductivity  of  substructure  of  geological  formations 
all  convoluted  to  represent  surface  impedances  along  each  path.  The  resulting 
secondary  phase  corrections  should  be  used  to  derive  time  difference  values 
at  every  location.  The  time  difference  values  will  of  necessity  relate  to 
the  fact  that  the  time  difference  numbers  at  the  area  monitor  are  held  at  a 
particular  value  and  are  held  constant.  On-site  comparisons  between  computed 
time  difference  numbers  and  measured  time  difference  numbers  at  the  various 
locations  should  provide  a  means  for  interpreting  the  validity  of  the  measure¬ 
ment  data.  They  should  also  determine  if  more  measurements  should  be  made, 
or  if  the  measurements  made  at  the  time  are  adequate. 

An  ECD  or  dispersion  prediction  should  also  be  made  for  each  of 
the  14  propagation  paths.  This  prediction  should  be  made  using  techniques 
similar  to  those  used  in  the  U.S.  Coast  Guard  transient  model  program  (Re¬ 
ference  15).  The  actual  predictions  should  yield  each  zero  crossing  and 
each  half  cycle  amplitude  for  the  first  6  cycles  of  the  Loran-C  pulse.  The 
initial  prediction  for  the  area  monitor  should  be  made  using  the  trans¬ 
mitter  pulse  of  known  ECD  as  measured  at  the  transmitting  station.  (That 
is,  that  the  third  cycle  axis  crossing  will  occur  at  a  precisely  known  point 
relative  to  the  third  cycle  crossing  derived  from  the  transient  solution 


utilizing  the  standard  mathematically  specified  pulse).  Subsequently, 
a  prediction  should  be  made  for  the  transmitters  assuming  the  pulse  re¬ 
ceived  at  the  area  monitor  has  the  chain  specified  ECD  value  and  that 
this  pulse  conforms  to  the  standard  mathematically  specified  pulse.  As 
was  pointed  out  by  M.  Dlshal  (Reference  10),  the  pulse  that  would  be 
radiated  by  the  transmitter  could  not  precisely  conform  to  the  mathema¬ 
tically  specified  pulse.  An  analysis  of  the  error  terms  discussed  by 
Dlshal  convoluted  with  the  degree  to  which  a  measured  pulse,  l.e.,  4  points 
per  cycle  for  the  first  6  cycles,  at  both  the  transmitters  and  at  the  area 
monitors  does  not  conform  to  the  mathematically  specified  pulse  will 
demonstrate  the  ability  to  make  predictions  and  measurements  agree.  This 
comparison  might  influence  the  final  operational  recommendations  that  would 
result  from  this  experimental  program. 

The  transient  model  program  Is  capable  of  predicting  propagation 
for  the  mathematically  defined  pulse.  It  will  also  be  capable  of  predicting 
a  pulse  that  would  yield  the  mathematically  specified  pulse  after  traversing 
a  described  propagation  path.  Pulse  shapes  other  than  the  specified  pulse 
could  also  be  predicted  using  the  transient  model  theory  (see  Reference  11, 
for  example).  However,  this  would  require  an  additional  research  effort. 
Furthermore,  predictions  using  the  specified  pulse  should  be  utilized  fully 
to  determine  if  predictions  of  other  pulse  shapes  are  required. 

The  ECD  spatial  predictions  should  also  be  made  over  the  trans¬ 
mitter  baselines.  In  order  to  make  these  predictions  the  pulse  used  for 
the  transmitter  pulse  should  be  the  one  derived  from  the  prediction  assum¬ 
ing  the  specified  pulse  and  the  ECD  values  assigned  for  the  monitor.  A 
comparison  of  these  predictions  for  both  ends  of  the  baseline  paths  with 
measurements  made  at  both  ends  of  the  baseline  paths  will  enhance  confidence 
of  the  prediction  techniques,  or  demonstrate  the  limitations  of  the  assump¬ 
tions  utilized  to  make  the  predictions. 

Finally,  all  of  the  ECD  and  secondary  phase  correction  predictions 
of  the  24  points  on  the  leading  edge  of  the  pulse  should  be  analyzed  to 
produce  ECDs  that  would  be  obtained  by  the  different  receiver's  measuring 
techniques.  In  this  study,  a  careful  evaluation  should  be  made  of  the 
effect  on  the  receivers  due  to  pulses  whose  pulse  shapes  are  not  the 
mathematically  specified  pulse  shape.  If  a  particular  receiver  design 
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appears  to  be  critically  affected  by  the  pulse  shape,  a  recommendation  may 
be  made  to  add  this  receiver  type  to  the  measurement  program  outlined  in 
the  next  section. 

A  plan  to  measure  the  ECD  (that  is  the  24  points  on  the  leading 
edge  of  the  pulse)  using  a  programmable  computer  oriented  receiver  such  as 
the  Austron-5000  should  be  devised  for  the  Area  Monitor.  This  measurement 
would  be  made  on  only  one  transmitter  signal  at  a  time,  and  it  would  be 
made  on  the  transmitter  being  visited  by  the  mobile  van  during  the  period 
the  van  is  at  that  particular  transmitter,  if  possible.  It  is  anticipated 
that  a  plan  might  be  worked  out  using  some  computer  time  on  the  resident 
Austron-5000.  Otherwise,  a  second  Austron-5000  should  be  provided,  if  one 
could  be  made  available  for  these  measurements,  so  that  these  measurements 
would  not  interfere  with  the  area  monitors  normal  operational  program.  These 
measurements  would  only  need  to  be  made  a  few  times  per  day  for  one  or  two 
days  for  each  transmitter  over  a  total  period  of  approximately  three  weeks 
while  the  mobile  van  was  measuring  on  the  various  transmitter  baseline 
extensions.  However,  if  a  reasonably  simple  automated  sequential  program 
for  making  these  measurements  could  be  devised,  it  would  be  desirable  to 
make  these  measurements  at  the  area  monitors  during  the  entire  period  when 
long-term  U.S.  Coast  Guard  observed  ECD  and  secondary  phase  measurements  are 
made  at  the  transmitter  locations  (see  Sections  III. 5  and  III. 6). 

111,3  Task  3:  Spatial  Feasurements  and  Required  Equipment 

This  section  of  the  experimental  plan  provides  for  the  require¬ 
ments  of  a  mobile  research  van  to  make  the  necessary  ECD  and  secondary  phase 
measurements.  Also  addressed  in  this  section  is  the  design  and  implementa¬ 
tion  of  the  spatial  measurement  program. 

A  four-wheel  drive  vehicle  is  required  to  house  the  measurement 
equipment  and  thus  permit  access  to  baseline  extensions  not  otherwise 
accessible  to  two-wheel  drive  vehicles.  Principle  equipment  required  for 
the  mobite  van  includes: 

(1)  A  programmable  computer  oriented  Loran-C  receiver 
programmed  for  this  experiment; 

(2)  Hewlett  Packard  Cesium  Standard  Oscillator; 

(3)  Heavy  duty  Alternator; 

(4)  Extra  battery  capacity. 


The  alternator  and  other  electrical  power  supply  capability  must  be  suffi¬ 
cient  to  operate  the  measurement  equipment  without  interruption.  Accordingly, 
extra  battery  capacity  is  required  and  facilities  for  recharging  these 
batteries  from  commercial  power  sources  is  necessary.  Also,  the  van  should 
be  equiped  to  convert  from  battery  operation  to  commercial  power  operation 
with  no  expected  loss  or  interruption  to  equipment  operation.  The  equipment 
should  be  mounted  in  the  van  in  a  manner  to  facilitate  easy  removal  and, 
also,  configured  in  a  way  that  permits  storage  and  operation  in  the  cargo 
space  of  a  helicopter.  Although  specific  helicopter  applications  are  not 
spelled  out  in  this  experimental  plan,  previous  experiences  with  mobile  van 
type  measurements  have  demonstrated  the  potential  advantages  that  may  occur 
if  the  equipment  is  so  configured.  In  addition  to  the  above  described  equip¬ 
ment,  the  van  should  have  adequate  test  equipment  to  determine  proper  opera¬ 
tion  for  the  Loran-C  measurements.  This  test  equipment  should  be  sufficient 
to  diagnose  problems  and  to  repair  most  equipment  malfunctions. 

The  design  and  development  of  this  mobile  research  van  should  also 
include  equipment  modifications  that  might  be  necessary  or  desirable  to  best 
perform  the  spatial  measurements.  The  functions  of  the  programmable  computer 
oriented  receiver  should  be  understood  and  documented  to  such  a  degree  to 
which  the  measurements  could  be  automotized.  Perhaps  this  program  would  be 
better  served  by  a  receiver  oriented  programmable  computer,  since  the  computer 
capabilities  and  software  development  may  well  be  crutial  to  the  success  of 
the  overall  experiment.  The  techniques  for  making  the  measurements  along  with 
the  determination  of  the  potential  automation  of  the  measurements  should  be 
well  documented  for  use  by  the  U.S.  Coast  Guard  for  making  long-term  temporal 
measurements  at  fixed  transmitter  sites  with  similar  equipment.  These  temporal 
measurements  would  be  needed  for  the  latter  parts  of  the  over-all  experimental 
plan. 

The  actual  measurement  program  for  the  spatial  effects  should  be 
thoroughly  planned.  Primarily  this  program  should  consist  of  making  secondary 
phase  and  ECD  measurements  on  the  baseline  extensions  of  each  transmitter 
as  shown  by  the  dotted  paths  in  Figure  1.  Measurements  of  both  the  close 
and  distant  transmitter  should  be  obtained  at  each  location.  For  the  West 
Coast  chain  there  are  actually  twelve  baseline  extensions  where  these 
measurements  should  be  made.  Careful  analysis  of  the  area  maps  for  the 
vicinity  of  the  transmitting  stations  would  allow  the  measurement  program  to 
be  planned  in  advance.  Several  locations  on  each  baseline  extensions  should 


be  chosen,  where  it  appears  possible,  to  drive  a  four-wheel  drive  vehicle. 

When  the  measurements  are  actually  made,  these  would  be  the  starting  points 
for  the  measurement  program.  The  number  of  sites  that  might  be  occupied  on 
each  baseline  extension  would  depend  upon  insitu  progress  of  the  .measurement 
program.  A  continuous  evaluation  of  measured  data  during  the  data  collec¬ 
tion  effort  is  essential.  Also  very  thorough  records  are  required  during 
the  measurement  program  so  that  suspect  measurements  could  be  repeated  and 
so  that  all  measured  data  is  carefully  documented  as  to  location  and  time 
of  measurement. 

Making  measurements  only  on  the  baseline  extensions  of  the  trans¬ 
mitter  pairs  somewhat  minimizes  the  need  for  precise  location  of  the  measuring 
point,  however  the  best  source  of  local  maps  (7*2  or  15  minute  U$GS  topo¬ 
graphical  maps)  should  be  used  to  locate  the  van  as  accurately  as  possible 
for  each  measurement  point.  All  measured  data  should  be  preliminarily 
evaluated  by  making  comparisons  with  predictions  and  by  using  small  pro¬ 
grammable  calculator  programs  to  check  the  data  as  it  is  being  measured 
similar  to  the  after  the  fact  analysis  described  in  Reference  1. 

Any  major  difference  between  what  is  measured  and  what  was  pre¬ 
dicted  for  a  particular  path  should  provide  incentive  to  investigate  that 
particular  path  or  location  more  thoroughly.  This  type  of  flexibility  and 
on-site  evaluation  in  a  measurement  program  can  often  determine  the  success 
of  the  experimental  program.  Frequently  poorly  obtained  measurements  lead 
to  experimental  programs  that  produce  more  questions  than  they  are  capable 
of  answering. 

Ill Task  Spatial  Verification  and  Modification 

Following  a  comprehensive  spatial  variation  measurement  program, 
the  initial  predictions  must  be  evaluated  and  prediction  parameters  should 
be  modified  where  the  predictions  did  not  meet  expected  accuracy.  Compari¬ 
sons  of  the  measured  and  the  predicted  results  may  be  evaluated  over  the 
baseline  paths  where  measurements  at  both  ends  of  the  paths  allow  a  unique 
deduction  of  both  the  secondary  phase  and  ECD  variations  over  the  path. 
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Parameters  that  can  be  used  to  correctly  predict  these  baseline  paths  must 
also  work  to  predict  the  paths  from  the  transmitters  to  the  area  monitors. 

The  agreement  between  the  spatial  predictions  and  the  spatial 
measurements  can  be  no  better  than  the  stability  of  the  spatial  measurements 
during  the  measurement  period.  That  is,  in  order  to  determine  a  fiducial 
value  for  the  propagation  over  a  particular  path,  the  temporal  variations 
that  may  have  occurred  during  the  measurement  period  must  be  ignored.  There¬ 
fore,  during  the  period  the  baseline  extension  measurements  are  being  made, 
careful  records  should  be  maintained  for  all  corrections  introduced  by  the 
area  monitors.  Also,  rather  comprehensive  weather  records  should  be  main¬ 
tained  for  that  period  so  that  estimates  of  possible  temporal  variations 
can  be  made.  The  secondary  phase  correction  and  the  ECD  effects  over  a 
particular  path  should  be  nearly  one  order  of  magnitude  greater  for  the 
spatial  effects  than  for  the  temporal  effects  (Reference  12).  Therefore, 
there  should  be  no  problem  in  deducing  a  fiducial  value  for  each  path  that 
predicts  the  spatial  variations  and  that  can  be  used  as  the  starting  point 
for  predicting  temporal  variations  over  that  particular  path. 

The  actual  measured  ECD  data  for  each  path  should  be  analyzed 
using  computer  simulation  techniques  equivalent  to  those  used  by  various 
receiver  manufacturers.  Equivalent  ECDs  for  different  receiver  types  will 
indicate  how  variable  the  different  manufacturers  equipment  may  be  for 
actual  pulse  dispersion  changes  that  may  occur. 

Once  allowance  has  been  made  for  the  temporal  variations,  a  com¬ 
parison  between  the  original  propagation  prediction  defined  in  III. 2  and 
the  actual  measured  data  defined  in  III. 3  will  allow  one  to  evaluate  how 
much  the  original  estimates  of  the  soil  or  geologic  conductivities  were  in 
error.  By  simultaneously  examining  both  the  ECD  and  secondary  phase  correc¬ 
tion  errors  for  all  14  paths  it  should  be  possible  to  determine  where  the 
major  errors^/ere  introduced. 

Af^ter  careful  simultaneous  study  of  both  secondary  phase  errors 
and  ECD  errors  for  all  14  paths,  judiciously  chosen  adjustments  of  the 
conductivity  of  certain  soil  types,  certain  geologic  rock  types  or  certain 
soil  depths  or  hydrology  should  be  made  and  repredictions  of  both  secondary 
phase  and  ECDs  for  all  paths  should  be  made.  This  process  should  be 
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repeated  until  agreement  between  predicted  and  measured  data  approaches  the 
limit  allowed  by  temporal  variations.  During  this  analysis  process,  parti¬ 
cular  attention  should  be  attributed  to  the  sensitivity  of  ECD  predictions 
to  the  various  parameters  under  consideration. 

The  final  selection  of  conductivity  values  for  the  various  sub¬ 
surface  geoparameters  will  provide  the  fiducial  information  necessary  as 
Input  information  for  propagation  predictives  of  the  temporal  variations. 
Thus,  the  geoparametrlc  variations  will  permit  a  fit  to  spatially  distributed 
observations  which  determines  the  proper  conductivity  values.  The  primary 
benefit  of  this  entire  effort  will  be  the  validation  of  secondary  phase 
corrections  using  the  geoparametric  parameters  and  integral  equation  solu¬ 
tions.  The  use  of  the  transient  solution  which  includes  the  integral  equa¬ 
tion  allows  a  further  evaluation  of  the  ECD  variations.  If  the  ECD  problem 
is  not  important,  this  study  will  still  be  very  important  for  establishing 
the  methods  for  predicting  the  secondary  phase  corrections. 

I I I. 5  Task  5:  Temporal  Propagation  Predictions 

Utilizing  results  from  a  dynamic  model  study,  the  temporal  varia¬ 
tions  for  the  secondary  phase  changes  can  be  evaluated  for  a  selected  path 
that  could  be  monitored  for  temporal  variations.  This  experimental  plan  does 
not  include  man  hours  or  equipment  for  monitoring  long-term  temporal  varia¬ 
tions.  Recommendations  for  the  use  of  programmable  computer  oriented  re¬ 
ceivers  such  as  Auctron-5000  receivers  at  both  ends  of  the  path  to  be 
monitored  for  temporal  variations  would  be  made.  It  is  anticipated  that 
these  measurements  could  be  made  over  an  existing  baseline  path  and  could 
be  made  by  employing  the  special  receivers  at  two  existing  transmitter  sites. 
These  measurements  of  the  temporal  secondary  phase  changes  and  ECD  changes 
would  be  made  in  accordance  with  recommendations  made  as  a  result  of  the 
spatial  measurement  information  and  instrumentation  design  work  gained  in 
Section  I II. 3. 

The  results  from  the  d  ::iic  model  study  should  be  used  to  predict 
temporal  changes  in  the  secondary  phase  as  a  function  of  estimated  impedance 
changes  due  to  weather  effects,  estimated  changes  in  the  surface  refractive 
index  and  also  estimated  changes  in  the  lapse  rate  of  the  refractive  index 
(that  is,  change  of  refractive  index  with  height  above  the  surface).  These 
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changes  should  be  estimated  based  on  the  types  of  terrain  the  path  traverses, 
l.e.,  desert  conditions,  typical  farm  land  conditions,  mountainous  terrain 
and/or  coastal  peneplains. 

The  above  predictions  of  secondary  phase  corrections  should  be 
extended  to  an  ECO  type  prediction  using  a  transient  propagation  solution 
such  as  that  described  in  the  final  report  of  the  transient  model  U.S.  Coast 
Guard  study  (reference  15).  Here  again,  predictions  should  be  based  upon 
anticipated  changes  in  surface  impedance,  the  surface  refractive  index,  and 
the  refractive  index  lapse  rate. 

All  of  these  temporal  variation  predictions  should  be  made  for 
long  period  durations  and  for  the  paths  where  the  temporal  measurements  are 
scheduled  to  be  made.  For  these  predictions  anticipated  weather  variations 
for  the  particular  chosen  paths  should  be  used  for  a  prediction  basis.  Based 
upon  the  results  of  previous  studies  (Reference  13),  particular  emphasis  on 
the  lapse  rate  predictions  should  be  made  if  any  of  the  selected  paths  are 
subjected  to  below  freezing  temperatures  for  long  periods  of  time. 

All  of  the  temporal  predictions  for  the  baseline  paths  over  which 
the  measurements  would  be  made  must  take  into  account  that  the  system  is 
held  constant  at  the  area  monitor  and  not  at  the  transmitter  sites.  There¬ 
fore,  the  paths  from  the  transmitters  to  the  area  monitors  must  be  predicted 
and  these  predictions  used  to  remove  the  effect  of  the  area  monitor  control 
as  seen  over  the  baselines.  The  true  variations  over  the  baselines  can  be 
deduced  from  measurements  made  at  both  ends  of  the  baseline  paths.  The 
effects  occurring  over  the  paths  from  the  transmitter  to  the  area  monitor 
can  be  evaluated  from  the  difference  between  the  apparent  baseline  varia¬ 
tions  as  seen  from  the  measurements  and  the  true  baseline  variations  as 
deduced  from  a  mathematical  treatment  of  the  measurements  obtained  at  both 
ends  of  the  propagation  path. 

The  effect  of  the  area  monitor  control  process  is  extremely  im¬ 
portant  to  the  interpretation  of  propagation  measurements.  This  subject  has 
been  treated  somewhat  in  References  2  and  3,  but  a  thorough  understanding  of 
all  of  the  ramifications  of  this  effect  is  mandatory  to  make  the  present 
experiment  successful. 

III. 6  Task  6:  Temporal  Validation  and  Edification 

The  objective  here  is  to  compare  measured  temporal  variations 
with  predictions  of  an  estimated  range  of  variations  that  might  be  anticipated. 


-22- 


Here  again  it  is  necessary  to  discuss  the  effect  of  an  area  monitor.  It 
is  necessary  to  see  that  the  effect  on  the  secondary  phase  and  the  effect 
on  the  ECD  is  somewhat  different.  The  area  monitor  maintains  a  constant 
secondary  phase  correction  at  the  area  monitor.  The  value  for  this  correction 
has  been  derived  from  a  system  calibration  involving  measurements  on  a 
number  of  points  throughout  the  service  area.  The  transmitted  signals  are 
adjusted  to  keep  the  area  monitor  time  difference  fixed.  Therefore,  tem¬ 
poral  changes  over  the  path  from  the  Secondary  transmitter  to  the  area 
monitor  minus  the  temporal  changes  over  the  path  from  the  Master  to  the  area 
monitor  must  be  removed  in  order  to  determine  the  actual  radiation  time 
at  the  transmitting  station.  To  a  much  lesser  extent  a  similar  condition 
exists  for  the  ECD,  since  an  assigned  ECD  number  is  held  for  each  transmitter 
at  the  area  monitor.  Since  ECD  values  may  be  held  with  larger  tolerance 
levels  than  for  the  phase  values,  the  ECD  changes  will  be  much  less  frequent. 
Also,  ECD  changes  will  occur  separately  for  the  Master  or  Secondary  path, 
thereby  allowing  separation  of  the  path  changes. 

It  would  be  most  desirable  for  a  temporal  analysis  study  to  be 
able  to  isolate  secondary  phase  adjustments  and  ECD  adjustments  that  are 
propagation  oriented  from  those  that  are  operator  or  equipment  oriented. 

Every  effort  would  be  made  to  design  the  measuring  equipment  in  a  manner 
that  would  help  resolve  this  issue.  This  discourse  is  meant  to  point  out 
the  many  ramifications  that  become  important  in  order  to  try  to  interpret 
propagation  measurements  or  data  that  supposedly  will  yeild  propagation 
information. 

From  measurements  made  at  both  ends  of  a  baseline  path  and  the 
corrections  introduced  by  the  area  monitor,  it  should  be  possible  to 
determine  both  the  ECD  and  secondary  phase  variations  that  occurred  over 
the  baseline  path  and  over  the  paths  from  the  transmitters  to  the  area 
monitor.  These  variations  should  be  compared  with  the  predicted  variations 
for  each  of  the  paths.  The  predicted  temporal  variations  for  changes  in 
surface  impedance  and  changes  in  surface  refractive  index  have  a  different 
ECD  to  secondary  phase  relationship  than  the  predicted  variations  due  to 
lapse  rate  changes.  Therefore,  comparing  the  measured  temporal  variations 
with  the  predicted  temporal  variations  for  these  paths  should  indicate 
whether  the  changes  were  due  to  surface  parameter  changes  or  were  due  to 
lapse  rate  changes. 

A  comprehensive  analysis  of  all  of  the  measured  temporal  variation 
data  for  all  of  the  monitored  paths  should  allow  one  to  provide  a  maximum 
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error  level  estimated  for  the  particular  chain  and  for  the  existing  operating 
procedures  for  that  chain.  This  error  model  could  be  provided  for  all  areas 
served  by  that  particular  chain.  An  error  model  map  of  this  nature  would 
provide  a  means  for  establishing  maximum  errors  that  could  be  expected  at 
any  given  location  in  the  coverage  area  due  to  temporal  variations  in  the 
propagated  signals. 

The  measured  temporal  variations  for  the  leading  edge  of  the  pulse 
(commonly  referred  to  as  ECD  measurements)  should  be  subjected  to  computer 
analysis  simulating  various  receiver  techniques  for  deriving  ECD  from  the 
pulse.  If  different  receiver  techniques  give  different  ECD  changes,  these 
should  be  studied  to  determine  if  any  particular  receiver  design  is  apt  to 
have  trouble  sooner  than  other  receiver  designs. 

111,7  Task  7:  Derive  the  ECD  to  Secondary  Phase  Relationship 
for  Temporal  Variations 

Utilizing  the  true  variations  in  secondary  phase  and  ECD  that  were 
deduced  after  removing  the  effect  of  the  area  monitor  controlling  the  chain, 
an  effective  mathematical  relation  should  be  derived  from  secondary  phase 
changes  relative  to  ECD  changes.  This  same  relationship  should  be  derived 
theoretically  for  surface  impedances  or  surface  refractive  index  changes 
and  also  for  lapse  rate  changes.  The  theory  suggests  that  a  strong  relation¬ 
ship  should  exist  for  lapse  rate  changes  but  that  a  less  well-defined 
relationship  should  exist  for  surface  parameter  changes. 

The  above  derived  relationship  should  be  used  in  a  computer  simula¬ 
tion  program  designed  to  allow  various  receivers  to  calculate  ECD  to 
demonstrate  a  possible  method  for  receiver  compensation  for  the  temporal 
variation.  During  this  analysis  the  different  ECD  derivation  techniques 
should  be  evaluated  independently  relative  to  the  potential  for  utilizing  a 
relationship  for  deemphasizing  temporal  variations. 

Also,  from  the  observed  temporal  variations,  the  present  opera¬ 
tional  techniques  used  in  the  Loran-C  systems  should  be  evaluated.  Re¬ 
commendations  for  changes  in  the  present  operating  techniques  may  be  forth¬ 
coming  from  this  study.  A  complete  computer  analysis  of  the  temporal 
variations  and  the  use  of  spatially  distributed  area  monitors  should  be  a 
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part  of  the  analysis  program.  The  possible  ramification  on  the  overall 
Loran-C  system  that  would  result  from  holding  the  ECD  constant  at  each 
transmitter  rather  than  constant  at  the  area  monitor  should  be  thoroughly 
discussed. 

A  final  analysis  of  ECD's  and  secondary  phase  corrections  for 
both  spatially  distributed  and  time  distributed  measuring  locations  should 
be  delivered  on  this  project.  As  a  result  of  this  analysis,  definitive 
ratios  of  errors  due  to  spatial  and  temporal  changes  should  be  defined  for 
both  ECD  and  secondary  phase  corrections.  Recommendations  on  best  methods 
for  minimizing  the  effects  of  both  of  these  error  sources  should  be  forth¬ 
coming  from  this  study. 
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IV.  ESTIMATED  TASK  TIME  SCHEDULE  ANALYSIS 


Figure  3  depicts  the  program  in  terms  of  individual  tasks  and 
the  expected  effort  in  man  months  required  to  complete  each  task.  The 
efforts  within  the  tasks  include  both  field  experiments  and  theoretical 
development  as  defined  in  the  task  descriptions  of  Section  III. 

Task  1  is  expected  to  require  a  significant  literature  search 
effort  to  acquire  pertinent  data  for  terrain,  geological  substructure, 
soil,  and  hydrology  in  appropriate  geographical  areas.  Contributing  to 
the  effort  expenditure  are  subtasks  to  place  the  data  in  an  appropriate 
format  to  permit  reliable  modeling  along  the  propagation  paths.  Accordingly, 
this  task  is  foreseen  to  require  the  greatest  expenditure  of  labor  among  all 
the  tasks. 

Task  2  is  anticipated  to  need  a  mixture  of  computer  programming 
expertise  and  novel  theoretical  development  to  convolute  assembled  data 
bases  with  developed  pulse  transient  theory  for  well-founded  ECD  and  secon¬ 
dary  phase  correction  predictions. 

Task  3  will  require  labor  for  van  outfitting  and  modification  and 
approximately  a  nine  man  week  on-site  effort.  The  dotted  line  represents 
expected  U.S.  Coast  Guard  long-term  monitoring  effort. 

Task  4  represents  an  analysis  effort  to  resolve  discrepancies  be¬ 
tween  theoretical  and  experimental  efforts. 

The  purpose  of  the  proposed  efforts  for  Task  5  is  to  bring  into 
concert  the  results  from  the  preceding  tasks  and  those  resulting  from  the 
U.S.  Coast  Guard  dynamic  model  study  to  provide  temporal  propagation  predic¬ 
tion  capability.  To  derive  the  necessary  algorithms  requires  significant 
theoretical  and  computer  discipline  efforts. 

Task  6  is  expected  to  consume  several  man  months  of  labor  in 
field  work  and  measurement  acquisition.  Also,  a  large  effort  is  anticipated 
for  algorithm  modification  to  provide  for  adequate  description  of  the  service 


iqgram  FIGURE  3  EXPECTED  TIME  TASK  SCHEDUL £  ANALYSIS 

TASK 


PROGRAM  TIME  IN  MONTHS 


V,  CONCLUSIONS  AND  RECOfTtNDATIONS 


This  program  is  designed  to  determine  the  predictability  of 
Loran-C  as  an  electromagnetic  ground  wave  propagating  within  geoparametric 
boundary  conditions  and  under  the  influence  of  temporal  variations.  The 
experiment  should  relate  the  sensitivity  of  Loran-C  spatial  variations  to 
surface  impedances  as  derived  from  soil  conductivity,  soil  depth  and 
hydrology,  and  geological  substructure  conductivities.  The  program  should 
demonstrate  and  validate  the  application  of  Pulse  Transient  Solutions  in 
loran  analyses. 

Pursuant  to  the  terms  of  this  contract  and  in  accordance 
with  the  attendant  statement  of  work,  this  report  defines  an  experiment  and 
appropriate  ECD  measurements  to  test  the  pulse  transient  solution.  A  tran¬ 
sient  solution  cannot  be  tested  unless  the  pulse  and  pulse  pharacteristics 
are  considered.  In  this  experiment,  however,  the  varification  of  predicted 
secondary  phase  corrections  are  felt  to  be  more  significant  than  the  ECO 
predictions  and  measurements.  An  Addendum  to  this  report  outlines  a  modest 
program  using  existing  data  that  would  emphasize  the  verification  of  pre¬ 
dicted  secondary  phase  corrections.  This  preliminary  program  would  lay  the 
foundation  for  the  major  program  described  in  this  report.  Additionally, 
this  preliminary  program  will  serve  to  aid  the  Coast  Guard  in  definitivizing 
goals  for  the  overall  program. 

This  program  is  intended  to  resolve  questions  regarding  the 
extent  of  geophysical  data  bases  required  for  high  accuracy  wave  propaga¬ 
tion  characterization.  For  example,  if  the  loran  signal  is  relatively 
insensitive  to  geophysical  conductivities,  the  prediction  problem  is 
ameliorated  by  rather  gross  estimates  of  individual  conductivities.  In 
the  other  extreme,  if  the  loran  is  very  sensitive  to  small  changes  in  con¬ 
ductivities,  prediction  problems  are  increased  by  requirements  to  construct 
adequate  geoparametric  data  bases.  In  the  latter  case,  several  advantages 
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can  be  listed  that  tend  to  offset  the  effort  required  to  assemble  geo¬ 
physical  data  bases.  Firstly,  the  expenditure  of  effort  for  data  collec¬ 
tion  is  expected  to  be  highest  initially  and  decrease  dramatically  as  the 
number  of  geographical  areas  that  are  characterized  increase.  This  expecta¬ 
tion  is  based  upon  the  technique's  universality  wherein  geologic  underlay- 
ment  characteristics  can  be  extrapolated  to  other  geographical  areas  of 
the  earth  surface.  Secondly,  in  the  event  of  high  correlation  between 
loran  signal  propagation  and  geoparametric  variation,  potentials  emerge 
for  other  uses  of  the  loran  signal.  For  example,  the  loran  becomes  a 
potential  tool  for  use  in  geophysical  prospecting  ventures  while  still 
providing  high  accuracy  navigation/positioning  capability.  Accordingly, 
difficulty  in  solving  the  prediction  problem  could  easily  be  offset  by 
interest  and  funding  from  the  geophysical  prospecting  community.  However, 
an  initial  investigation  of  the  type  outlined  in  this  report  is  necessary 
to  establish  the  fundamental  question  of  sensitivity. 

This  experimental  plan  should  develop  suggestions  for  improved 
monitoring  techniques  to  be  used  in  Loran-C  chain  operations.  It  is 
anticipated  that  an  assessment  of  methodologies  for  reduction  of  errors  due 
to  temporal  variations  should  be  forthcomir  from  this  experimental  plan. 
These  methodologies  should  include  suggest!,  .s  for  transfer  of  technology 
necessary  to  incorporate  improvements  into  Loran-C  receivers.  Evaluation 
of  the  possibility  of  incorporating  these  improvements  into  low  cost 
receivers  should  be  discussed. 

As  a  result  of  the  equipment  design  and  measurement  portion  of 
this  experimental  plan,  recommendations  should  be  made  for  equipment  modi¬ 
fication  necessary  for  purposes  of  precise  monitoring  of  Loran-C  signals 
anywhere  within  the  service  area.  Recommendations  for  automating  these 
precise  monitoring  functions  should  also  be  available. 

Finally,  a  definitive  expression  of  the  magnitude  of  the  spatial 
variations  of  both  the  secondary  correction  and  the  ECD  relative  to  the 
temporal  variations  of  these  quantities  for  all  possible  Loran-C  chain 
conditions  should  be  available  from  this  study. 

To  provide  for  timeliness  and  contunity  within  the  program, 
periodic  letter  reports  should  be  provided  which  delineate  program  progress 
and  define  expected  program  activity  for  ensuing  time  periods.  Definitive 
program  reports  should  be  provided  on  a  quarterly  basis  and  the  final  re- 
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port  should  provide  proposed  methods  for  enhancing  calibration  techniques 
through  secondary  phase  and  ECD  measurements  with  considerations  of  both 
spatial  and  temporal  variations.  The  final  report  should  also  include  the 
evaluation  outlined  in  Item  9  of  the  Summary,  In  addition  to  this, 
specific  program  reports  should  be  provided  to  the  Coast  Guard  with  the 
predictions  outlined  in  Item  3  above  and  with  the  results  of  the  comparison 
between  the  predicted  and  measured  ECD  and  secondary  phase  on  all  of  the 

baseline  extensions  as  outlined  in  Item  6. 
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INTERPRETATION  AND  PROJECTED  ANALYSIS  OF  EXISTING 
DATA  APPLIED  TO  BASELINE  PROPAGATION  PATHS 

ADDENDUM  TO  CRPL^  REPORT  78-10 

This  addendum  describes  further  analysis  and  interpretations  si¬ 
milar  to  those  described  in  References  1,  2  and  3.  Interperetations  of  the 
G.E.  Tempo's  worst  case  path  data  and  stability  experiment  results  beyond 
the  preliminary  analysis  mentioned  in  Reference  1  have  been  delayed  until 
the  precise  location  of  the  monitored  points  could  be  obtained.  Table  1 
lists  the  latitude  and  longitude  of  these  points  and  the  distances  to  the 
M,  X  and  Y  transmitters.  As  mentioned  in  Reference  1  and  demonstrated  in 
the  analysis  of  Reference  2,  the  observed  data  would  have  been  greatly  en¬ 
hanced  if  a  TD^  reading  had  been  obtained  at  every  location  visited.  Since 
such  data  is  not  available  the  present  discussions  relate  only  to  TDx  and 
TDy  readings  in  addition  to  the  T0Ay  observations. 

Table  2  gives  the  mean  values  of  the  Y  secondary  phase  corrections 

given  by  G.E.  Tempo,  (Reference  4).  In  addition,  M  and  X  secondary  phase 

corrections  have  been  derived  from  the  mean  TD  and  TD  measurements  also 

y  * 

given  in  Reference  4.  All  of  the  values  listed  in  Table  2  under  the  G.E. 
heading  are  dependent  upon  the  arbitrary  fixing  of  the  Tecopa  TOA  reading  to 
give  a  secondary  phase  correction  of  1224  nanoseconds.  All  of  the  values 
listed  in  Table  2  under  the  CRPL^  heading  are  based  on  fixing  the  value  to 
give  a  secondary  phase  correction  of  1624  nanoseconds.  This  value  is  much 
closer  to  the  value  derived  in  Reference  5. 

Table  3  lists  effective  impedances  to  each  transmitter  for  each 
location  using  the  two  sets  of  data  given  in  Table  2  and  techniques  described 
In  References  1,  2  and  3.  We  believe  the  results  listed  under  CRPL-j  in  Table 
3  are  in  better  agreement  with  the  Coast  Guard  and  NOAA  measurements  pre¬ 
sented  in  Reference  2.  For  example,  the  last  4  entries  of  Table  5,  (Reference 
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2),  all  had  X  paths  propagating  over  the  San  Joaquin  Valley.  Also,  NOAA 
ship  measureme.it  points  1390,  1411,  1430  and  1716  shown  in  Figure  1  and 
Table  8  of  Reference  2  had  paths  propagating  over  the  Valley.  The  X  paths 
to  Friant,  Merced,  Crows  Landing  and  to  a  lesser  extent  Livermore,  pro¬ 
pagated  over  the  San  Joaquin  Valley.  It  would  seem  reasonable  to  expect 
impedances  from  .0148  to  .0288  for  these  paths,  but  unreasonably  low  for 
impedances  from  .0043  to  .0165  for  the  paths. 

Figure  1  illustrates  the  geometrical  relations  between  the  Loran- 
C  transmitters  and  the  Worst  Case  Path  measurement  locations  occupied  by 
G.E.  Tempo.  Figure  2  shows  the  terrain  plot  from  Searchlight,  Nevada  to 
Fort  Cronkhite,  California,  along  with  the  approximate  locations  and  deduced 
impedance  values  for  the  intermediate  points  along  this  path.  It  is  inter¬ 
esting  to  note  that  the  impedance  values  consistently  increase  from  Tecopa 
to  Friant  where  the  path  is  very  irregular  and  the  values  tend  to  decrease 
after  Friant  as  the  path  crosses  the  San  Joaquin  Valley.  However,  the  slightly 
rugged  terrain  between  Crows  Landing  and  Livermore  apparently  increases  the 
impedance  again  slightly  at  Livermore.  Figures  3  and  4,  terrain  plots  from 
Searchlight,  Nevada  to  Death  Valley,  Calif.,  and  Searchlight,  Nevada  to 
Friant,  Calif.,  respectively  demonstrate  terrain  differences  from  the  path 
in  Figure  2  due  to  the  fact  that  intermediate  points  were  not  precisely  on 
the  geodesic  path  between  Searchlight  and  Fort  Cronkhite. 

A  complete  analysis  over  the  total  path  in  Figure  2  and  the  actual 
intermediate  paths  to  each  point  similar  to  the  Searchlight  to  Tecopa  analysis 
given  in  Reference  5  would  demonstrate  several  factors.  One,  it  would  demon¬ 
strate  the  potential  advantage  of  a  complete  calibration  by  prediction  enter¬ 
prise.  Two,  it  would  demonstrate  the  sensitivity  of  the  Loran-C  changes 
with  slight  changes  in  path  geometry.  Three,  it  would  tend  to  verify  the 
deduced  impedance  values  for  the  Y  path  and  consequently  tie  down  the  im¬ 
pedance  values  for  both  the  M  and  X  paths  that  are  discussed  in  the  next 
section. 

In  Figures  5  through  12  the  terrain  from  each  point  on  the  Worst 
Case  Path  to  the  Master  is  plotted.  Also  given  on  each  figure  is  the  deduced 
impedance  values.  As  would  be  expected  from  looking  at  Figure  1,  the  terrain 
is  extremely  variable  from  path  to  path  and  only  very  close  to  the  transmitter 
is  there  any  similarity  at  all.  The  deduced  impedance  values  appear  to  be 
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reasonably  related  to  the  terrain  variations  with  the  exception  of  the 
impedances  to  Death  Valley  and  Darwin.  If  the  impedances  are  primarily 
determined  by  terrain,  it  would  seem  that  these  two  impedances  should  be 
interchanged.  A  complete  integral  equation  analysis  similar  to  that  in¬ 
dicated  in  Reference  5  might  resolve  this  question.  In  a  complete  prediction 
analysis  the  differences  in  soil  and  sub  soil  parameters  between  the  Master 
to  Death  Valley  and  the  Master  to  Darwin  paths  could  be  studied  closely  to 
determine  the  reason  why  the  impedance  was  higher  to  Darwin  than  it  was  to 
Death  Valley. 

Figures  13  through  20  show  the  terrain  plots  from  the  X  transmitter 
to  each  point  and  also  the  deduced  impedance  values  for  these  paths.  Again 
only  the  portions  of  the  paths  close  to  the  transmitter  station  are  similar 
as  would  be  expected.  The  impedances  for  the  first  three  points  are  no 
doubt  primarily  determined  by  the  rough  terrain.  The  specific  values  for 
each  of  these  three  points  may  also  be  somewhat  affected  by  ground  conducti¬ 
vities  along  the  paths.  The  impedance  values  for  the  next  four  points  are 
no  doubt  primarily  influenced  by  the  San  Joaquin  Valley  and  the  lack  of 
terrain  features.  The  deduced  impedance  from  the  X  transmitter  to  Fort 
Cronkhite  is  surprisingly  high.  A  study  using  the  integral  equation  techniques 
of  Reference  5  could  help  to  understand  these  anomolous  points. 

In  addition  to  the  above  analysis  a  similar  analysis  is  presented 

for  the  baseline  paths  and  the  near  baseline  paths  to  Silyer  Springs, 

Nevada,  Jean,  Nevada,  and  Arbuckle,  California  (see  Figure  21).  The  impedances 

deduced  for  the  X  Secondary  to  the  Master  (Figure  22)  and  the  Master  to  Y 

Secondary  (Figure  23)  from  baseline  extension  measurements  given  in  Reference 

2  was  A  =  .049  and  A  =  .0423  respectively.  The  data  given  in  Table  4  A 
x  y 

is  the  deduced  secondary  phase  corrections  for  each  path  to  Silver  Springs, 

Jean  and  Arbuckle.  This  data  was  deduced  from  the  distances  given  in  Table  1 
and  the  time  difference  averages  given  in  Reference  6.  Table  4  B  gives  the 
effective  impedances  to  each  transmitter  associated  with  the  secondary  phase 
corrections  given  in  Table  4  A.  The  three  impedances  underlined  in  the 
table  represent  the  short  path  impedances  for  paths  shown  in  Figures  24,  25 
and  26  and  cannot  be  independently  checked.  However,  these  impedances  may 
be  indicative  of  surface  conductivities  in  the  vicinity  of  the  transmitting 
stations.  The  paths  marked  path  la  and  b  (Figures  27  and  28),  path  2a  and  b 
(Figures  29  and  30),  and  path  3a  and  b  (Figures  31  and  32)  are  nearly 


reciprocal.  Also,  paths  la  and  b  are  close  to  the  Master  -  X  Secondary 
(Figure  22)  and  paths  2a  and  b  are  close  to  the  Master  -  Y  Secondary  (Figure 
23).  Paths  la  and  lb  at  .0490  and  .0495  agree  well  with  the  A^  of  .049. 

Paths  3a  and  b  at  .0458  and  .0450  agree  well  with  each  other.  Also,  path  2b 
at  .0420  agrees  well  with  Ay  =  .0423.  However,  path  2a  at  .0476  does  not 
agree  well  with  =  .0423.  Jean,  Nevada  is  very  close  to  the  baseline  from 
Master  to  Searchlight,  and  the  terrain  for  that  path  is  presented  in  Figure 
23.  It  can  be  seen  that  Jean,  Nevada  falls  at  a  point  on  this  path  directly 
following  the  highest  terrain  on  the  entire  path.  From  previous  integral 
equation  runs,  see  Reference  5,  for  example,  the  maximum  impedance  values 
always  follow  the  greater  terrain  peaks.  Therefore  it  could  be  anticipated 
from  this  terrain  plot  alone  that  the  impedance  from  the  Master  to  Jean, 

Nevada  should  be  higher  than  that  from  the  Master  to  Searchlight,  Nevada. 

Figure  33  is  the  baseline  terrain  between  the  X  and  Y  transmitters. 
This  path  is  shown  in  comparison  to  paths  3a  and  3b  that  is  from  Y  to  Ar buckle 
and  from  X  to  Jean.  An  estimated  impedance  of  .0454  was  assigned  to  this  path 
only  as  an  average  of  the  deduced  impedances  for  paths  3a  and  3b.  This  im¬ 
pedance  could  be  independently  determined  by  making  integral  equation  calcula¬ 
tions  similar  to  those  given  in  Reference  5.  This  impedance  could  also  be 
independently  measured  by  making  X  -  Y  and  Y  -  X  baseline  extension  measure¬ 
ments  and  analyzing  the  data  as  was  done  for  the  M  -  X  and  M  -  Y  baseline 
extension  measurements  as  reported  in  Reference  2. 

Ground  structure  analysis  and  integral  equation  computations  to 
all  of  these  stations  would  demonstrate  more  fully  the  potential  of  Loran-C 
navigation.  Furthermore,  computations  using  the  full  wave  transient  analysis 
algorithm  for  all  of  these  paths  could  help  the  Coast  Guard  determine  the 
true  extent  of  the  envelope  to  cycle  problem  over  irregular  nonhomogeneous 
terrain. 
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LOCATIONS,  DISTANCE  TO  TRANSMITTERS  AND  PREDICTED  TDs  WITHOUT  SECONDARY 
PHASE  CORRECTION  FOR  POINTS  MONITORED  BY  G.E.  TEMPO 


TABLE  2 

SECONDARY  CORRECTIONS  TO  TRANSMITTERS  FOR  WORST  CASE  PATH  POINTS 


To  Y  Secondary 

GE  CRPLi 

To  Master 

GE  CRPLi 

To  X  Secondary 
GE  CRPLi 

Tecopa 

1.224 

1.624 

2.39 

2.79 

4.33 

4.73 

Death  Valley 

1.829 

2.229 

2.55 

2.95 

4.66 

5.06 

Darwin 

2.463 

2.863 

2.85 

3.25 

3.11 

3.51 

Frlant 

3.960 

4.360 

3.16 

3.56 

.93 

1.33 

Merced 

3.475 

3.875 

2.28 

2.68 

.62 

1.02 

Crows  Landing 

3.406 

3.806 

2.01 

2.41 

.38 

.78 

Livermore 

4.178 

4.578 

2.04 

2.44 

00 

VO 

• 

1.08 

Ft.  Cronkhlte 

4.080 

4.480 

2.01 

2.41 

1.05 

1.45 

TABLE  3 

CALCULATED  EFFECTIVE  IMPEDANCES  FOR  WORST  CASE  PATH  POINTS 


Y  Secondary 

GE  CRPLi 

Master 

GE  CRPLi 

X  Secondary 
GE  CRPLi 

Tecopa 

.0344 

.0468 

.0258 

.0317 

.0419 

.0468 

Death  Valley 

.0413 

.0506 

.0318 

.0382 

.0504 

.0554 

Darwin 

.0446 

.0516 

.0415 

.0483 

.0349 

.0407 

Frlant 

.0484 

.0537 

.0542 

.0591 

.0105 

.0183 

Merced 

.0378 

.0433 

.0394 

.0473 

.0069 

.0158 

Crows  Landing 

.0323 

.0376 

.0324 

.0404 

.0043 

.0148 

Livermore 

.0385 

.0433 

.0306 

.0382 

.0165 

.0288 

Ft.  Cronkhlte 

.0331 

.0378 

.0267 

.0337 

.0342 

.0482 
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TABLE  4  A 

SECONDARY  PHASE  CORRECTIONS  FOR  STABILITY 
EXPERIMENT  MONITORING  POINTS 


I 


Master 

X  Secondary 

Y  Secondary 

Silver  Springs 

.45  y  sec. 

2.76  y  sec. 

4.09  y  sec. 

Arbuckle 

2.70  y  sec. 

.97  y  sec. 

5.26  y  sec. 

Jean 

4.09  y  sec. 

4.95  y  sec. 

.83  y  sec. 

TABLE  4  B 

EFFECTIVE  IMPEDANCES  FOR  STABILITY 
EXPERIMENT  MONITORING  POINTS 

Master  X  Secondary 

Silver  Springs  .0215  .0495  (path 

Arbuckle  .0490  (path  la)  .0459 

.0476  (path  2a)  .0450  (path 


Y  Secondary 
lb)  .0420  (path  2b) 
.0458  (path  3a) 
3b)  .0329 


Oean 
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FIGURE  2.  Terrain  Path  from  Searchlight,  Nev.  to  Fort  Cronkhite,  Calif,  with  approximate  locations 
of  intermediate  points  and  derived  impedances  for  these  points  from  the  CRPLj  COMPUTER I Z1 
Continental  United  States  terrain  data  base. 


TRANSMITTER  TO  DEATH  VALLEY 
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FIGURE  3.  Precise  Terrain  Path  from  Searchlight.,  Nevada  to  Death  Valley,  Calif.,  i 

DERIVED  IMPEDANCE  FOR  THIS  PATH  FROM  THE  CRPLj  COMPUTERIZED  CONTINENTAL 

United  States  terrain  data  base. 


Y  TRANSMITTER  TO  FRIANT 


Precise  Terrain  Path  from  Searchlight,  Nevada  to  Friant,  Calif.,  and  derived 

IMPEDANCE  FOR  THIS  PATH  FROM  CRPLj  COMPUTERIZED  CONTINENTAL  UNITED  STATES 
TERRAIN  DATA  BASE. 
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FIGURE  6.  Terrain  Path  from  Fallon 

FOR  THIS  PATH  FROM  CRPL 


FIGURE  7,  Terrain  Path  from  Fallon,  Nevada  to  Darwin,  Calif,  and  derived  impedance  for 

THIS  PATH  FROM  CF^PLj  COMPUTERIZED  CONTINENTAL  UNITED  STATES  TERRAIN  DATA  BASE 


i zed  Continental 


TO  MERCED 


Lon,  Nevada  to  Merced,  Calif,  and  Derived  Impedance  for 
Computerized  Continental  United  States  terrain  data  base 


FIGURE  ID.  Terrain  Path  from  Fallon,,  Nevada  to  Crows  Landing,  Calif,  and  Derived 
Impedance  for  this  path  from  CRPLj  Computerized  Continental  United 
States  Terrain  Data  Base. 
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FIGURE  11.  Terrain  Path  from  Fallon,  Nevada  to  Livermore,  Calif,  and  Derived 
Impedance  for  this  path  from  CRPLi  Computerized  Continental  United 
States  Terrain  Data  Base. 


X  TRANSMITTER  TO  TECOPA 
A  =  .0468 


FIGURE  13.  Terrain  Path  from  Middletovn,  Calif,  to  Tecopa,  Calif,  and  Derived 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental  United 
States  Terrain  Data  Base. 


X  TRANSMITTER  TO  DEATH  VALLEY 
A  =  .0554 


Terrain  Path  from  Middletown.,  Calif,  to  Death  Valley,  Calif,  and  Derived 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental  United  States 
Data  Base. 


X  TRANSMITTER  TO  DARWIN 


FIGURE  15.  Terrain  Path  from  Middletown,  Calif,  to  Darwin,  Calif,  and  Derived 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental  United 
States  Terrain  Data  Base. 
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FIGURE  16.  Terrain  Path  from  Middletown,  Calif,  to  Friant,  Calif,  and  Derived 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental  United 
States  Terrain  Data  Base. 


TRANSMITTER  TO  MERCED 
A  =  .0158 


Terrain  Path  from  Middletown  Calif,  to  Forced,  Calif,  and  Derived 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental  United 
States  Terrain  Data  Base. 


Path  from  Middletown,  Calif,  to  Fort  Cronkhite,  Calif,  and 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental 
tates  Terrain  Data  Base. 


LONGITUDE 

Transmitter  Paths  to  Silver  Springs,,  Nevada;  Jean,  Nevada  and 
Arbuckle,  California. 
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Terrain  Path  from  Middletown,  Calif,  to  Arbuckle,  Calif,  and  Derived 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental  United 
States  Terrain  Data  Base. 


Figure  27.  Path  1a  Terrain  Path  from  Fallon.,  Nevada  to  Arbuckle,  Calif,  and 

Derived  Impedance  for  this  Path  from  CRPLj  Computerized  Continental 
United  States  Terrain  Data  Base. 
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FIGURE  28.  Path  Ib  Terrain  Path  from  Middletown,  Calif,  to  Silver  Springs,  Nevada 
and  Deduced  Impedance  for  this  Path  from  CRPLj  Computerized  Continental 
United  States  Terrain  Data  Base. 


ASTER  TO  JEAN 
A  *  .0476 


FIGURE  29.  Path  2a  Terrain  Path  from  Fallon.,  Nevada  to  Jean,  Nevada  and  Deduced 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental  United 
States  Terrain  Data  Base. 
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United  States  Terrain  Data  Base. 


PATH  3a  Y  TRANSMITTER  TO  ARBUCKLE 
A  -  .0458 


Terrain  Path  from  Searchlight,  Nevada  to  Arbuckle,  Calif,  ai 
Impedance  for  this  Path  from  CRPLj  Computerized  Continental 
tates  Terrain  Data  Base. 


PATH  3b  X  TRANSMITTER  TO  JEAN 
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FIGURE  32.  Path  3b  Terrain  Path  from  Middletow,  Calif,  to  Jean,  Nevada  and 

Deduced  Impedance  for  this  Path  from  CRPLx  Computerized  Continental 
United  States  Terrain  Data  Base, 
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FIGURE  33.  Terrain  Path  from  Middletown,  Calif,  to  Searchlight,  Nevada  and  an 
Estimated  Impedance  from  Paths  3a  and  3b  for  this  path  from  CRPLj 
Computerized  Continental  United  States  Terrain  Data  Base. 


